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

end;
procedure valvula;
var
xl,ax,
vl,hl : double;
j : integer;
begin
if t > tf
then
qvs[2] := 0.0e0
else
qvs[2] := (tf-t)*q0/tf;
xl := lt[2] - a[2]*dt;
j := trunc(xl/dx[2]);
ax := (xl - j*dx[2])/dx[2];
vl := v[2,1]^[j] + (v[2,1]^[j+1] - v[2,1]^[j])*ax;
hl := h[2,1]^[j] + (h[2,1]^[j+1] - h[2,1]^[j])*ax;
v[2,2]^[nx[2]] := qvs[2]/at[2];
h[2,2]^[nx[2]] := hl + a[2]*((vl - v[2,2]^[nx[2]])/g -
dt*a1[2]*vl*abs(vl))
end;
begin
for i := 1 to 2 do
for j := 1 to 2 do
begin
new(v[i,j]);
new(h[i,j]);
end;
inic_dad;
assign (arq,'lax_f02.txt');
rewrite(arq);
writeln(arq,0.0e0:10:5,' ',h[1,2]^[0],' ',hch[2],
139