UNIT REWL;
INTERFACE
uses rewob,pechob,rewlibl;
procedure rew;
IMPLEMENTATION
procedure rew;
label 2,6,8,11,13,15,16,17,22,23,26,29,30,31,32,33,34,35;
VAR e1,e,e1c,ec,de,ft,fto,zbn1,mcm,mm,s,rcm,
q1,dqk,dqp,mq,ta,k2,rcp,moct1,k1,tb,d15,tn,tn1,ko1,n,
tan,qcm,vn,vbn,zn,tkp,pkp,r15,tpr,ppr,al1,r4t,xm,xr,ek,al,be,
sf,dn2,de1,de2,qbp,vcm,rocm,stp,sq,sn1,sn2,qn,zb1,m1:real;
k3,m,x:array[1..ifm]of real;b:array[1..5]of real;
it,mo,i,i8,losn:integer;
TS:ARRAY[1..7] OF STRING[98];
procedure pp6 ;
VAR s,m1:real;
i:integer;
begin
ta:=0;s:=0;m1:=0.;
for i:=1 to if1 do
begin s:=s+x[i]*m[i]/r420[i];
m1:=m1+x[i]*m[i];
ta:=ta+x[i]*toi[i];
end;
rcp:=m1/s;
rcm:=0.994*(rcp)+0.0093;
k2:=1.216*pow(ta+273.16,0.333333)/rcm;
k2m[jb]:=k2;
end; (* < ------------------------ end of pp6 -----------------------------> *)
BEGIN
nnm:=1;
{writeln(outfl,'fb=',fb:8:1,' dbn=',dbn:8:2,' zbn=',zbn:8:2,' moct=',moct:9:1,' t=',t:8:2,' ',p:8:2);}
rpri3 (t,p);
{for i:=1 to if1 do writeln(outfl,f[i]:12:4,' ',toi[i]:12:2,' ',r420[i]:12:4);}
rpri3 (t,p);
it:=0;xb:=0;
for i:=iy+1 to n2 do xb:=xb+f[i];
for i:=1 to if1 do f1[i]:=f[i];
tn:=moct;
ko1:=0.1;
mo:=0;
for i:=1 to 18 do
begin d15:=0.994*r420[if1]+0.0093;
tb:=tn+273.16;
k1:=1.216*pow(tb,0.333333)/d15;
moct1:=160-5*k1-0.075*tn+0.000156*k1*tn*tn;
s:=moct1/moct;
if(abs(s-1)<1e-04) then goto 29;
if(mo=0) then goto 30;
if(abs(tn1-tn)<1) then goto 30;
ko1:=(ln(tn1)-ln(tn))/(ln(s)-ln(mm));
30: mm:=s; tn1:=tn;
tn:=tn1*exp(ko1*ln(s));
mo:=1;
end;
29: toi[if1]:=tn;s:=0;
for i:=1 to if1 do
begin if (i<=iy) then begin m[i]:=my[i]; goto 2; end;
d15:=0.994*r420[i]+0.0093;
tb:=toi[i]+273.16;
k3[i]:=1.216*pow(tb,0.333333)/d15;
m[i]:=160-5*k3[i]-0.075*toi[i]+0.000156*k3[i]*toi[i]*toi[i];
2: rmy[i]:=m[i];
s:=s+f[i]/m[i];
end;
for i:=1 to if1 do f[i]:=f[i]/(m[i]*s);
mcm:=100/s;
if (zbn>0.001) then begin zbn1:=100*zbn*mcm/(fb*18);end
else zbn1:=0;
if (nnm=1) then begin
ft:=37*t/(t+273.16);
for i:=1 to if1 do
begin fto:=37*toi[i]/(toi[i]+273.16);
k[i]:=exp(0.7944*(10.4*(ft-fto)/(31.6-fto)+0.1))/p;
if (k[i]<1e-17) then k[i]:=1e-10;
end; end; |