非稳态传热的时间步和迭代次数输出不对，请大家帮助

楚香饭

在我这里，你的代码执行到 Simpler.f 的 147 行
BLC=BLC+CON(I,J)+AIP(I,J)*F(I+1,J,N)+AIM(I,J)*F(I-1,J,N)
时，发生越界。
此时，N = NF = 11，而 F 定义为 F(NI,NJ,NFMAX), NFMAX=10。所以越界

另外，时间步长和迭代次数的输出不对，这是你自己才知道的说法，我不清楚什么是时间步长，什么是迭代次数。

你这代码 1000 多行，没有注释，别人很难看懂你的代码，你又不做任何介绍。

还有，什么是“输出不对”？应该是多少，实际是多少？提问时，一定要站在别人的角度想想，别人能否理解你的意思？