avaliações numéricas de chaminés de equilíbrio - ppgerha ...

egin
for i := 1 to 2 do
for j := 1 to nx[i]-1 do
begin
vm := 0.5e0*(v[i,1]^[j-1] + v[i,1]^[j+1]);
hm := 0.5e0*(h[i,1]^[j-1] + h[i,1]^[j+1]);
dvx := a3[i]*(v[i,1]^[j+1] - v[i,1]^[j-1]);
dhx := a3[i]*(h[i,1]^[j+1] - h[i,1]^[j-1]);
v[i,2]^[j] := vm - dt*(g*dhx + vm*dvx +
a1[i]*vm*abs(vm));
h[i,2]^[j] := hm - dt*(a2[i]*dvx + vm*dhx +
vm*st[i]);
end;
end;
procedure reservatorio;
begin
v[1,2]^[0] := v[1,2]^[1];
h[1,2]^[0] := zr - b1*v[1,2]^[0]*abs(v[1,2]^[0]);
end;
procedure chamine;
var
aux : double;
begin
v[1,2]^[nx[1]] := v[1,2]^[nx[1]-1];
v[2,2]^[0] := v[2,2]^[1];
qch[2] := v[1,2]^[nx[1]]*at[1] - v[2,2]^[0]*at[2];
hch[2] := hch[1] + b2*(qch[1] + qch[2]);
aux := b3*qch[2]*abs(qch[2]);
h[1,2]^[nx[1]] := hch[2] + aux;
h[2,2]^[0] := hch[2] - aux;
138