版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (3016)
>虫友互识 (732)
>文献求助 (167)
>休闲灌水 (134)
>导师招生 (133)
>考博 (124)
>教师之家 (80)
>找工作 (73)
>硕博家园 (72)
>博后之家 (71)
>论文投稿 (70)
>论文道贺祈福 (65)
>基金申请 (61)
>考研 (55)
>公派出国 (49)
>招聘信息布告栏 (31)

返回列表

jzq19841102

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1656.5
帖子: 36
在线: 47.9小时
虫号: 445771
注册: 2007-10-28
专业: 金属材料的微观结构

[求助] 请Fortran高手指点！！！

用Absoft70编了一个小程序，在赋值那两行一直报错，不知道怎么解决，请高手指点！！！谢谢！！！

程序如下：

Real::aa
aa=((-5.6409215)*(10**(-13))*(TEMPE**5)
   &+(0.01374899)*(EFEPSE**5)+(0.0038608)*
   &TEMPE*(EFEPSE**4)-(2.57162431)*(10**(-5))*
   &(TEMPE**2)*(EFEPSE**3)-(1.0327334)*(10**(-7))*
   &(TEMPE**3)*(EFEPSE**2)+(6.5978419)*(10**(-10))*
   &(TEMPE**4)*EFEPSE+(6.99708629)*(10**(-10))*
   &(TEMPE**4)-(3.776971)*(EFEPSE**4)+(0.0069343)*
   &TEMPE*(EFEPSE**3)+(0.00058326)*(TEMPE**2)*
   &(EFEPSE**2)-(1.534483)*(10**(-6))*(TEMPE**3)*
   &EFEPSE+(3.5831657)*(10**(-6))*(TEMPE**3)+
   &(16.316593)*(EFEPSE**3)-(0.79593988)*TEMPE*
   &(EFEPSE**2)+(0.000486652)*(TEMPE**2)*
   &EFEPSE-(0.00921802)*(TEMPE**2)+(318.481403)*
   &(EFEPSE**2)+(0.9722918)*TEMPE*EFEPSE+(7.7097759)*
   &TEMPE-(584.3833107)*EFEPSE-2210.448129)
IF(ISTATUS.EQ.1.AND.KELE.GT.0) THEN
      IF (TEPSE.GT.aa) THEN
      USRE2(1)=1.0-EXP(-((-4.24296403)*(10**(-12))*
   &(TEMPE**5)-(1.107242)*(EFEPSE**5)+(0.003646453)*
   &TEMPE*(EFEPSE**4)-(9.831506)*(10**(-5))*
   &(TEMPE**2)*(EFEPSE**3)+(1.301407)*(10**(-7))*
   &(TEMPE**3)*(EFEPSE**2)+(3.3214738)*(10**(-10))*
   &(TEMPE**4)*EFEPSE+(2.3691111)*(10**(-8))*(TEMPE**4)-
   &(8.739411)*(EFEPSE**4)+0.147491*TEMPE*(EFEPSE**3)+
   &0.0007152154*(TEMPE**2)*(EFEPSE**2)-(2.9709698)*
   &(10**(-6))*(TEMPE**3)*EFEPSE-(4.105068)*(10**(-5))*
   &(TEMPE**3)+131.416869*(EFEPSE**3)-(1.6836256)*TEMPE*
   &(EFEPSE**2)+0.00564391*(TEMPE**2)*EFEPSE+0.0248135*
   &(TEMPE**2)+647.23637*(EFEPSE**2)-(3.810692)*TEMPE*
   &EFEPSE+0.032757*TEMPE+821.969366*EFEPSE-3242.5681412)*
   &((TEPSE-((-5.6409215)*(10**(-13))*(TEMPE**5)+
   &(0.01374899)*(EFEPSE**5)+(0.0038608)*TEMPE*
   &(EFEPSE**4)-(2.57162431)*(10**(-5))*(TEMPE**2)*
   &(EFEPSE**3)-(1.0327334)*(10**(-7))*(TEMPE**3)*
   &(EFEPSE**2)+(6.5978419)*(10**(-10))*(TEMPE**4)*
   &EFEPSE+(6.99708629)*(10**(-10))*(TEMPE**4)-
   &(3.776971)*(EFEPSE**4)+(0.0069343)*TEMPE*(EFEPSE**3)+
   &(0.00058326)*(TEMPE**2)*(EFEPSE**2)-(1.534483)*
   &(10**(-6))*(TEMPE**3)*EFEPSE+(3.5831657)*(10**(-6))*
   &(TEMPE**3)+(16.316593)*(EFEPSE**3)-(0.79593988)*TEMPE*
   &(EFEPSE**2)+(0.000486652)*(TEMPE**2)*EFEPSE-
   &(0.00921802)*(TEMPE**2)+(318.481403)*(EFEPSE**2)+
   &(0.9722918)*TEMPE*EFEPSE+(7.7097759)*TEMPE-
   &(584.3833107)*EFEPSE-2210.448129))**
   &((-7.098369)*(10**(-9))*(TEMPE**4)+0.388622*
   &(EFEPSE**4)+0.001273*TEMPE*(EFEPSE**3)+
   &(2.773053)*(10**(-6))*(TEMPE**2)*(EFEPSE**2)+
   &(1.016019)*(10**(-6))*(TEMPE**3)*EFEPSE+
   &(-0.0195863)*TEMPE*(EFEPSE**2)+(-0.0029781)*
   &(TEMPE**2)*EFEPSE+(1.841733)*(10**(-5))*
   &(TEMPE**3)+(-5.613127)*(EFEPSE**3)+(-0.0134321)*
   &(TEMPE**2)+21.4240294*(EFEPSE**2)+2.923792*TEMPE*
   &EFEPSE+(-0.0032076)*TEMPE+(-962.367796)*EFEPSE+
   &2118.236633)))
      ELSE
      USRE2(1) = 0.0
      ENDIF
   ENDIF
   END

