这个报错该怎么改呢，急求解答

时光深处

Do J=1,ILO
       Do I=M,1,-1
          Do K=1,N
                        LDONO = JC(J,I,K)
                        DO LDO = 1,LDONO
              if (NUTRINH(j,i,k) .le. 0.) NUTRINH(j,i,k)=eps
              if (NUTRINO3(j,i,k) .le. 0.) NUTRINO3(j,i,k)=eps
              if (NUTRIP(j,i,k) .le. 0.) NUTRIP(j,i,k)=eps
              if (DON(j,i,k) .le. 0.) DON(j,i,k)=eps
              if (DOP(j,i,k) .le. 0.) DOP(j,i,k)=eps
              if (PHYTOS(j,i,k) .le. 0.) PHYTOS(j,i,k)=eps
              if (PHYTOL(j,i,k) .le. 0.) PHYTOL(j,i,k)=eps
              if (ZOOS(j,i,k) .le. 0.) ZOOS(j,i,k)=eps
              if (ZOOL(j,i,k) .le. 0.) ZOOL(j,i,k)=eps
              if (POP(j,i,k) .le. 0.) POP(j,i,k)=eps
              if (PON(j,i,k) .le. 0.) PON(j,i,k)=eps
              if (bSi(j,i,k) .le. 0.) bSi(j,i,k)=eps
              if (BAC(j,i,k) .le. 0.) BAC(j,i,k)=eps
              if (cod1(j,i,k) .le. 0.) cod1(j,i,k)=eps                         !2.0
              if (cod2(j,i,k) .le. 0.) cod2(j,i,k)=eps                         !2.0
              if (cod3(j,i,k) .le. 0.) cod3(j,i,k)=eps                         !2.0
              if (codm(j,i,k) .le. 0.) codm(j,i,k)=eps


这几个的报错都是这样的，实在不知道问题出在哪了
Error: The shapes of the array expressions do not conform.   [NUTRINH]
            if (NUTRINH(j,i,k) .le. 0.) NUTRINH(j,i,k)=eps

时光深处

从定义上看是没什么问题的，NUTRINH定义的话是实测值，肯定不可能小于0，eps的定义是1.e-8，这样看没什么问题呀

时光深处

时光深处 发表于 2023-5-11 19:55
从定义上看是没什么问题的，NUTRINH定义的话是实测值，肯定不可能小于0，eps的定义是1.e-8，这样看没什么问 ...

C:\Users\谷芳芳\Pictures\Screenshots\屏幕截图 2023-05-13 172112.png

时光深处

楚香饭 发表于 2023-5-12 08:31
有必要的话，把定义部分代码给出来。甚至全屏截图

C:\Users\谷芳芳\Pictures\Screenshots\屏幕截图 2023-05-13 172112.png 我还想再请教下，这个报错下面有个三角，这个到底是什么作用呢？