| 查看: 1387 | 回复: 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] |
» 猜你喜欢
Sci. Bull. 悲剧经验
已经有4人回复
找博士生导师
已经有7人回复
上海大学实验技术岗位非升即走
已经有11人回复
评审有感
已经有15人回复
26/27申博自荐-锂/钠电池方向
已经有4人回复
同样的基金本子,换个专家直接从C变A!
已经有3人回复
别被青基扩招骗了!26年科研内卷才刚刚开始
已经有4人回复
26/27博士推荐
已经有4人回复
2026博士还有哪些学校有名额
已经有8人回复
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有4人回复


2楼2015-05-03 10:03:09












回复此楼
30