Fortran编写反应动力学问题

18562317750

Fortran编写反应动力学时，整型变量与实型变量选择依据是什么，公共区选择哪些？？

forever999

统一选择成实型应该可以的吧。

jjzg2004

{:1106_362:}{:1106_382:}

18562317750

haha123456 发表于 2018-11-13 22:24
统一选择成实型应该可以的吧。

      非常感谢您的回复
IMPLICIT NONE
C
C     DECLARE VARIABLES USED IN DIMENSIONING
C
      INTEGER NSUBS, NINT,   NPO,   NIWORK, NWORK,
     +        NC,    NR,     NTCAT, NTSSAT, NCOMP,
     +        NRALL, NUSERV, NINTR, NREALR, NIWR,
     +        NWR
C
#iNClude "ppexec_user.cmn"
      EQUIVALENCE (RMISS, USER_RUMISS)
      EQUIVALENCE (IMISS, USER_IUMISS)
#iNClude "dms_NCOMP.cmn"
#iNClude "rxn_rcstrr.cmn"
#iNClude "rxn_rprops.cmn"
      EQUIVALENCE (TEMP,  RPROPS_UTEMP )
      EQUIVALENCE (PRES,  RPROPS_UPRES )
      EQUIVALENCE (VFRAC, RPROPS_UVFRAC)
      EQUIVALENCE (BETA,  RPROPS_UBETA )
      EQUIVALENCE (VVAP,  RPROPS_UVVAP )
      EQUIVALENCE (VLIQ,  RPROPS_UVLIQ )
      EQUIVALENCE (VLIQS, RPROPS_UVLIQS)
        EQUIVALENCE (B(1),  IB(1)        )
C
#iNClude "pputl_ppglob.cmn"
#iNClude "dms_maxWRt.cmn"
#iNClude "dms_plex.cmn"

C     DECLARE ARGUMENTS
C
      INTEGER IDXSUB(NSUBS), ITYPE(NSUBS),  INT(NINT),  IDS(2),
     +        NBOPST(6,NPO), IWORK(NIWORK), IDX(NCOMP), INTR(NINTR),
     +        IWR(NIWR),     NREAL,         KCALL,      KFAIL,   
     +        KFLASH,        NRL,           NRV,        I,   
     +        IMISS,         KDIAG,         KV, 纯组分汽相导热系数  KER错误返回代码,      
     +        DMS_IFCMNC偏移量,    LMW,           LMWI
C
      REAL*8  SOUT(1),       WORK(NWORK),   STOIC(NC,NSUBS,NR),   
     +        RATES(1),      FLUXM(1),      FLUXS(1),   RATCAT(NTCAT),
     +        RATSSA(NTSSAT),Y(NCOMP),      X(NCOMP),   X1(NCOMP),  
     +        X2(NCOMP)
C
      REAL*8  RATALL(NRALL), USERV(NUSERV), REALR(NREALR),
     +        WR(NWR),       RATEL(1),      RATEV(1),   XCURR,
     +        XMW(NCOMP_NCC),B(1),          TEMP,       PRES,     
     +        RGAS         
C
      REAL*8  REAL(NREAL),   RMISS, RGLOB_RMISS   XLEN,       DIAM,   rplg_rplugr.cmn
     +        VFRAC,         BETA,          VVAP,       VLIQ,   
     +        VLIQS,         VMXV,      DVMX 调用物性监视器    #include "rxn_rprops.cmn"
C
      REAL*8  NCARIN,        NCARGF,        NO2,        NCO,     
     +        NH2,           NCO2,          NH2O,       NCH4,      
     +        NN2,           NH2S,          NC6H6,      NTOTG,     
     +        NCARB,         NSULF,         MASH,       YO2,      
     +        YCO,           YH2,           YCO2,       YH2O,      
     +        YCH4,          YN2,           YH2S,       YC6H6,      
     +        YASH      
C     
      REAL*8  CCARB,         CO2,           CH2,        XC,        
     +        DP,            VOID,          RHOCOA,     VBED                       
C   
      REAL*8  PO2,           PCO,           PH2,        PCO2,     
     +        PH2O,          PCH4,          PN2,        PH2S,
     +        PC6H6,         PAMBI,         PT,         T,            
     +        PH2OEQ,        PCO2EQ,        PH2EQ        
C
        REAL*8  RCR,           Z,             FW,         E,   
     +        KFILM,         KASH,          KOVER,      K,         
     +        KH2O2,         RCARO2,        RCARH2O,    RCARCO2,      
     +        RCARH2,        RH2O2,         RCOH2O                          
这是煤气化的反应动力学例子，整型实型都有。

followdh

厉害 厉害 厉害