bearjrgm
Advanced Member | Редактировать | Профиль | Сообщение | ICQ | Цитировать | Сообщить модератору
C***********************************************************************
C ENVIRO.INI contains the common environment for all subroutines. This
C environment is parameters and common blocks. Finally an implicit none
C statement ensures that all variables have to be declared.
C***********************************************************************
C
CC IFILD Unit number for input file (system description).
CC IFILP Unit number for file SIMPAR.DAT (numerical parameters).
CC UDFIL Unit number for output file.
CC UDDAT Unit number for listing file.
CC UDMED Unit number for file holding fluid compositions.
CC FILUD File name for output file.
CC FILLS File name for listing file.
CC FILME File name for file holding fluid compositions.
CC W Work space for double precision variables.
CC IW Work space for integer variables.
CC AW Work space for character variables.
CC PWNEX Pointer to first empty space in W.
CC PIWNE Pointer to first empty space in IW.
CC PAWNE Pointer to first empty space in AW.
CC PSYSK Pointer to SYSKN in IW.
CC PSYSP Pointer to SYSPA in W.
CC PSPTA Pointer to SPTAB in IW.
CC PKMTA Pointer to KMTAB in AW.
CC PDYTA Pointer to DYTAB in AW.
CC PLITA Pointer to LITAB in IW.
CC PRETA Pointer to RETAB in IW.
CC PSYSR Pointer to SYSRE in IW.
CC PRE Pointer to RE in W.
CC PY Pointer to Y in W.
CC PRATA Pointer to RATAB in IW.
CC PRAPA Pointer to RAPAR in W.
CC PINTA Pointer to INTAB in IW.
CC PINPA Pointer to INPAR in W.
CC PIVTA Pointer to IVTAB in IW.
CC PIVPA Pointer to IVPAR in W.
CC PRADA Pointer to RADAT in IW.
CC POPPA Pointer to OPPAR in W.
CC PSHTA Pointer to SHTAB in IW.
CC PSHPA Pointer to SHPAR in W.
CC PTRTA Pointer to TRTAB in IW.
CC PX_J Pointer to X_J in W.
CC PBETI Pointer to BETIN in IW.
CC PMETA Pointer to METAB in IW.
CC PRW Pointer for double precision numerical work space.
CC PIWRK Pointer for integer numerical work space.
CC DSYSK Dimension of SYSKN in IW.
CC DSYSP Dimension of SYSPA in W.
CC DSPTA Dimension of SPTAB in IW.
CC DKMTA Dimension of KMTAB in AW.
CC DDYTA Dimension of DYTAB in AW.
CC DSHTA Dimension of SHTAB in IW.
CC DSHPA Dimension of SHPAR in W.
CC DLITA Dimension of LITAB in IW.
CC DRETA Dimension of RETAB in IW.
CC DSYSR Dimension of SYSRE in IW.
CC DRE Dimension of RE in W.
CC DY Dimension of Y in W.
CC DRATA Dimension of RATAB in IW.
CC DRAPA Dimension of RAPAR in W.
CC DINTA Dimension of INTAB in IW.
CC DINPA Dimension of INPAR in W.
CC DIVTA Dimension of IVTAB in IW.
CC DIVPA Dimension of IVPAR in W.
CC DRADA Dimension of RADAT in IW.
CC DOPPA Dimension of OPPAR in W.
CC DTRTA Dimension of TRTAB in IW.
CC DMETA Dimension of METAB in IW.
CC DX_J Dimension of X_J in W.
CC DRW Dimension of double precision numerical work space.
CC DIWRK Dimension of integer work space.
CC DBETI Dimension of BETIN in IW.
CC DIMSY One dimension of SYSRE.
CC DIMME One dimension of METAB.
CC ANTKM Number of components in actual system.
CC ANTTY Number of components defined.
CC ANTTR Number of nodes in actual system.
CC ANTBE Maximum number of branches for a single component in the actual
CC system.
CC ANTSH Number of dynamic shafts.
CC ANTPA Number of parameters for all components in the actual system.
CC ANTRE Number of residuals incl. "one-branch-node" residuals but excl.
CC additional conditions.
CC ANTVA Number of variables incl. repeated pressures.
CC ANTLI Number of equations.
CC ANTDI Number of differential equations.
CC ANTDY Number of dynamic variables.
CC ANTFR Number of residuals excl. "one-branch-node" residuals and additional
CC conditions.
CC ANTRA Number of additional conditions.
CC ANTIN Number of intial guesses in input file.
CC ANTIV Number of intial values for dynamic variables in input file.
CC ANTOP Number of additional conditions to be updated in each iteration
CC step.
CC ANTST Number of compounds.
CC ANTIL Number of closed internal loops (algebraic loops).
CC ANTSH Number of dynamic shafts.
CC PROKM Number of predefined components.
CC MAXPA Maximum number of component parameters.
CC MAXKM Maximum number of components in a system.
CC MAXTY Maximum number of different components defined.
CC MAXBE Maximum number of branches for a component.
CC MAXRA Maximum number of additional conditions.
CC MAXME Maximum number of fluids.
CC MAXLI Maximum number of equations.
CC MAXRE Maximum number of residuals incl. "one-branch-node" residuals but excl.
CC additional conditons.
CC MAXVA Maximum number of variables incl. repeated pressures.
CC MAXDY Maximum number of dynamic variables.
CC MAXW Maximum number of places in work space W.
CC MAXIW Maximum number of places in work space IW.
CC MAXAW Maximum number of places in work space AW.
CC FBETI Flag for checking conditions in each component :
CC TRUE Yes , FALSE No.
CC FINIT Flag for initializing user defined components :
CC TRUE : First user defined component to be initialized.
CC FALSE : Not the first.
CC FKOMP Flag with the value:
CC 1: Initialize the component.
CC 2: Initialize with actual system.
CC 3: Fluid composition calculation (constant).
CC 4: Find residuals.
CC 5: Find residuals and check variables.
CC 6: Output information about component.
CC ISMAX Maximum number of iterations in THERMO.
CC EPSST Stop criteria for iteration in THERMO.
CC IVMAX Maximum number of iterations in STETAB.
CC EPSVA Stop criteria for iteration in STETAB.
CC RESMAX Maksimum value for residuals in NOLEQ (divergence).
CC RESOK Residuals OK in NOLEQ.
CC DISOK Displacements OK in NOLEQ.
CC TOLINI Drop tolerance for elements in matrices.
CC RECINI Damping factor for displacement during iteration in NOLEQ.
CC RECCHA Factor for change of damping factor in each iteration step in
CC NOLEQ.
CC ITOPD Number of iterations between two updates of the Jacobian in
CC NOLEQ.
CC MAXIT Maksimum number of Newton iterations in NOLEQ.
CC PERTUW Weight used to pertubate variables for finding the Jacobian.
CC PERTU0 Pertubation of variables with the value zero.
CC PASC Pascal matrix (calculate predicted value of Nordsieck-vector).
CC LVEK L-vector (calculate corrected value of Nordsieck-vector).
CC FK Factor to find truncation error for pred.-corr. pair.
CC FKOS Change of order of BDF-method.
CC MSTEP Maximum number of timesteps.
CC VORD Chosen order of BDF-method.
CC H0 Initial step length.
CC HMAX Maximum step length.
CC MTIME Maximum actual time.
CC KONST Constant used for conservative choise of new step length.
CC EPSBDF Maximum allowed truncation error.
CC TIME Actual time.
CC STEP Step length.
CC ZNORDP Predicted value of the Nordsieck-vector.
CC ORDEN Order of the BDF method.
C DNA21 CC DIVMAX Maximum number of rejected step lengths.
CC DIVNEW Maximum error increase between two Newton steps.
CC FDIV Flag indicating divergence for chosen timestep
CC TRUE if divergence.
CC FPAU Flag indicating whether warnings should hold the screen
CC TRUE if hold screen.
C
CBE
CC PYPTA pointer for YPTAB in IW.
CC DYPTA dimension of YPTAB in IW.
CC ANTVP Number of variable parameters.
CC UDSTR Unit number for input storing/debugging file.
CC UDID Unit number for ID-file.
CC FILST Filename for input storing/debugging file.
CC FILID ID-file to use when converting dynamic results from machine
CC code to ASCII format.
CC ANTSOL Counter incremented if SOLVE command is given in input file.
CC LINIE Line number where SOLVE command is given in input file.
CC FSOLV Flag set when SOLVE command is given.
CC FFACT Flag set in case of error during factorization.
CC FOUT Flag set when dynamic output mode
C wanted, clear for static output mode
CC FTEM
CC PWNE0 Pointer for end of real workspace used in TABINP
CC PAWN0 Pointer for end of character workspace used in TABINP
CC PIWN0 Pointer for end of integer workspace used in TABINP
CC ANTAN Number of changes made in tables after a SOLVE
CC PANTA Pointer for ANTAB in IW
CC DANTA Dimension for ANTAB in IW
CC MAXAN Maximum number of changes in system
CC OUTS
CC OUTD
CC OUTM
CC SCRATS
CC SCRATM
CC SCRATD
CC UDEBUG
CC FWARN False turns warnings off in F. Used when perturbing.
C
C
C dna 21
CC MINUPD Number of Jacobain updates for divergence test.
CC FDISK Flag indicating discontinuity
CC TRUE if discontinuity.
CC TDISK Discontinuity time for nearest discontinuity.
C
CP Programmer : Bent Lorentzen 1995 (Claus Perstrup 1990)
CP Lab. for Energetics, DTU, Denmark.
CBE
CP Edited by Brian Elmegaard LfE 1995
C***********************************************************************
C
C All variables must be declared
C
IMPLICIT NONE
C
C Parameters setting maximum sizes of global arrays
C
C Maximum number of integers in work space (Length of IW)
INTEGER MAXIW
PARAMETER(MAXIW=4000000)
C Integer work space
INTEGER IW(MAXIW)
C Maximum number of double precisions in work space (Length of W)
INTEGER MAXW
PARAMETER(MAXW =4000000)
C Double precision work space
DOUBLE PRECISION W(MAXW)
C Maximum number of characters in work space (Length of AW)
INTEGER MAXAW
PARAMETER(MAXAW=200000)
c Character work space
CHARACTER*80 AW(MAXAW)
C Maximum number of component types (Length of KMTAB, SPTAB)
INTEGER MAXTY
PARAMETER(MAXTY=600)
C Maximum allowed number of components (Length of SYSKN)
INTEGER MAXKM
PARAMETER(MAXKM=200)
C Maximum allowed number of branches on a component (Used in SYSKN)
INTEGER MAXBE
PARAMETER(MAXBE=30)
C Maximum allowed number of parameters (length of SYSPA)
INTEGER MAXPA
PARAMETER(MAXPA=MAXKM*MAXBE*2)
C Maximum allowed number of additional conditions (length of RATAB,
C RAPAR)
INTEGER MAXRA
PARAMETER(MAXRA=100)
C Maximum number of media composed in the model, i.e., specficied by
C SOLID or FLUID keywords or in a component and by X_J guesses
C (Width of X_J)
INTEGER MAXME
PARAMETER(MAXME=80)
C Number of chemical elements defined in THERPROP (Length of X_J)
INTEGER ANTST
PARAMETER(ANTST=39)
C Maximum allowed number of variables in a system (Length of LITAB)
INTEGER MAXVA
PARAMETER(MAXVA=MAXKM*MAXBE*2)
C Maximum allowed number of equations in a system (Length of RETAB)
INTEGER MAXRE
PARAMETER(MAXRE=MAXKM*MAXBE*2)
C Maximum allowed number of equations needed to simulate a system
C (Length of Y)
INTEGER MAXLI
PARAMETER(MAXLI=MAXKM*MAXBE*2)
C Maximum allowed number of dynamic variables in a system
C (Length of ZNORDP, ZNORDC)
INTEGER MAXDY
PARAMETER(MAXDY=MAXKM*4)
C Maximum allowed number of variable parameters in a system (Length
C of YPTAB)
INTEGER MAXVP
PARAMETER(MAXVP=MAXPA)
C Maximum allowed number of changes made after to a model during
C parameter variations (Length of ANTAB)
INTEGER MAXAN
PARAMETER(MAXAN=50000)
C Maximum allowed number of SIMUL's allowed
INTEGER MAXSOL
PARAMETER(MAXSOL=10000)
C Maximum allowed number of state calls allowed (Length of SCTAB)
INTEGER MAXSC
PARAMETER(MAXSC=100)
C Maximum allowed number of state shifts allowed (Length of ISTAB)
INTEGER MAXIS
PARAMETER(MAXIS=100)
C Maximum allowed number of names (Length of NATAB, NAMES)
INTEGER MAXNA
PARAMETER(MAXNA=MAXKM+MAXME+MAXVA)
C Maximum allowed number of configuration parameters of a component
C (width of MMTAB)
INTEGER MAXMM
PARAMETER(MAXMM=1)
C Maximum number of refrigerants in use in a model
INTEGER MAXREF
PARAMETER(MAXREF=10)
C Maximum allowed loops in a model (Unused, intended for use in
C LOPTAB)
INTEGER MAXLOP
PARAMETER(MAXLOP=100)
CBE
C
C Maximum allowed values in each component
C
INTEGER COMPAR, COMZA, COMZD, COMZS, COMRES, COMZC
PARAMETER (COMPAR = 100)
PARAMETER (COMZA = 300)
PARAMETER (COMZD = 20)
PARAMETER (COMZS = 5)
PARAMETER (COMRES = 200)
PARAMETER (COMZC = 5)
C
C Size of work spaces
C
C
C Global variables
C
C BE
C
INTEGER IFILD
INTEGER IFILP, UDFILS, UDDAT, UDMED, DRAPA, ANTRA, ANTFR,
$ PSYSK, PSYSP, PIWNE, PSPTA, PWNEX, PKMTA, PAWNE, PLITA, PRETA
$ , PSYSR, PRE , PY , PRATA, PRAPA, PMETA, PIVTA, PINTA,
$ PINPA, PRADA, POPPA, PTRTA, PX_J , PBETI, PIVPA, PDYTA, PSHTA
$ , PSHPA, ANTIL, DIMSY, DIMME, DSHTA, DSHPA, DSYSP, DSPTA,
$ DSYSK, DKMTA, DLITA, DRETA, DSYSR, DRE, DINTA, DINPA, DY ,
$ DRADA, DOPPA, DTRTA, DRATA, DMETA, DX_J , DIVTA, DIVPA, DBETI
$ , DDYTA, PROKM, FKOMP, ANTTR, ANTKM, ANTTY, ANTRE, ANTVA,
$ ANTLI, ANTIN, ANTOP,
C DNA21
& ANTBE, ANTPA, PRW, DRW, PIWRK, DIWRK, MINUPD
C DNA21
INTEGER ISMAX, IVMAX, ITOPD, MAXIT, VORD, MSTEP, ANTIV,
& ORDEN, ANTDY, ANTZD, ANTSH,
CBE
& PYPTA,DYPTA,ANTVP,UDSTR,LINIE,ANTSOL,UDFILD,PIWN0,
& PWNE0,PAWN0,UDID,ANTAN,PANTA,DANTA,PANPA,DANPA,PTITL,
& DTITL,OUTS,OUTD,OUTM,SCRATS,SCRATM,SCRATD,UDEBUG,
& PYOLD,DYOLD,PPAR0,DPAR0,PWORK,DWORK,MAXHAL,PINEX,DINEX,
& ANTMD,PMDTA,DMDTA,ANTIE,ANTVM,PVMTA,DVMTA,PCP,DCP,
& PHV,DHV,PHF,DHF,UFILD,UFILP,SCRATU,ANTSC,PSCTA,DSCTA,
& PSCPA,DSCPA,PMFTA,DMFTA,PMMTA,DMMTA,ANTMM,PISTA,DISTA,
& ANTIS,PISPA,DISPA,PSWTA,DSWTA,ANTSD,PNATA,DNATA,ANTNA,
& PNAME,DNAME,ANTMF,MINITO,UFUPD,PLOPTA,DLOPTA,PLOPPA,
& DLOPPA,ANTLOP,PREFTA,DREFTA,ANTREF,ANTGP,PGENPA,DGENPA,
& ANTPP,PPAPOI,DPAPOI,ITER,ANTOLD,PWCHA,DWCHA,erpar1,
& erpar2
C
DOUBLE PRECISION EPSST,EPSVA,RESMAX,RESOK,TOLINI,DISOK,RECINI,
& RECCHA,PERTUW,PERTU0,PASC(5,5),LVEK(4,5),FK(4),
C DNA21
& FKOS(3),TIME,ZNORDP(MAXDY,5),STEP,TDISK,tdiskmex,
CBE
& H0,HMAX,MTIME,EPSBDF,KONST,HMIN,DIVNEW,HINC,
& RESINC,DISINC,BESTAC,BESTDS,BESTRS,STPRAT,
& TIMEOLD,PDEAD,TDEAD,RDISOK
C
C DNA21
LOGICAL FBETI, FINIT, FDIV, FPAU, FDISK,
CBE
& FSOLV,FFACT, FOUTD, FOUTS,FTEM, FLSET, FLSTO,
& FSTRU,FDFDE,FALLOW,FWARN,FEXS,FEXD,FEXM ,FSTS,
& FDEBUG,FXERGY,FITER,FCHA,FCOMLI,FMEX,FITER0,
& FIGNBT
C
CBE
CHARACTER*80 FILUDS, FILUDD, FILME, FILLS, FILST, FILID,
& FILPAR, FILUPD, FILIN, FILINC,DATVER
CBE 050816
CHARACTER*2 CLNEND
C
C Unit numbers
C
CBE
COMMON/NUNIT/IFILD,IFILP,UDFILS,UDFILD,UDDAT,UDMED,UDSTR,UDID,
& OUTS,OUTD,OUTM,SCRATS,SCRATM,SCRATD,UFILD,
& UFILP,SCRATU,UFUPD
C
C Unit names
C
CBE
COMMON/NAMUN/FILUDS,FILUDD,FILME,FILLS,FILST,FILID,FILPAR,FILUPD,
& FILIN,FILINC
CBE Data file version
C
PARAMETER(DATVER='060124')
C
CBE
C Solve and line number in input file
COMMON/INFIL/ANTSOL,LINIE
C
C Double precision work space
C
COMMON/WORKSP/ W
C
C Pointers to W
C
CBE
COMMON/PEGWSP/ PSYSP, PWNEX, PRE, PY, PRAPA, PRW, PINPA,
& POPPA, PSHPA, PX_J, PIVPA, PWNE0, PANPA, PYOLD,
& PPAR0, PWORK, PCP, PHV, PHF, PSCPA, PISPA, PLOPPA,
& PWCHA
C
C Integer work space
C
COMMON/IWORKS/ IW
C
C Pointers to IW
C
COMMON/PEGIWS/ PIWNE, PSYSK, PSPTA, PLITA, PRETA, PSYSR, PIVTA,
& PRATA, PIWRK, PINTA, PRADA, PTRTA, PBETI, PMETA,
CBE
& PSHTA, PYPTA, PIWN0, PANTA, PINEX, PMDTA, PVMTA,
& PSCTA, PMFTA, PMMTA, PISTA, PSWTA, PNATA, PLOPTA,
& PREFTA, PGENPA,PPAPOI
C Character work space
C
COMMON/AWORKS/ AW
C
C Pointers to AW
C
CBE
COMMON/PEGAWS/ PKMTA, PDYTA, PAWNE, PAWN0, PTITL, PNAME
C
C Dimensions of variables in W
C
COMMON/DIMWSP/ DSYSP, DRE, DY, DRAPA, DRW, DINPA, DOPPA,
& DIVPA, DSHPA, DX_J, DANPA, DYOLD, DPAR0,
& DWORK, DCP, DHV, DHF, DSCPA, DISPA, DLOPPA,
& DWCHA
C
C Dimensions of variables in IW
C
COMMON/DIMIWS/ DSYSK, DSPTA, DLITA, DRETA, DSYSR, DRATA, DIVTA,
CBE
& DIWRK, DINTA, DRADA, DTRTA, DBETI, DMETA, DSHTA,
& DYPTA, DANTA, DINEX, DMDTA, DVMTA, DSCTA, DMFTA,
& DMMTA, DISTA, DSWTA, DNATA, DLOPTA, DREFTA, DGENPA,
& DPAPOI
C
C Dimensions of variables in AW
C
CBE
COMMON/DIMAWS/ DKMTA, DDYTA, DTITL, DNAME
C
C Describing the system
C
COMMON/FORMAT/ ANTKM, ANTTY, ANTTR, ANTRE, ANTVA, ANTLI, ANTRA,
& ANTFR, PROKM, ANTIN, ANTOP, ANTPA, ANTBE, ANTIL,
CBE
C DNA21
& DIMSY, DIMME, ANTIV, ANTDY, ANTSH, ANTVP, ANTAN,
& ANTMD, ANTIE, ANTVM, ANTSC, ANTMM, ANTZD, ANTIS,
& ANTSD, ANTNA, ANTMF, ANTLOP, ANTREF, ANTGP, ANTPP,
& ANTOLD
C
C Logicals
C
CBE
C DNA21
COMMON/CHECK/ FBETI, FINIT, FPAU, FSOLV, FFACT,
& FOUTD,FOUTS, FTEM, FLSET, FLSTO,FSTRU,FDFDE,FIGNBT
COMMON/ERPAR/ ERPAR1,ERPAR2
COMMON/EXIST/FEXS,FEXD,FEXM
COMMON/DISK/ FDISK,FALLOW, FCHA
COMMON/WARN/ FWARN
COMMON/ITER/ FITER,FITER0
COMMON/DIV/ FDIV,FSTS
COMMON/DEBUG/ UDEBUG,FDEBUG
COMMON/EXERGY/ TDEAD,PDEAD, FXERGY
COMMON/COMLI/ FCOMLI
COMMON/FMEX/ FMEX
COMMON/CMEX/ clnend
C
C Controlling component calculation
C
COMMON/KOMPL/ FKOMP
C
C Maximum number of iterations
C
C DNA21
COMMON/MAXIT/ ISMAX, IVMAX, ITOPD, MAXIT, MINUPD, MAXHAL, MINITO,
: ITER
C
C Stop criterias for iterations
C
COMMON/STOP/ EPSST, EPSVA, RESMAX, RESOK, DISOK, DIVNEW,
: BESTAC, BESTDS, BESTRS, RDISOK
C
C Parameters to the modified Newton iteration
C
COMMON/TOLE/ TOLINI,RECINI,RECCHA,PERTUW,PERTU0,RESINC,DISINC
C
C Constants used for BDF-integration
C
COMMON/BDF1/ PASC,LVEK,FK,FKOS
C
C Variables used for BDF-integration
C
C DNA21
COMMON/BDF2/ TIME,STEP,ZNORDP,TDISK,STPRAT,TIMEOLD,tdiskmex,ORDEN
C
C Constants used for BDF-integration
C
C DNA21
COMMON/BDF3/ H0,HMAX,MTIME,EPSBDF,KONST,HMIN,HINC
C
C Constants used for BDF-integration
C
C DNA21
COMMON/BDF4/ VORD,MSTEP
CBE
C
C 031220 BE: Declaration of integer shortcuts for fluids
C
INTEGER H2$, O2$, N2$, CO$, NO$, CO2$ ,H2O_G$,NH3$,H2S$,
$ SO2$,CH4$,C2H6$,C3H8$,C4H10_N$,C4H10_I$,C5H12$,
$ C6H14$,C7H16$,C8H18$,C2H4$,C3H6$,C5H10$,C6H12_1$,
$ C7H14$,C2H2$,C6H6$,C6H12_C$,C$,S$,NO2$,HCN$,COS$,
$ N2O$,NO3$,SO3$,AR$,H2O_L$,ASH$,TAR$,H$,O$,N$,W$
PARAMETER(H2$ = 1)
PARAMETER(O2$ = 2)
PARAMETER(N2$ = 3)
PARAMETER(CO$ = 4)
PARAMETER(NO$ = 5)
PARAMETER(CO2$ = 6)
PARAMETER(H2O_G$ = 7)
PARAMETER(NH3$ = 8)
PARAMETER(H2S$ = 9)
PARAMETER(SO2$ = 10)
PARAMETER(CH4$ = 11)
PARAMETER(C2H6$ = 12)
PARAMETER(C3H8$ = 13)
PARAMETER(C4H10_N$ = 14)
PARAMETER(C4H10_I$ = 15)
PARAMETER(C5H12$ = 16)
PARAMETER(C6H14$ = 17)
PARAMETER(C7H16$ = 18)
PARAMETER(C8H18$ = 19)
PARAMETER(C2H4$ = 20)
PARAMETER(C3H6$ = 21)
PARAMETER(C5H10$ = 22)
PARAMETER(C6H12_1$ = 23)
PARAMETER(C7H14$ = 24)
PARAMETER(C2H2$ = 25)
PARAMETER(C6H6$ = 26)
PARAMETER(C6H12_C$ = 27)
PARAMETER(C$ = 28)
PARAMETER(S$ = 29)
PARAMETER(NO2$ = 30)
PARAMETER(HCN$ = 31)
PARAMETER(COS$ = 32)
PARAMETER(N2O$ = 33)
PARAMETER(NO3$ = 34)
PARAMETER(SO3$ = 35)
PARAMETER(AR$ = 36)
PARAMETER(H2O_L$ = 37)
PARAMETER(ASH$ = 38)
PARAMETER(TAR$ = 39)
PARAMETER(H$ = 1)
PARAMETER(O$ = 2)
PARAMETER(N$ = 3)
PARAMETER(W$ = 7)
INTEGER REALFL$,SOLUT$,CONSOL$,ANYSOL$,ANYGAS$,ANYFLU$,CONGAS$,
$ CONFLU$ ,SIMAIR$,COAL$,NATGAS$,METHANE$,NITROG$,OXYGEN$
$ ,GASWAT$ ,LIQWAT$,LIBRW$,WAT67$,WATHF$,REALCO2$,WATER$
$ ,MASSFLOW$, SHAFT$, DSHAFT$,DSHAFTDT$,POWER$,HEAT$,DHEAT$
$ ,MOLEFRAC$,MASSFRAC$ ,PRESSURE$,ENTHALPY$,CONC$,MASS$,MASS1$
$ ,MASS2$,MASSDT$,MASS1DT$,MASS2DT$ , MASSW$,MASS1W$ ,MASS2W$
$ ,INTENER$,INTENER1$,INTENER2$,INTENERDT$,INTENER1DT$
$ ,INTENER2DT$ ,INTENERW$,INTENER1W$,INTENER2W$,INTENERS$
$ ,INTENERSDT$,INTENERSW$ ,DYNVAR$,DYNVARDT$,DYNVARW$,ALGVAR$
$ ,CONTROL$,VARPAR$, STEAM$
PARAMETER(REALFL$=-9)
PARAMETER(SOLUT$=-8)
PARAMETER(CONSOL$=-6)
PARAMETER(ANYSOL$=-5)
PARAMETER(ANYGAS$=-4)
PARAMETER(ANYFLU$=-3)
PARAMETER(CONGAS$=-2)
PARAMETER(CONFLU$=-1)
PARAMETER(SIMAIR$=1)
PARAMETER(COAL$=2)
PARAMETER(NATGAS$=3)
PARAMETER(METHANE$=4)
PARAMETER(NITROG$=5)
PARAMETER(OXYGEN$=6)
PARAMETER(GASWAT$=7)
PARAMETER(LIQWAT$=8)
PARAMETER(LIBRW$=95)
PARAMETER(WAT67$=96)
PARAMETER(WATHF$=97)
PARAMETER(REALCO2$=98)
PARAMETER(WATER$=99)
PARAMETER(STEAM$=99)
PARAMETER(MASSFLOW$=99)
PARAMETER(SHAFT$=100)
PARAMETER(DSHAFT$=110)
PARAMETER(DSHAFTDT$=120)
PARAMETER(POWER$=200)
PARAMETER(HEAT$=300)
PARAMETER(DHEAT$=301)
PARAMETER(MOLEFRAC$=400)
PARAMETER(MASSFRAC$=450)
PARAMETER(PRESSURE$=500)
PARAMETER(ENTHALPY$=600)
PARAMETER(CONC$=650)
PARAMETER(MASS$=710)
PARAMETER(MASS1$=710)
PARAMETER(MASS2$=710)
PARAMETER(MASSDT$=720)
PARAMETER(MASS1DT$=720)
PARAMETER(MASS2DT$=720)
PARAMETER(MASSW$=730)
PARAMETER(MASS1W$=730)
PARAMETER(MASS2W$=730)
PARAMETER(INTENER$=810)
PARAMETER(INTENER1$=810)
PARAMETER(INTENER2$=810)
PARAMETER(INTENERDT$=820)
PARAMETER(INTENER1DT$=820)
PARAMETER(INTENER2DT$=820)
PARAMETER(INTENERW$=830)
PARAMETER(INTENER1W$=830)
PARAMETER(INTENER2W$=830)
PARAMETER(INTENERS$=850)
PARAMETER(INTENERSDT$=860)
PARAMETER(INTENERSW$=870)
PARAMETER(DYNVAR$=910)
PARAMETER(DYNVARDT$=920)
PARAMETER(DYNVARW$=930)
PARAMETER(ALGVAR$=980)
PARAMETER(CONTROL$=999)
PARAMETER(VARPAR$=1000)
C
INTEGER NODE1$,NODE2$,NODE3$,NODE4$,NODE5$,NODE6$,NODE7$,NODE8$
$ ,NODE9$,NODE10$
PARAMETER(NODE1$=-1)
PARAMETER(NODE2$=-2)
PARAMETER(NODE3$=-3)
PARAMETER(NODE4$=-4)
PARAMETER(NODE5$=-5)
PARAMETER(NODE6$=-6)
PARAMETER(NODE7$=-7)
PARAMETER(NODE8$=-8)
PARAMETER(NODE9$=-9)
PARAMETER(NODE10$=-10)
C
C Convenient shortcuts for identifiers in IVTAB
C
c$$$ INTEGER IVZD$,IVM1$,IVM2$,IVU1$,IVU2$,IVT1$,IVT2$,IVH1$,IVH2$
c$$$ $ ,IVS1$,IVS2$,IVX1$,IVX2$,IVP1$,IVP2$,IVUS$,IVTS$,IVDZDDT$
c$$$ $ ,IVDM1DT$,IVDM2DT$,IVDU1DT$,IVDU2DT$,IVDUSDT$,IVSIDE1$,
c$$$ $ IVSIDE2$
c$$$ PARAMETER(IVZD$=17)
c$$$ PARAMETER(IVM1$=18)
c$$$ PARAMETER(IVM2$=19)
c$$$ PARAMETER(IVU1$=20)
c$$$ PARAMETER(IVU2$=21)
c$$$ PARAMETER(IVT1$=22)
c$$$ PARAMETER(IVT2$=23)
c$$$ PARAMETER(IVH1$=24)
c$$$ PARAMETER(IVH2$=25)
c$$$ PARAMETER(IVS1$=26)
c$$$ PARAMETER(IVS2$=27)
c$$$ PARAMETER(IVX1$=28)
c$$$ PARAMETER(IVX2$=29)
c$$$ PARAMETER(IVP1$=30)
c$$$ PARAMETER(IVP2$=31)
c$$$ PARAMETER(IVUS$=32)
c$$$ PARAMETER(IVTS$=33)
c$$$ PARAMETER(IVDZDDT$=34)
c$$$ PARAMETER(IVDM1DT$=35)
c$$$ PARAMETER(IVDM2DT$=36)
c$$$ PARAMETER(IVDU1DT$=37)
c$$$ PARAMETER(IVDU2DT$=38)
c$$$ PARAMETER(IVDUSDT$=39)
INTEGER IVSIDE1$,IVSIDE2$
PARAMETER(IVSIDE1$=1)
PARAMETER(IVSIDE2$=2)
C
C Convenient shortcuts for identifiers in TITAB
C Numbers are half of the position in the TITAB array as this is
C used in the arrays INTAB, RATAB, IVTAB...
C BE 050221
INTEGER TIM$,TIP$,TIH$,TIE$,TIW$,TIQ$,TIT$,TIV$,TIS$,TIX$,TIU$,
$ TIY$,
$ TIX_J$,TIZC$,TIN$,TIZA$,TIY_J$,TIZD$,TIM1$,TIM2$,TIU1$,TIU2$
$ ,TIT1$,TIT2$,TIH1$,TIH2$ ,TIS1$,TIS2$,TIX1$,TIX2$,TIP1$,TIP2$
$ ,TIUS$,TITS$,TIDZDDT$ ,TIDM1DT$,TIDM2DT$,TIDU1DT$,TIDU2DT$
$ ,TIDUSDT$
PARAMETER(TIM$=1)
PARAMETER(TIP$=2)
PARAMETER(TIH$=3)
PARAMETER(TIE$=4)
PARAMETER(TIW$=5)
PARAMETER(TIQ$=6)
PARAMETER(TIT$=7)
PARAMETER(TIV$=8)
PARAMETER(TIS$=9)
PARAMETER(TIX$=10)
PARAMETER(TIU$=11)
PARAMETER(TIY$=12)
PARAMETER(TIX_J$=13)
PARAMETER(TIZC$=14)
PARAMETER(TIN$=15)
PARAMETER(TIZA$=16)
PARAMETER(TIY_J$=17)
PARAMETER(TIZD$=18)
PARAMETER(TIM1$=19)
PARAMETER(TIM2$=20)
PARAMETER(TIU1$=21)
PARAMETER(TIU2$=22)
PARAMETER(TIT1$=23)
PARAMETER(TIT2$=24)
PARAMETER(TIH1$=25)
PARAMETER(TIH2$=26)
PARAMETER(TIS1$=27)
PARAMETER(TIS2$=28)
PARAMETER(TIX1$=29)
PARAMETER(TIX2$=30)
PARAMETER(TIP1$=31)
PARAMETER(TIP2$=32)
PARAMETER(TIUS$=33)
PARAMETER(TITS$=34)
PARAMETER(TIDZDDT$=35)
PARAMETER(TIDM1DT$=36)
PARAMETER(TIDM2DT$=37)
PARAMETER(TIDU1DT$=38)
PARAMETER(TIDU2DT$=39)
PARAMETER(TIDUSDT$=40)
C
C Convenient shortcuts for identifiers in ANTAB
C BE 050221
C
INTEGER ANPAR$,ANRA$,ANIN$,ANIV$
PARAMETER(ANPAR$=1)
PARAMETER(ANRA$=2)
PARAMETER(ANIN$=3)
PARAMETER(ANIV$=4)
C=======================================================================
|
