程序运行有问题，有关子程序C06LAF

zjl88915 · 发表于 2014-2-27 11:50:08

chuxf 发表于 2014-2-26 20:39
发你修改后的程序。这个错误就比较简单一些了

您好，昨晚学校断网，没办法给您回复了。多谢帮助，附件是我修改后的文件，请您再帮我看一下吧。

楚香饭 · 发表于 2014-2-27 12:14:49

本帖最后由 chuxf 于 2014-2-27 12:16 编辑

嗯，这个错误就很明了了。

Ua,Uw,Set 三个函数中，没有对 m1a,m2w,m2a 的定义，于是这三个量默认是整型的。

而我猜，你是想让他们是双精度。因为计算时他们的值小于1，整型的话就=0了，而后被 Cw=(m1w-m2w)/m2w 作为分母，出现了除0错误。

在这三个函数中，添加对 m1a,m2w,m2a 的定义。

[Fortran] 纯文本查看 复制代码

double precision m1a,m2w,m2a

得到如下计算结果，不知是否正确：
------------------------------------------------------------------
t(s) Pa(kPa) Pw(kPa) Set(cm)
------------------------------------------------------------------
0.1000E+05 0.43307 0.63307 0.00000

zjl88915 · 发表于 2014-2-27 13:08:57

chuxf 发表于 2014-2-27 12:14
嗯，这个错误就很明了了。

Ua,Uw,Set 三个函数中，没有对 m1a,m2w,m2a 的定义，于是这三个量默认是整型的 ...

太感谢您了！就是我想要的结果。您的指导真的帮了我很大的忙。

楚香饭 · 发表于 2014-2-27 13:17:00

zjl88915 发表于 2014-2-27 13:08
太感谢您了！就是我想要的结果。您的指导真的帮了我很大的忙。

不客气，常来论坛看看。

zjl88915 · 发表于 2014-2-27 14:30:41

chuxf 发表于 2014-2-27 13:17
不客气，常来论坛看看。

我按照你的指导修改后，可以运行了。但是当我读取1D_Unsat.dat这个文件中的T(1)这个数时，当T（1）<1e7时，很正常，一旦大于就又出现以前的问题了（正常的话我需要算到T(1)=1e11)。我按照C06LAF子程序的帮助文件的错误提示进行修改，问题依然存在。现在附件里是程序和问题截图，有时间您再帮我看下？我是小白。。。但是这个程序我需要尽快弄出来，好无力的赶脚。。。

楚香饭 · 发表于 2014-2-27 14:54:27

这个恐怕就是管理在 4 楼说的那种情况了。

T(1) 增大以后，要达到相同的精度也就越困难，你需要同步调整其他的参数，例如 RELERR，TFAC，ALPHAB（根据你的具体情况）

zjl88915 · 发表于 2014-2-27 14:59:34

chuxf 发表于 2014-2-27 14:54
这个恐怕就是管理在 4 楼说的那种情况了。

T(1) 增大以后，要达到相同的精度也就越困难，你需要同步调整其 ...

恩。。。那我只有再看看吧。。。但是我觉得问题不一定就是出在提示的地方。因为我这里类似这个计算的其他几个程序我计算时即使T很大也没有出现这个问题。唉。。。

zjl88915 · 发表于 2014-2-27 16:23:44

chuxf 发表于 2014-2-27 14:54
这个恐怕就是管理在 4 楼说的那种情况了。

T(1) 增大以后，要达到相同的精度也就越困难，你需要同步调整其 ...

您好，在您方便的时候，可不可以帮我再看一下？附件里是我另外的一个程序以及C06LAF的帮助文件，那个运行没有任何问题，两个程序的代码差不多的。。。只是输入的数有个别不同。

schoolhui · 发表于 2014-9-23 15:34:52

你好，请问你有完整的NAG fortran library 源码没有？

zjl88915 · 发表于 2014-9-29 13:44:24

schoolhui 发表于 2014-9-23 15:34
你好，请问你有完整的NAG fortran library 源码没有？

不好意思，这个真没有。。。

[子程序] 程序运行有问题，有关子程序C06LAF

评分

爱心勋章

新人勋章

水王勋章

元老勋章

热心勋章