qq593369792 发表于 2020-1-5 16:52:52

关于动力学FORTRAN的问题

C   KEROGEN FLOW RATE (KG/S) AND CONCENTRATION (KG/M**3 SHALE) IN ORIGINAL SHALE
        FKO      = 0.0192
      CKO      = 323.66-----------------为什么这变成了CK0而不是图片中的FK0
C-----------------------------------------------------------------------------------


C   RETRIEVE REACTION TEMPERATURE (K) AND LEFT KEROGEN FLOW RATE (KG/S)
        T      = SOUT(IDXSUB(1)-1+NCOMP_NCC+2)   -----------这是什么意思
        FK       = SOUT(IDXSUB(3)-1+1)                        ------------这是什么意思


C   RETRIEVE VOID FRACTION AND REACTOR VOLUME (M**3)
      VBED   = RCSTRR_VFRRC
      VOLR   = RCSTRR_VOLRC


C        RETRIVE MOLECULAR WEIGHT OF EACH COMPONENT (KG/KMOL)
        LMW      =DMS_IFCMNC('MW')
        DO I   =1,NCOMP_NCC
          LMWI   =LMW+I
          XMW(I) =B(LMWI)
        END DO


C   TOTAL PYROLYSIS RATE OF KERGOEN (KG KEROGEN/M**3 SHALE/S)
      K      = 6.9E10*EXP(-21790.0/T)
        RKEROGEN = K * CKO * (FK/FKO)**1.4 ---------------这的式子为什么和图片中的不一样


C   REACTION RATE OF EACH COMPONENT (CONVENTIONAL: KMOL/S; NONCONVENTIONAL: KG/S)
        RH2      =RKEROGEN * FACTH2    / XMW(3)* (1.0-VBED) * VOLR
        RH2O   =RKEROGEN * FACTH2O   / XMW(4)* (1.0-VBED) * VOLR
        RH2S   =RKEROGEN * FACTH2S   / XMW(6)* (1.0-VBED) * VOLR
        RNH3   =RKEROGEN * FACTNH3   / XMW(7)* (1.0-VBED) * VOLR
        RCO      =RKEROGEN * FACTCO    / XMW(9)* (1.0-VBED) * VOLR
        RCO2   =RKEROGEN * FACTCO2   / XMW(10) * (1.0-VBED) * VOLR
        RCH4   =RKEROGEN * FACTCH4   / XMW(11) * (1.0-VBED) * VOLR
        RC2H6    =RKEROGEN * FACTC2H6/ XMW(12) * (1.0-VBED) * VOLR
        RC3H8    =RKEROGEN * FACTC3H8/ XMW(13) * (1.0-VBED) * VOLR
        RC4H10   =RKEROGEN * FACTC4H10 / XMW(14) * (1.0-VBED) * VOLR
        ROIL   =RKEROGEN * FACTOIL   / XMW(15) * (1.0-VBED) * VOLR
        RCHAR    =RKEROGEN * FACTCHAR * (1.0-VBED) * VOLR
        RKEROGEN = -RKEROGEN * (1.0-VBED) * VOLR

C        WRITE(MAXWRT_MAXBUF(1),200) XMW(15)
C 200        FORMAT(1X,"XMW=",F11.5)
C        CALL DMS_WRTTRM(1)


C   INITIALIZE RATES
      DO 100 I = 1, NC
      RATES(I) = 0D0

qyf 发表于 2020-1-5 16:52:52

复杂，看不懂啊

lsd_jc 发表于 2020-1-5 16:52:52

T = SOUT(IDXSUB(1)-1+NCOMP_NCC+2) 是指温度变量的意思，下面那个不清楚。。

xsahh 发表于 2020-1-5 16:52:52

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

田旭明 发表于 2020-1-5 16:52:52

帮顶一下，Fortran这一块好难啊感觉

zhangshengli 发表于 2020-1-5 16:52:52

SOUT是两维数组，怎么括号里只有一个值？

ken6666 发表于 2020-1-5 16:52:52

看起來像打錯
页: [1]
查看完整版本: 关于动力学FORTRAN的问题