fortran热源程序指导调试

云一林

请问谁计算过大气热源？这里有一个原始热源程序，但是运行时零错误，零警告，运行结果如下：
Begin time =            1
end TIME=           1
Begin time =            2
end TIME=           2
Begin time =            3
end TIME=           3
Begin time =            4
end TIME=           4
Begin time =            5
end TIME=           5
Begin time =            6
end TIME=           6
Begin time =            7
end TIME=           7
Begin time =            8
end TIME=           8
Begin time =            9
end TIME=           9
Begin time =           10
end TIME=          10
Begin time =           11

forrtl: severe (36): attempt to access non-existent record, unit 15, file E:\nc
p\fortran\q1\air.grd
Image              PC        Routine            Line        Source
q1.exe             0041D579  Unknown               Unknown  Unknown

Incrementally linked image--PC correlation disabled.
Press any key to continue

总共12时长，11时次时出现问题，我看不出来是我输入数据问题还是输入数据的单位问题，还是程序哪里需要修改我没看出来，请问谁懂的或者计算过热源出来结果的能给我指导指导吗？我这里有30年的正确热源结果可以对照，程序附件已经上传，谁能帮我看看？

清宵逝水

从报错信息来看，是读取数据文件时出了问题，建议检查grd文件。
另外，所用程序在计算平流项时，差分格式写错了，需要修改