Design og modellering af metanolanlæg til VEnzin-visionen Bilag

methanol.for
d:/DTU/Eksamensprojekt/bilag/
equal=.false.
2 if (k.lt.i) then
k=k+1
if ((abs((p−par_file(1,k))/p).lt.1.d−1).and.
$ (abs((h−par_file(2,k))/h).lt.1.d−2)) then
equal=.true.
endif
goto 2
endif
if ((.not.equal).and.(p.lt.500)) then
if (i.eq.(file_size−1)) then
write(UNIT=2, REC=file_size+1) file_size
rec_nr=i+1
elseif (i.eq.file_size) then
read(UNIT=2, REC=file_size+1) rec_nr
if (rec_nr.eq.file_size) then
rec_nr=1
else
rec_nr=rec_nr+1
endif
write(UNIT=2, REC=file_size+1) rec_nr
else
rec_nr=i+1
endif
c if (fiters) print*,’rec_nr: ’,rec_nr
write(UNIT=2,REC=rec_nr)P,h,t,v,s,x,u
endif
endif
3 continue
c if (fiters) print*,in(1),f(1),in(2),f(2)
c if (fiters) print*,’result’,t,v
return
end
subroutine iterer_omega_tau(f1,f2,tau,omega,in1,in2,
$ error)
integer in1,in2,i,error,IMAX,phase,step
double precision f1,f1_old,f1_new,f2,f2_old,f2_new,tau,tau_old,
$ tau_new,omega,omega_old,omega_new,dtau,domega,df1,df2,df
$ ,f1_new_omega,f1_new_tau,f2_new_omega,f2_new_tau,df1domega
$ ,df2domega,df1dtau,df2dtau,factor,rho_star
$ ,omega_l,omega_g,EPSV,R,P_c,T_c,T_ref,rho_ref,M,T_0C,const
logical found,fiters
common/constants/EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C,IMAX
$ ,fiters
omega_old=omega
tau_old=tau
call calculate(f1_old,tau_old,omega_old,in1)
call calculate(f2_old,tau_old,omega_old,in2)
omega_new=omega_old*(1.d0+1d−11)
tau_new=tau_old*(1.d0+1d−11)
call calculate(f1_new_omega,tau_old,omega_new,in1)
call calculate(f1_new_tau,tau_new,omega_old,in1)
call calculate(f2_new_omega,tau_old,omega_new,in2)
call calculate(f2_new_tau,tau_new,omega_old,in2)
i=0
factor=1
phase=0
step=3
found=.false.
5 if ((i .lt. IMAX) .and.(.not.found)) then
call calculate(f1_new,tau_new,omega_new,in1)
call calculate(f2_new,tau_new,omega_new,in2)
df1=f1_new−f1
df2=f2_new−f2
df=df1*df1+df2*df2
if ((abs(df1/f1).le.EPSV).and.(abs(df2/f2).le.EPSV)) then
found=.true.
goto 10
endif
if ((tau_new.eq.tau_old).and.(omega_new.eq.omega_old)) then
c if (fiters) print*
c $ ,’tau_new=tau_old and omega_new=omega_old ’,abs(df1/f1)
c $ ,abs(df2/f2),factor,f1,f2,phase
if ((abs(df1/f1).le.1.d−8).and.(abs(df2/f2).le.
$ 1.d−8)) then
found=.true.
6/19
19−03−2007