Fortran Coder

标题: mona表示月份,这个循环如何修改 [打印本页]

作者: 时光深处    时间: 2023-9-1 18:20
标题: 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
[Fortran] 纯文本查看 复制代码
if(any(mona==[5,6,7,8]))then
end if





欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2