fcode 发表于 2014-2-27 21:47
so，可喜可贺。感谢楼上的捐赠。

我是真觉得咱们这个论坛很不错，才会支持的。真心谢谢管理员、版主和各位大神。

so，可喜可贺。感谢楼上的捐赠。

fcode 发表于 2014-2-27 19:47
楼主你别这样，确定问题解决了再采纳。

To:珊瑚虫，这个程序里大多数是 NAG Library 里的，所以函数库有问 ...

的确是数据的问题，子程序都是ok的。我按照楼上的建议调整了下一个参数，程序确实可以运行了，而且得出的数据与之前的比对了一下，结果还是很理想的。。。

是的

楼主你别这样，确定问题解决了再采纳。

To:珊瑚虫，这个程序里大多数是 NAG Library 里的，所以函数库有问题的可能性很小。应该还是数据的事儿。

珊瑚虫 发表于 2014-2-27 18:43
试了一下，我的看法是
程序退出的地方在第387行（第一次调用C06LAF子程序的时候）：
必须同时满足：RELERX ...

好的，谢谢您的建议！

btx97 发表于 2014-2-27 16:21
我是来灌水的。　不学习fortran，都敢用得这么专业。　称赞下伟大的冒险精神。 ...

有时候真的迫不得已。。。唉。。。

我是来灌水的。　不学习fortran，都敢用得这么专业。　称赞下伟大的冒险精神。

试了一下，我的看法是
程序退出的地方在第387行（第一次调用C06LAF子程序的时候）：
必须同时满足：RELERX大于0.1D0*ERREST(NCURR)以及LINC不等于NCURR
我把上述4个量打印出来的结果是：
RELERX=  1.000000000000000E-003
0.1D0*ERREST(NCURR)=1.869595362724266E-003
LINC=0
NCURR=1
这说明你的第一个大于没有满足，至于为什么没有满足 你就顺着往上查找，可能是你的某个参数输入错误了，如果是用过的程序，本身应该没什么问题。