| 查看: 1245 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿985，本科211，0817化学工程与技术319求调剂
已经有4人回复
-
311求调剂
已经有6人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有4人回复
-
化学工程321分求调剂
已经有6人回复
-
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 菜鸟求助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人回复
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
- 专业: 金属非晶态、准晶和纳米晶
3楼2014-10-16 18:04:40
okagh
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 31
- 帖子: 18
- 在线: 24.9小时
- 虫号: 2819104
- 注册: 2013-11-22
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
4楼2014-10-16 21:25:16
