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

落月黎明 · 发表于 2019-12-1 11:47:23

请求帮助！最近为了能够使用chemkinIII开源软件，学习了Fortran语言，但遇到如下问题。所用为intel Fortran+IVS
forrtl: severe (408): fort: (3): Subscript #1 of the array SWORK has value -858993460 which is less than the lower bound of 1
调试过程中，发现LDTEM，LEMAT始终为-858993460.
看到别人的贴子，采用SAVE命令，该值变为0，就出现 Value 0 which is less than the lower bound of 1。
从初始为-858993460后，在运行过程中，这两个值始终未改变。在调试的时候我将他们的值赋为1，又出现其他问题。
由于该源程序是从网上下载，不清楚别人有没有调整过。请大家提提意见，谢谢。

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

01SUBROUTINE DDASAC (RES,NSYS,T,Y,YPRIME,TOUT,
02* INFO,ISEN,RTOL,ATOL,IDID,SWORK,LSW,
03* RWORK,LRW,IWORK,LIW,RPAR,IPAR,
04* JAC,DRES,DFDYP)
05IMPLICIT DOUBLE PRECISION (A-H,O-Z), INTEGER (I-N)
06C*****END precision > double
07C-------------------------------------------------------------------
08C This code solves a system of differential/algebraic
09C equations of the form f(t,y,yprime;p)=0. where p is a
10C parameter vector of time independent constants.
11C If requested the code will also solve the following system
12C of equations for the parametric sensitivity functions W(t) :
13C
14C (df/dyprime)*W'(t)+(df/dy)*W(t)+df/dp=0.
15C
16C where W(t):=dy/dp
17C--------------------------------------------------------------------
18LOGICAL DONE
19EXTERNAL RES,JAC,DRES,DFDYP
20DOUBLE PRECISION SWORK
21DIMENSION Y(*), YPRIME(*), INFO(15),ISEN(5)
22DIMENSION RWORK(*), IWORK(*), SWORK(*),RTOL(*), ATOL(*)
23DIMENSION RPAR(*), IPAR(*)
24COMMON/DDA001/NPD,NTEMP,
25* LML,LMU,LMXORD,LMTYPE,
26* LNST,LNRE,LNJE,LETF,LCTF,LIPVT
27COMMON/DDA002/INDEX,NALG,IDFDP,ICALC,NPAR
28DATA LTSTOP,LHMAX,LH,LTN,
29* LCJ,LCJOLD,LHOLD,LS,LROUND,
30* LALPHA,LBETA,LGAMMA,
31* LPSI,LSIGMA,LDELTA
32* /1,2,3,4,
33* 5,6,7,8,9,
34* 11,17,23,
35* 29,35,41/
36 
37C
38C First see if parametric sensitivity calculations
39C have been requested by the user.
40C
41ICALC=ISEN(1)

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

01CALL DDASTP(TN,Y,YPRIME,NEQ,NSYS,
02   *   RES,JAC,H,RWORK(LWT),INFO(1),IDID,RPAR,IPAR,
03   *   RWORK(LPHI),RWORK(LDELTA),RWORK(LE),
04   *   RWORK(LWM),IWORK(LIWM),
05   *   RWORK(LALPHA),RWORK(LBETA),RWORK(LGAMMA),
06   *   RWORK(LPSI),RWORK(LSIGMA),
07   *   RWORK(LCJ),RWORK(LCJOLD),RWORK(LHOLD),
08   *   RWORK(LS),HMIN,RWORK(LROUND),
09   *   IWORK(LPHASE),IWORK(LJCALC),IWORK(LK),
10   *   IWORK(LKOLD),IWORK(LNS),INFO(10),
11   *   SWORK(LDTEM),SWORK(LEMAT),DRES,DFDYP)

在调试后，发现此行出现问题。

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

01SUBROUTINE DDASTP(X,Y,YPRIME,NEQ,NSYS,
02     *  RES,JAC,H,WT,JSTART,IDID,RPAR,IPAR,
03     *  PHI,DELTA,E,WM,IWM,
04     *  ALPHA,BETA,GAMMA,PSI,SIGMA,
05     *  CJ,CJOLD,HOLD,S,HMIN,UROUND,
06     *  IPHASE,JCALC,K,KOLD,NS,NONNEG,
07     *  DTEM,EMAT,DRES,DFDYP)
08C-----------------------------------------------------------------
09C     DDASTP solves a system of differential and algebraic
10C     equations for one step, normally from t to t+h.
11C-----------------------------------------------------------------
12C*****precision > double
13      IMPLICIT DOUBLE PRECISION(A-H,O-Z), INTEGER (I-N)
14C*****END precision > double
15      LOGICAL CONVGD
16      DIMENSION Y(*), YPRIME(*), WT(*), PHI(NEQ,*), DELTA(*), E(*)
17      DIMENSION DTEM(*), EMAT(*), WM(*), IWM(*)
18      DIMENSION PSI(*), ALPHA(*), BETA(*), GAMMA(*), SIGMA(*)
19      DIMENSION RPAR(*), IPAR(*)
20      EXTERNAL RES,JAC,DRES,DFDYP
21      COMMON/DDA001/NPD,NTEMP,
22     *   LML,LMU,LMXORD,LMTYPE,
23     *   LNST,LNRE,LNJE,LETF,LCTF,LIPVT
24      COMMON/DDA002/INDEX,NALG,IDFDP,ICALC,NPAR
25      DATA MAXIT/4/
26      DATA XRATE/0.25E0/
27      DATA ZERO/0.0E0/, PT25/0.25E0/, PT5/0.5E0/, PT9/0.9E0/