fortran浮点数异常

查看数: 14427 | 评论数: 6 | 收藏 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2018-5-16 20:45

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

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

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

向着光亮 发表于 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] 范围内

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

向着光亮 发表于 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] 范围内

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

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

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

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

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

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

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

回复