| 查看: 495 | 回复: 1 | ||
[求助]
fortran 中DOwhile的失效问题 已有1人参与
|
|
我的程序中有这样一段 但是不知怎么回事，当Y0（2）大于等于1.0时它不进入DO WHILE的循环 很是郁闷 有知道的可以帮忙解决下吗 j=1 p=P_exp write(*,*) "p_exp = ", p err1=5.0e-2 err2=8.0e-2 yn(1)=P_exp yn(2)=y(1) i=0 35 j=j+1 r=1.0 call Adf(z,x,y,T,P,AA) write(*,*) AA(1,3), AA(2,3) call GAUSS(2,1,AA) y0(1)=yn(1)-r*AA(1,3) y0(2)=yn(2)-R*AA(2,3) c write(*,*)'y0(2)',y0(2) do 36, while(((Y0(1).le.0.0).or.(Y0(2).le.0.0).or.(Y0(2).ge.1.0)) & .and.(i.lt.10)) i=i+1 R = R/10.0 Y0(1) = yn(1)- R*AA(1,3) Y0(2) = yn(2)- R*AA(2,3) 36 continue k=0 do 37, while(Y0(2).ge.1.0) k=k+1 R0 = R/10.0 Y0(2) = yn(2)- R0*AA(2,3) 37 continue write(*,*)'y0(2)',y0(2) yn(1)=y0(1) yn(2)=y0(2) y(1)=yn(2) y(2)= 1-yn(2) P = yn(1) write(*,*) "i =", i, "p = ", P,AA(1,3)/yn(1),AA(2,3)/yn(2) if( (abs(AA(1,3)/yn(1)).gt.ERR1) & .or.(abs(AA(2,3)/yn(2)).gt.ERR2) ) goto 35 call AROL(z,x,T,RO_l,p) call AROV(z,y,T,RO_v,p) P_cal = P rou1_l=ro_l rou1_v=ro_v write(*,*) "P=",P_cal,"Y(1)=" ,y(1) write(*,*) "rou_cal = ",ro_l,ro_V |
回复此楼» 猜你喜欢
上海大学实验技术岗位非升即走
已经有8人回复
-
今年审到国自然15份，谈谈感受
已经有14人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：VASP编译问题
已经有10人回复
- Fortran如何读取不确定行数的数据
已经有3人回复
- vsap 5.2 环境和编译问题
已经有14人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2015-06-22 08:53:29
20