| 查看: 1161 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
okagh铁虫 (初入文坛)
|
[求助]
菜鸟求助关于fortran if循环的问题，经验丰富的前辈肯定一眼能看出问题出在哪
|
|
|
！ 我把源程序出问题的那部分发了上来，这个主要是一个if循环，clock在0~100和100~200区间程序一切正常debug也没问题，但是在200~300区间Ac的输出就不正常了，后来我debug模式调试，在200步之后发现循环是正常，但是在第十一行单步调试提示“无可用源”，程序运行到第十二行时查看局部变量中的blf ，看到blf的值经过十一行的运算没有改变，值还是上一个循环结尾的blf的值，因此在200~300区间程序把第十一行跳了过去，运行时又不显示出错，但结果是错的，请教大神能不能指出我的错误在哪？ 一 DO 1000 CLOCK=1,10000 二 DO 100 i=((CLOCK-1)*10+1),CLOCK*10 三 if ((CLOCK.LT.100 .AND. CLOCK.GT.0) .OR. (CLOCK.LT.300 .AND. CLOCK.GT.200)) then！ 四 DO 110 j=1,mm 五 Aa(j)=Aa0*Ac(j) 六 ok(j)=Ke0*Ac(j)/Ad(j+50) 七 f(j)=ok(j)*dt1/(Aa(j)*dx*dx) 八 g(j)=g0*Ad(j+50)/Ac(j) 九 alf(j)=g(j)**2*dt1*dt1/(Ab*Aa(j)*(1+g(j)*dt1/Ab))-1-2*f(j)-g(j)*dt1/Aa(j) 十 alf(j)=alf(j)/f(j) 十一 blf(j)=-Ac(j)-(W/Aa(j))*exp(-alpha*j*dx)*exp(-4*0.69315*((i-11000)/1000.0)**2)-g(j)*dt1*Ad(j+50)/Aa(j) ！出问题的地方 出问题的地方 出问题的地方 十二 blf(j)=blf(j)/f(j) 110 CONTINUE blf(1)=blf(1)-Ac(1) blf(mm)=blf(mm)-Ac(mm) DO 120 k=2,mm alf(k)=alf(k)-1/alf(k-1) blf(k)=blf(k)-blf(k-1)/alf(k-1) 120 CONTINUE Ac(mm)=blf(mm)/alf(mm) DO 130 k=mm-1,1,-1 Ac(k)=(blf(k)-Ac(k+1))/alf(k) 130 CONTINUE else DO 140 j=1,mm Aa(j)=Aa0*Ac(j) ok(j)=Ke0*Ac(j)/Ad(j+50) f(j)=ok(j)*dt1/(Aa(j)*dx*dx) g(j)=g0*Ad(j+50)/Ac(j) alf(j)=g(j)**2*dt1*dt1/(Ab*Aa(j)*(1+g(j)*dt1/Ab))-1-2*f(j)-g(j)*dt1/Aa(j) alf(j)=alf(j)/f(j) blf(j)=-Ac(j)-g(j)*dt1*Ad(j+50)/Aa(j) blf(j)=blf(j)/f(j) 140 CONTINUE blf(1)=blf(1)-Ac(1) blf(mm)=blf(mm)-Ac(mm) DO 150 k=2,mm alf(k)=alf(k)-1/alf(k-1) blf(k)=blf(k)-blf(k-1)/alf(k-1) 150 CONTINUE Ac(mm)=blf(mm)/alf(mm) DO 160 k=mm-1,1,-1 Ac(k)=(blf(k)-Ac(k+1))/alf(k) 160 CONTINUE end if 100 CONTINUE write(*,*)Ac(5) 1000 CONTINUE |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有18人回复
-
垃圾破二本职称评审标准
已经有12人回复
-
职称评审没过，求安慰
已经有17人回复
-
EST投稿状态问题
已经有7人回复
-
谈谈两天一夜的“延安行”
已经有15人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 菜鸟求助fortran数组声明及子程序的一些问题啊
已经有4人回复
- Fortran输出二进制数据，用matlab读出，结果出错，求助问题何在？
已经有4人回复
- 求助：运行Fortran程序 提示出现的问题?
已经有16人回复
- 菜鸟求助fortran编写二阶微分方程
已经有5人回复
- fortran循环求助
已经有3人回复
- 求助：intel visual fortran的格式问题
已经有14人回复
- 求助：win7装fortran软件需要注意哪些问题，有没有针对性的软件？
已经有6人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- 用fortran遇到了一个问题，来求助了，谢谢先
已经有3人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】一点FORTRAN的问题【已完结】
已经有3人回复
- 【求助】Fortran编译目标文件后的链接问题
已经有4人回复
okagh
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 31
- 帖子: 18
- 在线: 24.9小时
- 虫号: 2819104
- 注册: 2013-11-22
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
3楼2014-10-16 18:04:40
qchem
铁杆木虫 (著名写手)
- 应助: 280 (大学生)
- 金币: 6871.6
- 散金: 1411
- 红花: 43
- 帖子: 2586
- 在线: 669.2小时
- 虫号: 183011
- 注册: 2006-02-12
- 专业: 理论和计算化学
2楼2014-10-16 13:45:31
okagh
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 31
- 帖子: 18
- 在线: 24.9小时
- 虫号: 2819104
- 注册: 2013-11-22
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
4楼2014-10-16 21:25:16