if K[I,J]<=0.6 then
   Begin
               A1:=-0.100021*K[I,J]+8.887012;
               A2:=0.279999*K[I,J]+8.672999;
               B1:=859.999756*K[I,J]+1884.000122;
               B2:=1159.999756*K[I,J]+1762.000122;
               K1:=K[I,J]-5.000001e-2;
               K2:=K[I,J];
               P6:=0; goto C1;
   End;

 if K[I,J]<=0.65 then
   Begin
               A1:=0.279999*K[I,J]+8.658999;
               A2:=0.240002*K[I,J]+8.696997;
               B1:=1160.001221*K[I,J]+1703.999268;
               B2:=1500.001343*K[I,J]+1557.999146;
               K1:=1.000001*K[I,J]-0.050000727;
               K2:=K[I,J];
               P6:=0; goto C1;
   End;

if K[I,J]<=0.7 then
    Begin
              A1:=0.240002*K[I,J]+8.684998;
              A2:=3.579997*K[I,J]+6.526001;
              B1:=1499.999634*K[I,J]+1483.000244;
              B2:=3099.999268*K[I,J]+518.000549;
              K1:=0.999999*K[I,J]-0.0499918;
              K2:=K[I,J];
              P6:=0;
              goto C1;
    End;

if K[I,J]<=0.75 then
    Begin
              A1:=3.579997*K[I,J]+6.347001;
              A2:=4.005432e-002*K[I,J]+9.003961;
              B1:=3099.999268*K[I,J]+363.000549;
              B2:=2439.999512*K[I,J]+980.000427;
              K1:=K[I,J]-0.05000001;
              K2:=K[I,J];
              P6:=0;goto C1;
    End;

if K[I,J]<=0.8 then
   Begin
              A1:=4.005432e-002*K[I,J]+9.0019579;
              A2:=5.179976*K[I,J]+5.149019;
              B1:=2439.999512*K[I,J]+858.000427;
              B2:=4599.999023*K[I,J]-639.999207;
              K1:=K[I,J]-0.0500002;
              K2:=K[I,J];
              P6:=0; goto C1;
   End;

if K[I,J]<=0.85 then
   begin
               A1:=5.179976*K[I,J]+4.890020;
               A2:=-1.079979*K[I,J]+10.156983;
               B1:=4599.999023*K[I,J]-869.999206;
               B2:=2699.9992068*K[I,J]+880.000488;
               K1:=K[I,J]-0.05002;
               K2:=K[I,J];
               P6:=0; goto C1;
   End;

if K[I,J]<=0.9 then
  Begin
               A1:=-1.07998*K[I,J]+10.21983;
               A2:=-0.319976*K[I,J]+8.967020;
               A3:=7.751;
               A4:=7.897;
               B1:=2700.002686*K[I,J]+744.997742;
               B2:=4300.003906*K[I,J]-480.003601;
               B3:=3742;
               B4:=3685;
               K1:=K[I,J]-0.05000102;
               K2:=K[I,J];
               K3:=0.85;
               K4:=0.9; goto Ce;
  End;

if K[I,J]<=0.93 then
  Begin
               A1:=0.533293*K[I,J]+8.759037;
               A2:=17.833309*K[I,J]-6.794978;
               A3:=4.866659*K[I,J]+3.371007;
               A4:=9.099998*K[I,J]-0.292999;
               B1:=499.99948*K[I,J]+2925.0004;
               B2:=16599.982422*K[I,J]-11549.984375;
               B3:=-1899.997925*K[I,J]+5451.998047;
               B4:=-966.665588*K[I,J]+4554.999023;
               K1:=1.666663*K[I,J]-0.6499999;
               K2:=1.666665*K[I,J]-0.599999;
               K3:=1.666663*K[I,J]-0.649997;
               K4:=K[I,J]; goto Ce;
  End;

if K[I,J]<=0.95 then
   Begin
               A1:=9.2550;
               A2:=9.789999;
               A3:=-36.350025*K[I,J]+41.702523;
               A4:=7.300003*K[I,J]+1.380998;
               B1:=3390;
               B2:=3888;
               B3:=-1450.001587*K[I,J]+5033.501465;
               B4:=1550.001587*K[I,J]+2214.498535;
               K1:=0.89999;
               K2:=0.94999;
               K3:=1.5*K[I,J]-0.495003;
               K4:=K[I,J];goto Ce;
  End;

if K[I,J]<=0.983 then
   Begin
               A1:=16.212120*K[I,J]-6.146514;
               A2:=-0.303037*K[I,J]+10.077886;
               A3:=34.727276*K[I,J]-25.820913;
               A4:=4.666677*K[I,J]+3.882658;
               B1:=15090.913086*K[I,J]-10946.367188;
               B2:=9787.880859*K[I,J]-5410.485816;
               B3:=939.394043*K[I,J]+2763.575684;
               B4:=-2848.485352*K[I,J]+6393.061035;
               K1:=1.515152*K[I,J]-0.539395;
               K2:=K[I,J];
               K3:=0.60606*K[I,J]+0.354243;
               K4:=0.909092*K[I,J]+0.08636254;
   End;


Ce:
 A6:=A3+(A4-A3)/(K4-K3)*(K[I,J]-K3);
 B6:=B3+(B4-B3)/(K4-K3)*(K[I,J]-K3);
 P6:=133.3*exp(ln(10)*(A6-(B6/Q[I,J])));
C1:
 A0:=A1+(A2-A1)/(K2-K1)*(K[I,J]-K1);
 B0:=B1+(B2-B1)/(K2-K1)*(K[I,J]-K1);
 P7:=133.3*exp(ln(10)*(A0-(B0/Q[I,J])));
 P9:=P7-P6;
 P8:=P0*Z[I,J]/(18/29+Z[I,J]);