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

fortran浮点数异常

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

fcode 发表于 2018-5-17 09:09:36

注意第 21 行
ws(n)=acos(-tan(lat)*tan(dr(n)))
此处的 acos 很可能不在定义域 [-1,1] 之间，这里你确认一下。

向着光亮 发表于 2018-5-17 15:52:17

什么意思啊，它在定义域，我实在找不到问题，您能帮忙看看不，文本我也上传了，

fcode 发表于 2018-5-18 07:48:41

ws(n)=acos(-tan(lat)*tan(dr(n)))
也就是
-tan(lat)*tan(dr(n)) 不在 [-1,1] 范围内

向着光亮 发表于 2018-5-25 08:34:38

fcode 发表于 2018-5-18 07:48
ws(n)=acos(-tan(lat)*tan(dr(n)))
也就是
-tan(lat)*tan(dr(n)) 不在 [-1,1] 范围内

好像不在，但是算出来确定是这样，您能帮忙看看吗

向着光亮 发表于 2018-5-25 11:15:56

fcode 发表于 2018-5-18 07:48
ws(n)=acos(-tan(lat)*tan(dr(n)))
也就是
-tan(lat)*tan(dr(n)) 不在 [-1,1] 范围内

是超过了定义域的范围，但是它算出来就是这样

kyra 发表于 2018-5-31 08:20:17

如果你确信超过了定义域，并且代码逻辑是符合算法逻辑的，那么你应该质疑算法的正确性。
页: [1]
查看完整版本: fortran浮点数异常