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

methanol.for
d:/DTU/Eksamensprojekt/bilag/
if (omega_new*rho_ref.lt.rho_star) then
omega_old=omega_new*0.999
else
omega_old=omega_new*1.001
endif
call calculate(f_old,tau,omega_old,in)
error=0
endif
df=f_new−f
if ((abs(df/f).le.EPSV).or.(abs(df/f).le.1.d−6.and.
$ omega_old.eq.omega_new)) then
found=.true.
goto 10
end if
if ((f_new−f_old).eq.0) then
dfdomega=(1D−12)/(omega_new−omega_old)
else
dfdomega=(f_new−f_old)/(omega_new−omega_old)
endif
c if (fiters) print*,in
c $ ,f,f_new,f_old,abs(df/f),1/(omega_old*rho_ref*M)*1.d−6,1
c $ /(omega_new*rho_ref*M)*1.d−6,T_ref/tau−T_0C,i
omega_old=omega_new
domega=−(df/dfdomega)
omega_new=omega_old+domega
if (omega_old.lt.1) then
if (domega.lt.(−(omega_old*.9))) then
omega_new=omega_old−omega_old*.9
elseif (domega.gt.(omega_old*9)) then
omega_new=omega_old+omega_old*9
endif
else
if ((abs(domega)).gt.(abs(omega_old*.05)))
$ omega_new=omega_old+domega/abs(domega)*omega_old*.05
endif
f_old=f_new
i=i+1
goto 5
else
goto 10
endif
10 if (.not.found) then
error=5
if (fiters) print*,’Warning: No convergence in iterer_omega’,in
$ ,f,f_new,abs(df/f),1/(omega_old*rho_ref*M)*1.d−6,1
$ /(omega_new*rho_ref*M)*1.d−6,T_ref/tau−T_0C,i
end if
omega=omega_new
c if (fiters) print*,’iterer_omega. Antal iterationer=’,i
return
end
subroutine iterer_tau(f,tau,omega,in,error)
integer in,i,error,IMAX
double precision f,tau,omega,tau_old,tau_new,dtau,f_new
$ ,f_old,df,dfdtau,EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C
logical found,fiters
common/constants/EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C,IMAX
$ ,fiters
tau_old=tau
call calculate(f_old,tau_old,omega,in)
tau_new=tau_old*1.01
i=0
found=.false.
5 if ((i.lt.IMAX).and.(.not.found)) then
call calculate(f_new,tau_new,omega,in)
df=f_new−f
if (abs(df/f).le.EPSV) then
found=.true.
goto 10
end if
dfdtau=(f_new−f_old)/(tau_new−tau_old)
tau_old=tau_new
dtau=−(df/dfdtau)
tau_new=tau_old+dtau
if ((abs(dtau)).gt.(abs(tau_old*.2))) tau_new=tau_old+dtau
$ /abs(dtau)*tau_old*.2
9/19
19−03−2007