| 查看: 1548 | 回复: 17 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
再问一段fortran代码
|
||
|
关于一段fortran的子例程子程序 SUBROUTINE UCP(DECRA,DESWA,STATEV,SERD,EC0,ESW0,P,QTILD, 1 TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND, 2 COORDS,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) ....... ....... ....... GOTO 100 ............ ......... 100 RETURN END 请问大家，这段程序中，goto语句 转到 标号100位置后，执行 returan语句， 程序返回开始位置 还是结束呢？？谢谢大家 |
回复此楼» 猜你喜欢
博后基金刷到的BUG,图片来的更直观
已经有9人回复
-
去英国的小伙伴儿都在哪儿租的房子呀？
已经有8人回复
-
太卷了
已经有22人回复
-
这个博士要读吗
已经有17人回复
-
论文已接收，但发现修改稿传的是旧版该怎么办？
已经有6人回复
-
有人中过人文社科类的博后特助吗？
已经有4人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有8人回复
-
最新消息：2024国自然下载文件名变了
已经有28人回复
-
审稿问题：为什么荧光激发波长和紫外吸收波长差的大？
已经有10人回复
-
初始合伙人来啦！（生物试剂耗材标准品）
已经有21人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran里一个子程序不运行是怎么回事
已经有13人回复
- fluent UDF代码调用预定义程序的问题（如fortran程序）
已经有15人回复
- fortran读写问题
已经有6人回复
- fortran输出求助
已经有6人回复
- 一百个FORTRAN任务怎么能一次提交SSH计算
已经有3人回复
- 想问一下影响fortran计算速度的因素都有哪些？
已经有12人回复
- 在Linux下的Fortran程序计算问题
已经有3人回复
- 急求fortran运行错误原因，在线等
已经有7人回复
- 请教 fortran 运行错误的原因
已经有13人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- Fortran的格式化输入输出问题
已经有14人回复
- fortran两个问题，文件中的空格，屏幕同一行覆盖输出
已经有23人回复
- Fortran被科学家抛弃了？
已经有49人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】fortran程序问题
已经有12人回复
- 【求助】初学fortran，
已经有15人回复
- 【转帖】DLA模型代码（VB，FORTRAN）
已经有7人回复
- 【求助】这样的公式，在fortran里该怎么写了，谢谢【求助成功】
已经有4人回复
- 【交流】Fortran语言答疑专帖
已经有157人回复
5楼2013-06-16 13:21:04
6楼2013-06-16 13:21:11
7楼2013-06-16 13:22:03
|
把源代码给您看看，这里goto语句用的真的不好吗？ SUBROUTINE UCP(DECRA,DESWA,STATEV,SERD,EC0,ESW0,P,QTILD, 1 TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND, 2 COORDS,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 CMNAME C DIMENSION DECRA(5),DESWA(5),STATEV(*),PREDEF(*),DPRED(*), 1 TIME(2),COORDS(*) C C DEFINE CONSTANTS C AA=0.009 DD=70.0 AIWA=0.0 BB=0.001 QQ=10000.0 RT=8.314 XX_sat=8.0 miu=7.0 omig=0.015 NN=1.2 YY=0.1 C c set initial RR, XX c IF(KSTEP.EQ.1.AND.KINC.EQ.1) THEN STATEV(4)=2.0 STATEV(3)=0.0 ENDIF IF(KSTEP.GE.2.AND.KINC.EQ.1) THEN STATEV(4)=STATEV(4)+YY/STATEV(3) ENDIF IF(KINC.EQ.1) THEN STATEV(1)=STATEV(4) STATEV(2)=0.0 ENDIF c c set FLAG_SIGN (elastic or inelastic) c IF(KINC.EQ.1) THEN FLAG_SIGN=0 ENDIF c RR=STATEV(1) XX=STATEV(2) c KK=1.5**0.5*(QTILD-AIWA)-RR IF(KK.LT.0) THEN DECRA(1)=0 DECRA(5)=0 FLAG_SIGN=1 GOTO 100 ELSE SV=KK DFSV=1.5**0.5 END IF c IF(FLAG_SIGN.EQ.0) THEN GOTO 100 ENDIF c c calculat inelastic strain rate c DECRA(1)=1.5**0.5*AA*(SV/DD)**NN*EXP(BB*(SV/DD)**(NN+1)) 1 *EXP(-QQ/(RT*TEMP))*DTIME IF(LEXIMP.EQ.1) THEN DECRA(5)=1.5**0.5*AA*NN*(SV/DD)**(NN-1)*(1/DD)*DFSV 2 *EXP(BB*(SV/DD)**(NN+1))*EXP(-QQ/RT*TEMP)*DTIME 3 +1.5**0.5*AA*(SV/DD)**NN*EXP(BB*(SV/DD)**(NN+1)) 4 *BB*(NN+1)*(SV/DD)**NN 5 *(1/DD)*DFSV*EXP(-QQ/(RT*TEMP))*DTIME END IF c c calculat XX c CR_rate=DECRA(1)/DTIME XX_rate=miu*CR_rate*(XX_sat-XX) XX=XX+XX_rate*DTIME STATEV(2)=XX STATEV(3)=STATEV(3)+XX_rate*DTIME STATEV(1)=STATEV(4)+XX_rate*(1-omig) c 100 RETURN END |
13楼2013-06-16 18:16:52
15楼2013-06-16 21:51:12
17楼2013-06-17 10:55:10
8