fortran浮点数异常

向着光亮 · 发表于 2018-5-16 20:45:29

编写了一个代码，通过光合有效辐射求日照时数，如图，但是总是显示NAN揭示浮点数异常，自己通过调试发现日落时角以后都是NaN，可是实在不知道原因，我把代码和工式发上来，希望大神指导！跪谢！

[Fortran] 纯文本查看 复制代码

program main
implicit none
real G, lat
real a,b
integer n
real,parameter::pi=3.1415926
real::ws(365),dr(365),angle(365),rizhao(365),q(365),qq(365),tianwen(365)
a=0.18
b=0.55
G=0.0820*24.0*60.0
lat=42.53*(pi/180.0)

open(10,file="D:\1998PAR.txt")
open(11,file="D:\1998.txt")
do n=1,365
    read(10,*)q(n)
    read(11,*)tianwen(n)
    q(n)=q(n)/0.5*0.0864
angle(n)=0.409*sin((2*pi/365.0)*n-1.39)
dr(n)=1+0.033*cos((2*pi/365.0)*n)
ws(n)=acos(-tan(lat)*tan(dr(n)))

qq(n)=24.0*(60.0/pi)*G*dr(n)*(ws(n)*sin(lat)*sin(dr(n))+cos(lat)*cos(dr(n))*sin(ws(n)))
rizhao(n)=(q(n)/qq(n)-a)/b*tianwen(n)
end do 
write(*,"(f8.2)")rizhao

stop
end