作者zfrankz (双栖动物)
看板Fortran
标题[请益] 帮忙debug P币500
时间Thu Apr 17 10:32:58 2014
不好意思 有人可以帮我把这个程式码debug吗?
我已经试过很多次都 都无法
我是用fortran95
如果OK P币500 但你要教我怎麽把P币给你
-----------程式码如下---------------------
Program Energy
IMPLICIT NONE
REAL::E1,E2,E3,a1,a2,a3,a4,b1,b2,b3,b4,c1,c2,c3,i,s,r
OPEN(Unit=11, file="Energy out")
Data E1,E2,E3/16000.0,0.1,0.001/
Data a1,a2,a3,a4/750.0,-0.047,0.0,-0.036/
Data b1,b2,b3,b4/0.00036,-3.612,-0.015,-5.0/
Data c1,c2,c3/1.0,-0.29,-100.0/
100 FORMAT(a4,2x,a15,2x,a15,2x,a15)
WRITE(11,100)"N","E1","E2","E3"
WRITE(11,100)"","kcal/m**2","kcal/m**2","kcal/m**2"
WRITE(11,200)0,E1,E2,E3
DO 300 i=1,2000,1
200 FORMAT(i4,2x,f15.3,2x,f15.5,2x,f15.6)
E1=E1+0.1*(a1+a2*E1+a3*E1**2+a4*E1*E2)
E2=E2+0.1*(b1*E1*E2+b2*E3+b3*E3**2+B4*E2*E3)
E3=E3+0.1*(C1*E2*E3+C2*E3+C3*E3**2)
r=real(i)/10
s=(r-aint(r))*10
IF (s.eq.0) WRITE(11,200) i/10,E1,E2,E3
300 continue
END Program Energy
thanks
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.45.29
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/Fortran/M.1397701981.A.1FC.html
1F：推 rex0707:i必须宣告成integer 04/17 10:50
2F：→ rex0707:或许你可以改用gfortran 他会告诉你语法哪里有错 04/17 10:50
3F：→ fcuiii:if(s .eq. 0.) 04/18 00:29
4F：→ fcuiii:10 要用浮点数10. 04/18 00:33
5F：→ egg9:i宣告成integer, 还有IF判断式你不能期望两个浮点相等 04/18 09:34