关于动力学FORTRAN的问题

qq593369792

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

复杂，看不懂啊

lsd_jc

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

xsahh

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

田旭明

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

zhangshengli

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

ken6666

看起來像打錯