Fortran Coder

标题: fortran浮点数异常 [打印本页]

作者: 向着光亮 时间: 2018-5-16 20:45
标题: fortran浮点数异常
编写了一个代码，通过光合有效辐射求日照时数，如图，但是总是显示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

22.png (255.23 KB, 下载次数: 365)

33.png (24.23 KB, 下载次数: 368)

44.png (58.48 KB, 下载次数: 389)

1998.txt

3.56 KB, 下载次数: 3

1998PAR.txt

2.77 KB, 下载次数: 1

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

作者: 向着光亮 时间: 2018-5-17 15:52
什么意思啊，它在定义域，我实在找不到问题，您能帮忙看看不，文本我也上传了，

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

作者: 向着光亮 时间: 2018-5-25 08:34

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

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
如果你确信超过了定义域，并且代码逻辑是符合算法逻辑的，那么你应该质疑算法的正确性。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)