| 查看: 1345 | 回复: 1 | |||
[交流]
我用“dimension"定义了数组，可是调试提示我没有声明此数组。
|
|
DIMENSION EELAS(6), EELASP(3),BBAR(6), BBARP(3), BBARN(3, 3), 1 DISTGR(3,3), GRF(3,3) PARAMETER(ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0, FOUR=4.D0, 1 SIX=6.D0) C---------------------------------------------------------------------------------------------- C------------------------------------------------------------------------------------------------ EMOD=PROPS(1) ENU=PROPS(2) LAME1=EMOD*ENU/((1+ENU)*(1-2*ENU)) LAME2=0.5*EMOD/(1+ENU) C---------------------------------------------------------------------------------------------- C JACOBIAN OR DET DET=DFGRD1(1, 1)*DFGRD1(2, 2)*DFGRD1(3, 3) 1 -DFGRD1(1, 2)*DFGRD1(2, 1)*DFGRD1(3, 3) IF(NSHR.EQ.3) THEN DET=DET+DFGRD1(1, 2)*DFGRD1(2, 3)*DFGRD1(3, 1) 1 +DFGRD1(1, 3)*DFGRD1(3, 2)*DFGRD1(2, 1) 2 -DFGRD1(1, 3)*DFGRD1(3,1)*DFGRD1(2, 2) 3 -DFGRD1(2, 3)*DFGRD1(3, 2)*DFGRD1(1, 1) END IF C------------------------------------------------------------------------------------------------- C CALCULATE THE GROWTH FACTOR IF TIME(1)=1 THEN STATEV(NTENS+1)=1.0005 GR0=STATEV(NTENS+1) GO TO 110 GR0=STATEV(NTENS+1) AB=DET/GR0**3 IF (AB.LE.1.0) THEN STATEV(NTENS+1)=GR0 GO TO 110 END IF C---------------------------------------------------------------------- EPS=1.0E-08 CALL DNEWT(GR0,EPS,RES,L) IF (L.NE.0) THEN WRITE(*,10) GR0 STATEV(NTENS+1) =GR0 END IF 10 FORMAT(001X,'GR0=',E15.6) END C----------------------------------------------------------------------- C CALCULATE GROTH TENSOR 110 DO K1=1,3 DO K2=1,3 GRF(K1,K2)=0 END DO END DO DO K3=1,3 GRF(K3,K3)=GR0 END DO C--------------------------------------------------------------------------- C CALCULATE ELASTIC TENSOR SCALE=ONE/GR0 DO K1=1, 3 DO K2=1, 3 DISTGR(K2,K1)=SCALE*DFGRD1(K2,K1) END DO END DO C---------------------------------------------------------------------------- C------------------------------------------------------------------------------ BBAR(1)=DISTGR(1, 1)**2+DISTGR(1, 2)**2+DISTGR(1, 3)**2 BBAR(2)=DISTGR(2, 1)**2+DISTGR(2, 2)**2+DISTGR(2, 3)**2 BBAR(3)=DISTGR(3, 3)**2+DISTGR(3, 1)**2+DISTGR(3, 2)**2 BBAR(4)=DISTGR(1, 1)*DISTGR(2, 1)+DISTGR(1, 2)*DISTGR(2, 2) 1 +DISTGR(1, 3)*DISTGR(2, 3) IF(NSHR.EQ.3) THEN BBAR(5)=DISTGR(1, 1)*DISTGR(3, 1)+DISTGR(1, 2)*DISTGR(3, 2) 1 +DISTGR(1, 3)*DISTGR(3, 3) BBAR(6)=DISTGR(2, 1)*DISTGR(3, 1)+DISTGR(2, 2)*DISTGR(3, 2) 1 +DISTGR(2, 3)*DISTGR(3, 3) 可是调试的时候老提示 1>G:\Fortran\Console22\Console22\Source1.for(70): error #6410: This name has not been declared as an array or a function. [GRF] 1>G:\Fortran\Console22\Console22\Source1.for(74): error #6410: This name has not been declared as an array or a function. [GRF] 1>G:\Fortran\Console22\Console22\Source1.for(81): error #6410: This name has not been declared as an array or a function. [DISTGR] 1>G:\Fortran\Console22\Console22\Source1.for(86): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(87): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(88): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(89): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(92): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(94): error #6410: This name has not been declared as an array or a function. [BBAR] 1>G:\Fortran\Console22\Console22\Source1.for(102): error #6410: This name has not been declared as an array or a function. [STRESS] 1>G:\Fortran\Console22\Console22\Source1.for(105): error #6410: This name has not been declared as an array or a function. [STRESS] |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
2楼2015-05-03 10:03:09