时光深处 发表于 2023-9-1 18:20:55

mona表示月份，这个循环如何修改

if (mona=5,6,7,8) then
                BSPHYTOS(j,i,k)=sPS(j,i,k)-rPHY(j,i,k)*
   *                  FpptPS(j,i,k)-mPS(j,i,k)*FpptPS(j,i,k)*FpptPS(j,i,k)
   *                -dPS(j,i,k)*FpptPS(j,i,k)-gPS(j,i,k)*FzptS(j,i,k)
                BSPHYTOL(j,i,k)=sPL(j,i,k)-rPHY(j,i,k)*
   *                FpptPL(j,i,k)-mPL(j,i,k)*FpptPL(j,i,k)*FpptPL(j,i,k)
   *                -dPL(j,i,k)*FpptPL(j,i,k)-gPL(j,i,k)*FzptL(j,i,k)
              else
                BSPHYTOS(j,i,k)=sPS(j,i,k)*FpptPS(j,i,k)-rPHY(j,i,k)*
   *                  FpptPS(j,i,k)-mPS(j,i,k)*FpptPS(j,i,k)*FpptPS(j,i,k)
   *                -dPS(j,i,k)*FpptPS(j,i,k)-gPS(j,i,k)*FzptS(j,i,k)
                BSPHYTOL(j,i,k)=sPL(j,i,k)*FpptPL(j,i,k)-rPHY(j,i,k)*
   *                FpptPL(j,i,k)-mPL(j,i,k)*FpptPL(j,i,k)*FpptPL(j,i,k)
   *                -dPL(j,i,k)*FpptPL(j,i,k)-gPL(j,i,k)*FzptL(j,i,k)
              endif

Transpose 发表于 2023-9-1 18:59:37

if(any(mona==))then
end if
页: [1]
查看完整版本: mona表示月份，这个循环如何修改