]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);:=2*pat*qm[i,j]*297/293-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);:=d-(g*pp^[i,j- 1] + c * pp ^ [i - 1, j] + a * pp ^ [i, j] + b * pp ^ [i +1, j] + f * pp ^ [i, j +1]) ; odd (count) then k:=k +1; k> 10 then k:=1; (i, j, g, c, a, b, f, R, b1 ^ [i, j - 1], f1 ^ [i - 1, j], al [k], V ^ [i - 1, j], V ^ [i, j - 1], bbb, fff, v1); ^ [i, j]:=v1 ; ^ [i, j]:=bbb; ^ [i, j]:=fff;;;:=true; i:=N - 1 downto 2 doj:=M - 1 downto 2 doh [i, j] < ;> 0 then ^ [i, j]:=V ^ [i, j]-b1 ^ [i, j] * dp ^ [i +1, j]-f1 ^ [i, j] * dp ^ [ i, j +1]; abs (dp ^ [i, j])> 1e - 4 then flag:=false; ^ [i, j]:=pp ^ [i, j] + dp ^ [i, j ];;
{} for i:=1 to N doj:=1 to M do begin b1 ^ [i, j]:=0; f1 ^ [i, j]:=0 end;; flag;:=0; Pcp:=0; Hcp:=0; i:=1 to N doj:=1 to M do [i, j]:=pp ^ [i, j];:=ZapT +293 / 297 * h [i, j] * dx * dy * kp * al1 * p [i, j] / pat / z (p [i, j] / 1e6) * z (pat/1e6);
Продовження додатка А
Qdob:=Qdob + qm [i, j] * dt;
if h [i, j] <> 0 then begin Pcp:=Pcp + pp ^ [i, j] * h [i, j]; Hcp:=Hcp + h [i, j] end; ;;
{if (ZapN-ZapT-Qdob) / Qdob * 100> delta then delta:=(ZapN-ZapT-Qdob)/Qdob*100;}:=Pcp/Hcp;:=Pcp;_2:=pp^[10,13];_2:=pp^[10,26];_2:=pp^[6,15];_2:=pp^[7,29];_2:=pp^[9,5];:=pp^[14,23];:=pp^[15,11];(kh,M*N*SizeOf(tipR));(b1,M*N*SizeOf(tipR));(pp,M*N*SizeOf(tipR));(dp,M*N*SizeOf(tipR));(V,M*N*SizeOf(tipR));(g1,M*N*SizeOf(tipR));(f1,M*N*SizeOf(tipR));;soprVII(gi:m;dt,t:real; var u5_2, u8_2, u10_2, u16_2, u75_2, u11, u15: real;
e5_1, e5_2, e8_1, e8_2, e10_1, e10_2, e16_1, e16_2, e75_1, e75_2, e11, e15: real);
const=19; {strok}=34; {stolbcov}
Ppl=7.6e6; {Pa} {!!!!!!}=0.661e - 12; {mkm2, IV block}=0.44e - 12; { mkm2, III block}=0.32;=0.77;=0.76;=0.011e - 3; {Pa * sec}
dx=100; {m}=100; {m}=100000;=86400;=real;=array [1 .. 20] of tipR;=array [1 .. N, 1 .. M] of tipR;
dmas=^ mas2;, dat1: text;
al:mas1;,j,count,k,cell:integer;,p,u:mas2;,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;,g1,f1,kh,V,dp,pp:dmas;:boolean;alf(var al: mas1);
var, j: integer;, par1, par11, par2: tipR;
begin:=100; i:=1 to N - 1 doj:=1 to M - 1 doh [i, j] <> 0 then (h [i, j +1] <> 0 ) and (h [i +1, j] <> 0) then:=dx * dx * (kh ^ [i, j +1] + kh ^ [i, j]) / (kh ^ [i + 1, j] + kh ^ [i, j]) / dy / dy;
par1:=pi * pi / (2 * N * N * (1 + par));:=pi * pi / (2 * M * M * (1 +1 / par));
if par1
{5} [10,13]:=e5_2 {/ 86.4};
{75}
qm [9,5]:=e75_2 {/ 86.4};
Продовження додатка А
{10}
qm [6,15]:=e10_2 {/ 86.4};
{11}
qm [15,11]:=e11 {/ 86.4};
{8}
qm [10,26]:=e8_2 {/ 86.4};
{16}
qm [7,29]:=e16_2 {/ 86.4};
{15}[14,23]:=e15{/86.4};;(kh,M*N*SizeOf(tipR));(b1,M*N*SizeOf(tipR));(pp,M*N*SizeOf(tipR));(dp,M*N*SizeOf(tipR));(V,M*N*SizeOf(tipR));(g1,M*N*SizeOf(tipR));(f1,M*N*SizeOf(tipR));i:=1 to N doj:=1 to M doh [i, j] <> 0 thenj <= 21 then kh^[i,j]:=k2*h[i,j]kh^[i,j]:=k1*h[i,j]kh^[i,j]:=0;[i,j]:=gi[i,j];[i,j]:=0;^[i,j]:=u[i,j];[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;;;(al);:=0;delta:=0;:=0;:=0;;:=count+1;odd(count) then {obratnoe yporyadochivanie} i:=2 to N - 1 doj:=M ...
