||
上一篇博文写到sh_bcfit存在错误,原因如下:
因此修改/opt/gamit10.71/gamit/lib/thdred.f源代码如下:
1)增加变量声明
2)参照写入时候的格式,读取按照相同样式;读取完成以后,把时间相应赋值给对应变量;
3)修改完成以后,重新编译gamit代码(最简单办法,重新安装)
但是发现问题还是无法解决,总是出现month = 0的情况,为了找到原因,在代码中添加了输出变量:
发现im0,imf都是0,说明这种读取方式不对,还是改回
read(iut) head,ite,tee,itb,tbb,itstp,tstp,sdelt,nepcht,ut1s,pole
采用数组方式读取,正确读取年月日信息;但是发现了另一个问题,写入和读取的时候,nics数值改变了!
NICS = number of initial conditions plus non-grav parameters
从中可以看出老版本T文件需要设置nics=6,本文使用10.71最新版本,应该不需要修改nics;注释掉这一句代码会后,重新编译(又要等待20分钟),解决问题!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 17:16
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社