请教一简单问题

iemcsl

我用VFortran编的一有限元程序,其中一段如下:
        DO 150 IGAUS=1,NGAUS
        EXISP=POSGP(IGAUS)
C
C*** EVALUATE THE SHAPE FUNCTIOS AT THE SAMPLING POINTS
C
        CALL SFR2(DERIV,ETASP,EXISP,SHAPE)
C
C*** CALCULATE COMPNENTS OF THE EQUIVALENT NODAL LOADS
C
...................
        DVOLU=WEIGP(IGAUS)
........
   150 CONTINUE

我用DEBUG调试程序,发现在DO 150 循环中,开始IGAUS=1,等运行完CALL SFR2(DERIV,ETASP,EXISP,SHAPE)后,IGAUS=0,导致 DVOLU=WEIGP(IGAUS)语句中WEIGP( )数组维数超限. 子程序SFR2中根本没涉及到IGAUS, 这是什么原因?请问谁遇到过类似现象吗?谢谢!