Fortran Coder

查看: 312|回复: 2

[求助] 新人求助,很多语法错误,自己试了试,但是太菜了

[复制链接]

3

帖子

1

主题

0

精华

新人

F 币
14 元
贡献
5 点
发表于 2024-5-23 13:47:40 | 显示全部楼层 |阅读模式
本帖最后由 kyra 于 2024-5-23 13:52 编辑

详细的代码在附件,主要是一下几个问题
这是我的维度定义
real(kind = 8) :: F(3,3), B(3,3), J, t1, t2, t3, mu, kappa
     1 U(3,3),BBNEW(nblock,3,3),S(3,3),eigVec(nblock,3,3)
      dimension stressBE(nblock,ndir+nshr),stressT(nblock,ndir+nshr)

      dimension  STRESSA(nblock,6),STRESSB(nblock,6),DVSB(nblock,6)
      integer i
C
      CHARACTER*80 cmname




error #6410: This name has not been declared as an array or a function.   [U]                        U(J,K) = 1.0d0


DO J = 1,3
            DO K = 1,3
                  IF (J==K) THEN
                        U(J,K) = 1.0d0
                  ELSE
                        U(J,K) = 0.0D0
                  END IF
            END DO
      END DO


error #6351: The number of subscripts is incorrect.
BBNEW(1,1) = U(1,1) + dt*AA*(LAMADA-1+EE)**CC*MS*2*S(1,1)
error #6410: This name has not been declared as an array or a function   [BENEW]
BENEW(3,3) =error #6366: The shapes of the array expressions do not conform.   [STRESSBE]
      stressBE(i,1) = tt2 * S*(BENEW(1,1) -t1) + t3 B(3,3)-BBNEW(3,3)

vumat_nh11.for

9.69 KB, 下载次数: 1

258

帖子

0

主题

0

精华

版主

World Analyser

F 币
713 元
贡献
508 点

新人勋章美女勋章元老勋章热心勋章规矩勋章管理勋章

QQ
发表于 2024-5-23 13:56:03 | 显示全部楼层
本帖最后由 kyra 于 2024-5-23 13:59 编辑

real(kind = 8) :: F(3,3), B(3,3), J, t1, t2, t3, mu, kappa ,
     1 U(3,3),BBNEW(nblock,3,3),S(3,3),eigVec(nblock,3,3)

此处第一行结尾处加一个逗号。

BBNEW(nblock,3,3) 是一个三维数组,所以你不能
BBNEW(1,1) = U(1,1) + dt*AA*(LAMADA-1+EE)**CC*MS*2*S(1,1)  
只给2个下标。

其他问题,例如 BENEW 没有定义等。
建议你写上 implicit none,这些莫名其妙的问题就会变得很清晰。


3

帖子

1

主题

0

精华

新人

F 币
14 元
贡献
5 点
 楼主| 发表于 2024-5-23 14:18:10 | 显示全部楼层
kyra 发表于 2024-5-23 13:56
real(kind = 8) :: F(3,3), B(3,3), J, t1, t2, t3, mu, kappa ,
     1 U(3,3),BBNEW(nblock,3,3),S(3,3), ...


非常感谢您!我下去试试
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-6-25 01:31

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表