Design og modellering af metanolanlæg til VEnzin-visionen Bilag
Design og modellering af metanolanlæg til VEnzin-visionen Bilag Design og modellering af metanolanlæg til VEnzin-visionen Bilag
methanol.for d:/DTU/Eksamensprojekt/bilag/ goto 1 endif else i=i−1 endif else OPEN(UNIT=2, FILE=filename, STATUS=’new’,ACCESS=’DIRECT’, RECL $ =8*7) i=0 endif omega=rho_start/rho_ref tau=T_ref/T_start c if (fiters) print*,’guess’,T_start−T_0C,1/(rho_start*M)*1.d−6 c $ ,value_old if ((in1.eq.1).and.(in2.eq.2).or.(in1.eq.2).and.(in2.eq.1)) $then c if (x_start.eq.1.d−7) then c print*,’x=0’ c omega=rho_star/rho_ref c call iterer_tau(P,tau,omega,1,error) c T=T_ref/tau c x=1.d−7 c call calculate_v(rho,T,x) c omega=rho/rho_ref c endif call iterer_omega_tau(P,h,tau,omega,1,2,error) T=T_ref/tau rho=omega*rho_ref call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.1).and.(in2.eq.3).or.(in1.eq.3).and.(in2.eq.1)) $then tau=T_ref/T call iterer_omega(P,tau,omega,1,error) call calculate_h(h,tau,omega) rho=omega*rho_ref call calculate_s(s,tau,omega) x=100 call calculate_u(u,tau,omega) elseif ((in1.eq.1).and.(in2.eq.4).or.(in1.eq.4).and.(in2.eq.1)) $then omega=rho/rho_ref call iterer_tau(P,tau,omega,1,error) call calculate_h(h,tau,omega) T=T_ref/tau call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.1).and.(in2.eq.5).or.(in1.eq.5).and.(in2.eq.1)) $then call iterer_omega_tau(P,s,tau,omega,1,5,error) call calculate_h(h,tau,omega) T=T_ref/tau rho=omega*rho_ref iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.1).and.(in2.eq.6).or.(in1.eq.6).and.(in2.eq.1)) $then omega=rho_star/rho_ref call iterer_tau(P,tau,omega,1,error) T=T_ref/tau call calculate_v(rho,T,x) omega=rho/rho_ref call calculate_h(h,tau,omega) call calculate_s(s,tau,omega) call calculate_u(u,tau,omega) elseif ((in1.eq.1).and.(in2.eq.7).or.(in1.eq.7).and.(in2.eq.1)) $then call iterer_omega_tau(P,u,tau,omega,1,7,error) call calculate_h(h,tau,omega) T=T_ref/tau rho=omega*rho_ref 3/19 19−03−2007
methanol.for d:/DTU/Eksamensprojekt/bilag/ call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) elseif ((in1.eq.2).and.(in2.eq.3).or.(in1.eq.3).and.(in2.eq.2)) $then tau=T_ref/T call iterer_omega(h,tau,omega,2,error) call calculate_P(P,tau,omega) rho=omega*rho_ref call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.2).and.(in2.eq.4).or.(in1.eq.4).and.(in2.eq.2)) $then omega=rho/rho_ref call iterer_tau(h,tau,omega,2,error) call calculate_P(P,tau,omega) T=T_ref/tau call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.2).and.(in2.eq.5).or.(in1.eq.5).and.(in2.eq.2)) $then call iterer_omega_tau(h,s,tau,omega,2,5,error) call calculate_P(P,tau,omega) T=T_ref/tau rho=omega*rho_ref iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.2).and.(in2.eq.6).or.(in1.eq.6).and.(in2.eq.2)) $then call iterer_omega_tau(h,x,tau,omega,2,6,error) call calculate_P(P,tau,omega) T=T_ref/tau rho=omega*rho_ref call calculate_s(s,tau,omega) call calculate_u(u,tau,omega) elseif ((in1.eq.2).and.(in2.eq.7).or.(in1.eq.7).and.(in2.eq.2)) $then call iterer_omega_tau(h,u,tau,omega,2,7,error) call calculate_P(P,tau,omega) T=T_ref/tau rho=omega*rho_ref call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) elseif ((in1.eq.3).and.(in2.eq.4).or.(in1.eq.4).and.(in2.eq.3)) $then omega=rho/rho_ref tau=T_ref/T call calculate_P(P,tau,omega) call calculate_h(h,tau,omega) call calculate_s(s,tau,omega) iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.3).and.(in2.eq.5).or.(in1.eq.5).and.(in2.eq.3)) $then tau=T_ref/T call iterer_omega(s,tau,omega,5,error) call calculate_P(P,tau,omega) call calculate_h(h,tau,omega) rho=omega*rho_ref iterate_rho_l=.true. call calculate_x(x,tau,omega,omega_g,omega_l) call calculate_u(u,tau,omega) elseif ((in1.eq.3).and.(in2.eq.6).or.(in1.eq.6).and.(in2.eq.3)) $then tau=T_ref/T call calculate_v(rho,T,x) omega=rho/rho_ref call calculate_P(P,tau,omega) call calculate_h(h,tau,omega) call calculate_s(s,tau,omega) call calculate_u(u,tau,omega) 4/19 19−03−2007
- Page 291 and 292: VEnzin.for c:/dna/source/ C Subrout
- Page 293 and 294: VEnzin.for c:/dna/source/ CALL STAT
- Page 295 and 296: VEnzin.for c:/dna/source/ CA 3: Flu
- Page 297 and 298: VEnzin.for c:/dna/source/ C C C M_B
- Page 299 and 300: VEnzin.for c:/dna/source/ C C SETFL
- Page 301 and 302: VEnzin.for c:/dna/source/ C C GASCO
- Page 303 and 304: VEnzin.for c:/dna/source/ RES(5) =
- Page 305 and 306: VEnzin.for c:/dna/source/ c c c E2=
- Page 307 and 308: VEnzin.for c:/dna/source/ C C HEATE
- Page 309 and 310: VEnzin.for c:/dna/source/ $ ’$\\d
- Page 311 and 312: VEnzin.for c:/dna/source/ c 1 = Wat
- Page 313 and 314: VEnzin.for c:/dna/source/ CA 4: Fin
- Page 315 and 316: VEnzin.for c:/dna/source/ CL K_MED
- Page 317 and 318: VEnzin.for c:/dna/source/ CA MEDIE
- Page 319 and 320: VEnzin.for c:/dna/source/ CA −4 :
- Page 321 and 322: VEnzin.for c:/dna/source/ C ENDDO E
- Page 323 and 324: VEnzin.for c:/dna/source/ MEDIE(2)
- Page 325 and 326: VEnzin.for c:/dna/source/ C C 400 C
- Page 327 and 328: VEnzin.for c:/dna/source/ GOTO 9999
- Page 329 and 330: VEnzin.for c:/dna/source/ MEDIE(1)
- Page 331 and 332: VEnzin.for c:/dna/source/ C−−
- Page 333 and 334: VEnzin.for c:/dna/source/ RES(1) =
- Page 335 and 336: VEnzin.for c:/dna/source/ C−−
- Page 337 and 338: VEnzin.for c:/dna/source/ C−−
- Page 339 and 340: 35. Metanol (fluid) til DNA - Fortr
- Page 341: methanol.for d:/DTU/Eksamensprojekt
- Page 345 and 346: methanol.for d:/DTU/Eksamensprojekt
- Page 347 and 348: methanol.for d:/DTU/Eksamensprojekt
- Page 349 and 350: methanol.for d:/DTU/Eksamensprojekt
- Page 351 and 352: methanol.for d:/DTU/Eksamensprojekt
- Page 353 and 354: methanol.for d:/DTU/Eksamensprojekt
- Page 355 and 356: methanol.for d:/DTU/Eksamensprojekt
- Page 357 and 358: methanol.for d:/DTU/Eksamensprojekt
methanol.for<br />
d:/DTU/Eksamensprojekt/bilag/<br />
call calculate_s(s,tau,omega)<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
elseif ((in1.eq.2).and.(in2.eq.3).or.(in1.eq.3).and.(in2.eq.2))<br />
$then<br />
tau=T_ref/T<br />
call iterer_omega(h,tau,omega,2,error)<br />
call calculate_P(P,tau,omega)<br />
rho=omega*rho_ref<br />
call calculate_s(s,tau,omega)<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.2).and.(in2.eq.4).or.(in1.eq.4).and.(in2.eq.2))<br />
$then<br />
omega=rho/rho_ref<br />
call iterer_tau(h,tau,omega,2,error)<br />
call calculate_P(P,tau,omega)<br />
T=T_ref/tau<br />
call calculate_s(s,tau,omega)<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.2).and.(in2.eq.5).or.(in1.eq.5).and.(in2.eq.2))<br />
$then<br />
call iterer_omega_tau(h,s,tau,omega,2,5,error)<br />
call calculate_P(P,tau,omega)<br />
T=T_ref/tau<br />
rho=omega*rho_ref<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.2).and.(in2.eq.6).or.(in1.eq.6).and.(in2.eq.2))<br />
$then<br />
call iterer_omega_tau(h,x,tau,omega,2,6,error)<br />
call calculate_P(P,tau,omega)<br />
T=T_ref/tau<br />
rho=omega*rho_ref<br />
call calculate_s(s,tau,omega)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.2).and.(in2.eq.7).or.(in1.eq.7).and.(in2.eq.2))<br />
$then<br />
call iterer_omega_tau(h,u,tau,omega,2,7,error)<br />
call calculate_P(P,tau,omega)<br />
T=T_ref/tau<br />
rho=omega*rho_ref<br />
call calculate_s(s,tau,omega)<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
elseif ((in1.eq.3).and.(in2.eq.4).or.(in1.eq.4).and.(in2.eq.3))<br />
$then<br />
omega=rho/rho_ref<br />
tau=T_ref/T<br />
call calculate_P(P,tau,omega)<br />
call calculate_h(h,tau,omega)<br />
call calculate_s(s,tau,omega)<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.3).and.(in2.eq.5).or.(in1.eq.5).and.(in2.eq.3))<br />
$then<br />
tau=T_ref/T<br />
call iterer_omega(s,tau,omega,5,error)<br />
call calculate_P(P,tau,omega)<br />
call calculate_h(h,tau,omega)<br />
rho=omega*rho_ref<br />
iterate_rho_l=.true.<br />
call calculate_x(x,tau,omega,omega_g,omega_l)<br />
call calculate_u(u,tau,omega)<br />
elseif ((in1.eq.3).and.(in2.eq.6).or.(in1.eq.6).and.(in2.eq.3))<br />
$then<br />
tau=T_ref/T<br />
call calculate_v(rho,T,x)<br />
omega=rho/rho_ref<br />
call calculate_P(P,tau,omega)<br />
call calculate_h(h,tau,omega)<br />
call calculate_s(s,tau,omega)<br />
call calculate_u(u,tau,omega)<br />
4/19<br />
19−03−2007