关于程序不进入循环的问题

kyra

昨天我就忍着你，发代码截图干嘛？
第一，别人不能复制粘贴，想帮你调试代码，就得照着图片重新打一遍代码。试问几个人有这耐心？
第二，浪费论坛资源，浪费浏览者的流量。
第三，自己还麻烦。

为啥不直接复制粘贴呢？
最后，给你介绍一个很好的检查错误方法，叫单步调试：
http://v.fcode.cn/video-debugger.html
http://debug.w.fcode.cn

kyra

你试试对

        REAL::T(16,12)
        REAL::TA(16,12)
        REAL::ET(16,12)

这三个数组进行初始化

也就是

        REAL::T(16,12)  = 0
        REAL::TA(16,12)= 0
        REAL::ET(16,12)= 0

此外，你可以学习一下数组整体操作和数组片段，可以省很多代码，例如
DO N=1,7
    T(6,N)=T1
END DO
可以简写为
T(6,1:7) = T1

再例如
DO M=2,5
  DO N=1,11
      T(M,N)=10
  END DO
END DO
可以简写为
T(2:5,1:11) = 10

而

DO M=1,16
    DO N=1,12
        ET(M,N)=ABS(TA(M,N)-T(M,N))
    END DO
END DO

则可以简写为
ET = abs( TA - T )

kyra

你用的什么编译器？这是我修改后的输出。（我没仔细看对不对，反正有变化了）

kyra

我用的gfortran，想起来了，你的 err 是不是也没初始化？

kyra

err 不能初始化为0呀，不然肯定不进入循环。
第一次判断就没有误差了，就不需要进入循环计算了。

REAL::ERR=999.0 试试