Fortran Coder

标题: 请教一下Fortran程序中子程序参数为何发生变化 [打印本页]

作者: elf唯sj 时间: 2018-11-20 19:52
标题: 请教一下Fortran程序中子程序参数为何发生变化

GREASELINETHERMEHL.f90 (15.61 KB, 下载次数: 9)
这个程序同样出现nan的问题经过我单步调试发现子程序EROEQ在结束第二轮的第一次循环后，参数B1的值正常，为1点多，为什么进入第二轮第二次循环时，B1的值自动变为了0，导致后面出现nan，刚刚接触Fortran不久，请各位大神赐教，感激不尽。

作者: vvt 时间: 2018-11-21 08:20
天哪，这是基础知识呀。子程序局部变量，在子程序返回之后，就被回收了。
下一次调用子程序，局部变量成为“未初始化状态”

如果你想保留上一次的值，可以定义其具有 save 属性：
real , save :: B1

此外，变量都不定义，是一个不好的习惯。

作者: elf唯sj 时间: 2018-11-21 09:58

vvt 发表于 2018-11-21 08:20
天哪，这是基础知识呀。子程序局部变量，在子程序返回之后，就被回收了。
下一次调用子程序，局部变量成为 ...

刚刚接触Fortran不就好多基础知识还没有完全掌握谢谢您啦
但是我加上save以后，发现nan的现象并没有改善所以可能我之前的判断是错的请您能不能帮我看一下这个程序的问题出在哪里了十分感谢！

作者: vvt 时间: 2018-11-21 10:53
加上 save 以后，是否 B1 的值就正常了？
NaN 可能出在别处，继续调试就好了。

没输入文件没法判断。而且我也不知道你的程序正常逻辑应该是怎样的，所以这种很难帮助你。

作者: elf唯sj 时间: 2018-11-21 19:09

vvt 发表于 2018-11-21 10:53
加上 save 以后，是否 B1 的值就正常了？
NaN 可能出在别处，继续调试就好了。

我是真的不知道该怎么继续调试了

可能还没有掌握单步调试的要领请您赐教，应该怎么继续找问题的所在
另外您说的输入文件是什么我运行这个程序里面的参数已经设定好了这是程序最后是想得出在考虑温度的情况下，油膜膜厚和压力随坐标的变化
请您帮忙指导吧十分感谢

作者: vvt 时间: 2018-11-21 19:53
就是你运行程序以后该输入什么？

READ(*,'(A)')S
READ(*,*)W,US,FF
READ(*,*)NZ,CU

这些

作者: elf唯sj 时间: 2018-11-21 20:19

vvt 发表于 2018-11-21 19:53
就是你运行程序以后该输入什么？

READ(*,'(A)')S

这个不需要输入在上一步'Show the example or not (Y or N) 这里输入y就可以了

作者: vvt 时间: 2018-11-21 21:49
这程序你也别改了，错误太多太多了。各种变量都需要 save
比如
UCAL 函数里的 A1 A2 A3
THERM 函数里的 T
TCAL 函数里的 A1 A2 A3 A4 A5 A6
等等

如果你使用的 IVF，你可以在工程属性里这样设置：

作者: elf唯sj 时间: 2018-11-22 10:28

vvt 发表于 2018-11-21 21:49
这程序你也别改了，错误太多太多了。各种变量都需要 save
比如
UCAL 函数里的 A1 A2 A3

啊这些都要save啊可是我用的是codeblocks 没找到在哪里设置所有变量save
所以您调试以后就是这个变量save的问题吗

作者: elf唯sj 时间: 2018-11-22 10:33

vvt 发表于 2018-11-21 21:49
这程序你也别改了，错误太多太多了。各种变量都需要 save
比如
UCAL 函数里的 A1 A2 A3

可否跟您加一下qq或者微信之类的

作者: vvt 时间: 2018-11-22 12:21
本帖最后由 vvt 于 2018-11-22 12:23 编辑

你可以在每个子程序里都加上
save

比如
subroutine sub()
save
end subroutine sub

然后编译选项加上
-finit-local-zero
QQ微信就不用加了，你可以加本论坛QQ群：2338021

作者: elf唯sj 时间: 2018-11-24 08:23

vvt 发表于 2018-11-22 12:21
你可以在每个子程序里都加上
save

好的十分感谢！请问编译选项这个怎么加啊我不大明白

作者: vvt 时间: 2018-11-26 08:32
Projects 菜单，Build Options，Other Options

作者: elf唯sj 时间: 2018-11-27 18:18

vvt 发表于 2018-11-26 08:32
Projects 菜单，Build Options，Other Options

十分感谢！

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