上一篇博文写到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分钟），解决问题！