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

[Fortran] 纯文本查看 复制代码

if(any(mona==[5,6,7,8]))then
end if