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

VEnzin.for
c:/dna/source/
K_MED Media description.
C
C Subroutines : COMINF
C MIXER
C
CP Programmer : Bent Lorentzen 1994 (Niels Emsholm 1991)
CP Lab. for Energetics, DTH, Denmark.
C***********************************************************************
C
C Including the common "environment"
C
INCLUDE ’ENVIRO.INI’
INCLUDE ’THERPROP.DEC’
C
C Parameter variables
C
INTEGER ANTLK, ANTKN, MEDIE(10),
: ANTM1, ANTME, VARME(1), ANTEL(1),
: VAREL(ANTST,1),MMVAR(MAXMM)
DOUBLE PRECISION X_J(MAXME,ANTST), RES(60), MDOT(10), P(10)
CHARACTER*80 KOMTY
C
C Local variables
C
INTEGER I,J,MAXKN
DOUBLE PRECISION XMIX(ANTST),T_SUM,M_BL(10),N(ANTST,10)
$ ,N_total(ANTST)
CHARACTER*500 K_LIG(40), K_BET
CHARACTER*1000 KOMDSC,K_INP
CHARACTER*100 KMEDDS(3)
EXTERNAL COMINF,MIXER
INCLUDE ’THERPROP.INI’
C=======================================================================
GOTO (100,200,1,400,400,200) FKOMP
1 RETURN
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component name
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
100 CONTINUE
KOMTY = ’MIXER_03’
MMVAR(1) = 10
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component characteristics
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
200 CONTINUE
KOMTY = ’MIXER_03’
ANTKN = MMVAR(1)
ANTLK = MMVAR(1)+37
ANTM1 = MMVAR(1)
DO I=1,MMVAR(1)
MEDIE(I) = ANYGAS$
ENDDO
ANTME = 1
VARME(1) = −MMVAR(1)
ANTEL(1) = 38
DO I=1,36
VAREL(I,1) = I
ENDDO
VAREL(37,1) = 38
VAREL(38,1) = 39
IF (FKOMP.EQ.6) GOTO 600
*** FKOMP = 3
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component equations. All in residual form.
C Do not include the conservation laws, since these are treated
C automatically by DNA.
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
400 CONTINUE
C
DO J=1, ANTKN−1
M_BL(J) = 0.D0
ENDDO
DO J=1, ANTKN−1
DO I=1, ANTST
M_BL(J) = M_BL(J) + X_J(MEDIE(J),I)*M_MOL(I)
N_total(I)=0
49/67
19−03−2007