- 1 downto 2 do:=0; c:=0; b:=0; f:=0; a:=0; d:=0 ; R:=0; h [i, j] <> 0 thenh [i, j - 1]=0 then g:=0g:=(dx / dy) * p [i, j] * (kh ^ [i, j - 1] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i - 1, j]=0 then c:=0c:= (dy / dx) * p [i, j] * (kh ^ [i - 1, j] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i +1, j]=0 then b:=0b:=(dy / dx) * p [i, j] * (kh ^ [i +1, j] + kh ^ [i, j]) / mu/z (p [i, j] / 1e6) / 2; h [i, j +1]=0 then f:=0f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];:=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, g1 ^ [i - 1, j], f1 ^ [i, j +1], al [k], V ^ [i - 1, j], V ^ [i, j +1], ggg, fff, v1); ^ [i, j]:=v1 ; ^ [i, j]:=ggg; ^ [i, j]:=fff;;;:=true; i:=N - 1 downto 2 doj:=2 to M - 1 doh [i, j] < ;> 0 then ^ [i, j]:=V ^ [i, j]-g1 ^ [i, j] * dp ^ [i, j - 1]-f1 ^ [i, j] * dp ^ [ i +1, j]; 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 g1 ^ [i, j]:=0; f1 ^ [i, j]:=0 end;
Продовження додатка А
enD
eLSe
begiN
for i:=2 to N - 1 doj:=2 to M - 1 do:=0; c:=0; b:=0; f:=0; a:=0; d:=0; R:=0; h [i, j] <> 0 thenh [i, j - 1]=0 then g:=0g:=(dx / dy) * p [i, j] * (kh ^ [ i, j - 1] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i - 1, j]=0 then c:=0c:=( dy / dx) * p [i, j] * (kh ^ [i - 1, j] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [ i +1, j]=0 then b:=0b:=(dy / dx) * p [i, j] * (kh ^ [i +1, j] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i, j +1]=0 then f:=0f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];:=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;_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));;VIIa(giv1:m1;var gi1: m1; dt, t: real; var PcpVIIa,p5_1,p8_1,p10_1,p16_1,p75_1:real;_1,q5_2,q8_1,q8_2,q10_1,q10_2,q16_1,q16_2,q75_1,q75_2:real);=16;{strok}=34;{stolbcov}
Ppl=7.6e6; {Pa}=1.877e - 12; {mkm2, IV block}=2.062e - 12; {mkm2, III block}=0.3;=0.78;=0.76;=0.011e - 3; {Pa * sec}
dx=100; {m}=100; {m}=100000;=86400;=real;=array [1 .. 10] of tipR;=array [1 .. N, 1 .. M] of tipR;
dmas=^ mas2;, dat1: text;
al:mas1;,j,count,k,cell:integer;,p,h:mas2;,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;
Продовження додатка А
b1, g1, f1, kh, V, dp, pp: dmas;
flag: boolean; alf (var...
