关于Chemkin III源程序中数组越界问题 - Fortran Coder

[Fortran] 纯文本查看 复制代码

 SUBROUTINE DDASAC (RES,NSYS,T,Y,YPRIME,TOUT,

* INFO,ISEN,RTOL,ATOL,IDID,SWORK,LSW,

* RWORK,LRW,IWORK,LIW,RPAR,IPAR,

* JAC,DRES,DFDYP)

IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER (I-N)

C*****END precision > double

C-------------------------------------------------------------------

C This code solves a system of differential/algebraic

C equations of the form f(t,y,yprime;p)=0. where p is a

C parameter vector of time independent constants.

C If requested the code will also solve the following system

C of equations for the parametric sensitivity functions W(t) :

C

C (df/dyprime)*W'(t)+(df/dy)*W(t)+df/dp=0.

C

C where W(t):=dy/dp

C--------------------------------------------------------------------

LOGICAL DONE

EXTERNAL RES,JAC,DRES,DFDYP

DOUBLE PRECISION SWORK

DIMENSION Y(*), YPRIME(*), INFO(15),ISEN(5)

DIMENSION RWORK(*), IWORK(*), SWORK(*),RTOL(*), ATOL(*)

DIMENSION RPAR(*), IPAR(*)

COMMON/DDA001/NPD,NTEMP,

* LML,LMU,LMXORD,LMTYPE,

* LNST,LNRE,LNJE,LETF,LCTF,LIPVT

COMMON/DDA002/INDEX,NALG,IDFDP,ICALC,NPAR

DATA LTSTOP,LHMAX,LH,LTN,

* LCJ,LCJOLD,LHOLD,LS,LROUND,

* LALPHA,LBETA,LGAMMA,

* LPSI,LSIGMA,LDELTA

* /1,2,3,4,

* 5,6,7,8,9,

* 11,17,23,

* 29,35,41/



C

C First see if parametric sensitivity calculations

C have been requested by the user.

C

ICALC=ISEN(1)

[Fortran] 纯文本查看 复制代码

  CALL DDASTP(TN,Y,YPRIME,NEQ,NSYS,

     *   RES,JAC,H,RWORK(LWT),INFO(1),IDID,RPAR,IPAR,

     *   RWORK(LPHI),RWORK(LDELTA),RWORK(LE),

     *   RWORK(LWM),IWORK(LIWM),

     *   RWORK(LALPHA),RWORK(LBETA),RWORK(LGAMMA),

     *   RWORK(LPSI),RWORK(LSIGMA),

     *   RWORK(LCJ),RWORK(LCJOLD),RWORK(LHOLD),

     *   RWORK(LS),HMIN,RWORK(LROUND),

     *   IWORK(LPHASE),IWORK(LJCALC),IWORK(LK),

     *   IWORK(LKOLD),IWORK(LNS),INFO(10),

     *   SWORK(LDTEM),SWORK(LEMAT),DRES,DFDYP)

[Fortran] 纯文本查看 复制代码

 SUBROUTINE DDASTP(X,Y,YPRIME,NEQ,NSYS,

     *  RES,JAC,H,WT,JSTART,IDID,RPAR,IPAR,

     *  PHI,DELTA,E,WM,IWM,

     *  ALPHA,BETA,GAMMA,PSI,SIGMA,

     *  CJ,CJOLD,HOLD,S,HMIN,UROUND,

     *  IPHASE,JCALC,K,KOLD,NS,NONNEG,

     *  DTEM,EMAT,DRES,DFDYP)

C-----------------------------------------------------------------

C     DDASTP solves a system of differential and algebraic

C     equations for one step, normally from t to t+h.

C-----------------------------------------------------------------

C*****precision > double

      IMPLICIT DOUBLE PRECISION(A-H,O-Z), INTEGER (I-N)

C*****END precision > double

      LOGICAL CONVGD

      DIMENSION Y(*), YPRIME(*), WT(*), PHI(NEQ,*), DELTA(*), E(*)

      DIMENSION DTEM(*), EMAT(*), WM(*), IWM(*)

      DIMENSION PSI(*), ALPHA(*), BETA(*), GAMMA(*), SIGMA(*)

      DIMENSION RPAR(*), IPAR(*)

      EXTERNAL RES,JAC,DRES,DFDYP

      COMMON/DDA001/NPD,NTEMP,

     *   LML,LMU,LMXORD,LMTYPE,

     *   LNST,LNRE,LNJE,LETF,LCTF,LIPVT

      COMMON/DDA002/INDEX,NALG,IDFDP,ICALC,NPAR

      DATA MAXIT/4/

      DATA XRATE/0.25E0/

      DATA ZERO/0.0E0/, PT25/0.25E0/, PT5/0.5E0/, PT9/0.9E0/