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

时光深处 · 发表于 2023-5-11 17:10:51

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