18562317750 发表于 2018-11-13 19:43:05

Fortran编写反应动力学问题

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

forever999 发表于 2018-11-13 19:43:05

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

jjzg2004 发表于 2018-11-13 19:43:05

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

18562317750 发表于 2018-11-13 19:43:05

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*8SOUT(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*8RATALL(NRALL), USERV(NUSERV), REALR(NREALR),
   +      WR(NWR),       RATEL(1),      RATEV(1),   XCURR,
   +      XMW(NCOMP_NCC),B(1),          TEMP,       PRES,   
   +      RGAS         
C
      REAL*8REAL(NREAL),   RMISS, RGLOB_RMISS   XLEN,       DIAM,   rplg_rplugr.cmn
   +      VFRAC,         BETA,          VVAP,       VLIQ,   
   +      VLIQS,         VMXV,      DVMX 调用物性监视器    #include "rxn_rprops.cmn"
C
      REAL*8NCARIN,      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*8CCARB,         CO2,         CH2,      XC,      
   +      DP,            VOID,          RHOCOA,   VBED                     
C   
      REAL*8PO2,         PCO,         PH2,      PCO2,   
   +      PH2O,          PCH4,          PN2,      PH2S,
   +      PC6H6,         PAMBI,         PT,         T,            
   +      PH2OEQ,      PCO2EQ,      PH2EQ      
C
        REAL*8RCR,         Z,             FW,         E,   
   +      KFILM,         KASH,          KOVER,      K,         
   +      KH2O2,         RCARO2,      RCARH2O,    RCARCO2,      
   +      RCARH2,      RH2O2,         RCOH2O                        
这是煤气化的反应动力学例子，整型实型都有。

followdh 发表于 2018-11-13 19:43:05

厉害 厉害 厉害
页: [1]
查看完整版本: Fortran编写反应动力学问题