求帮看NaN出错原因

youngstar

RT lz小白，没学过fortran，最近要用个程序，程序是某本书上的算例，我只是改了其中的一些参数，可能和例子相差几个数量级，然后编译未报错，但输出结果出现了NaN,不懂错在哪里，代码在附件，求大神帮忙，感激不尽！

vvt

NaN 是由于数学错误，比如除数为零，或负数开方导致的。

至于具体原因，请自行单步调试。

如何单步调试可以看：
http://debug.w.fcode.cn
http://v.fcode.cn/video-debugger.html

youngstar

vvt 发表于 2017-3-13 21:37
NaN 是由于数学错误，比如除数为零，或负数开方导致的。

至于具体原因，请自行单步调试。

感谢！！！