其中  USRE2  : USER ELEMENT STATE VARIABLES (OUTPUT: At the End of the STEP N)
      EFEPSE : EFFECTIVE STRAIN RATE
      TEMPE  : NODAL TEMPERATURES
      ISTATUS: 0 - the begain of the step
C                      1 - the end of the step
   KELE : ELEMENT NUMBER
   TEPSE  : TOTAL EFFECTIVE STRAIN

报错内容：cf90-324 error on line 1311 of .\def_usr.f: Assignment of a rank 1 expression to a rank 0 variable is not allowed.

cf90-324 error on line 1328 of .\def_usr.f: Assignment of a rank 1 expression to a rank 0 variable is not allowed.

报错的两行分别是：
aa=((-5.6409215)*(10**(-13))*(TEMPE**5)
USRE2(1)=1.0-EXP(-((-4.24296403)*(10**(-12))*

谢谢啦！！！

注：附件里有txt版程序

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 新建文本文档.txt

2012-08-10 17:19:53, 5.05 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

急求fortran运行错误原因，在线等已经有7人回复
如何改变fortran的浮点精度已经有28人回复
origin 曲线拟合求助，高手请进已经有11人回复
Fortran中关于数组的变化已经有5人回复
【求助】如何解这样的方程组【有效期至2011年4月30日】已经有9人回复
【求助】如何写FORTRAN程序实现求平均最近邻距离已经有34人回复
【求助】intel fortran for linux的问题已经有4人回复
【讨论】求助cpmd的编译问题已经有4人回复
【求助】求助，如何使linux能运行windows dos 的程序！已经有5人回复
【交流】Fortran语言答疑专帖已经有157人回复

1楼 2012-08-10 17:20:01

已阅关注TA 给TA发消息送TA红花 TA的回帖

smswt

铁杆木虫 (小有名气)

应助: 4 (幼儿园)
金币: 6049.6
帖子: 112
在线: 362.5小时
虫号: 905964
注册: 2009-11-18
专业: 光学

【答案】应助回帖

感谢参与，应助指数 +1

错误说出来了嘛
等号两边的rank（维度？）不一样

赞一下

回复此楼

2楼2012-08-10 22:08:26

已阅关注TA 给TA发消息送TA红花 TA的回帖

hitbenxun

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 2238.7
红花: 3
帖子: 233
在线: 63.8小时
虫号: 1437190
注册: 2011-10-11
专业: 工程热物理相关交叉领域

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢 2012-08-12 21:00:37

1、检查括号是否成对
2、检查变量是否都是纯变量，是否有数组混入其中

赞一下

回复此楼

3楼2012-08-10 22:22:33

已阅关注TA 给TA发消息送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

应助: 125 (高中生)
金币: 3820.1
散金: 5
红花: 9
帖子: 262
在线: 281.8小时
虫号: 1781694
注册: 2012-04-27
专业: 理论和计算化学

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2012-08-12 21:00:27

程序写出这样debug也难！

aa=((-5.6409215)*(10**(-13))*(TEMPE**5)
常数不要写成这样，(-5.6409215)*10**(-13)可以直接写成 -5.6409215E-13

USRE2(1)=1.0-EXP(-((-4.24296403)*(10**(-12))*
USRE2是一个数组，你没有先定义怎么可以使用？

赞一下

回复此楼

4楼2012-08-11 03:51:20

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 jzq19841102 的主题更新

返回列表