Fortran Coder

标题: 这块该如何修改呢 [打印本页]

作者: 时光深处 时间: 2023-8-19 00:00
标题: 这块该如何修改呢
if  (FNH(j,i,k)>200)  sPS(j,i,k)=0

   if (max(FNO3(j,i,k),FNH(j,i,k),Fdon(j,i,k))=FNH(j,i,k)) then
            sPS(j,i,k)=KpptgPS1 *exp(ktpptg*Ft(j,i,k))*
   *          Fpar(j,i,k)/parms * exp(1-Fpar(j,i,k)/parms)*
   *          min(SN(j,i,k),SP(j,i,k))
   else (max(FNO3(j,i,k),FNH(j,i,k),Fdon(j,i,k))=Fdon(j,i,k))
         sPS(j,i,k)= KpptgPS2 *exp(ktpptg*Ft(j,i,k))*
   *          Fpar(j,i,k)/parms * exp(1-Fpar(j,i,k)/parms)*
   *          min(SN(j,i,k),SP(j,i,k)
endif

作者: fcode 时间: 2023-8-19 09:44
if  (FNH(j,i,k)>200)  sPS(j,i,k)=0

   if (max(FNO3(j,i,k),FNH(j,i,k),Fdon(j,i,k))==FNH(j,i,k)) then
            sPS(j,i,k)=KpptgPS1 *exp(ktpptg*Ft(j,i,k))*
   *             Fpar(j,i,k)/parms * exp(1-Fpar(j,i,k)/parms)*
   *          min(SN(j,i,k),SP(j,i,k))
   else if (max(FNO3(j,i,k),FNH(j,i,k),Fdon(j,i,k))==Fdon(j,i,k)) then
                     sPS(j,i,k)= KpptgPS2 *exp(ktpptg*Ft(j,i,k))*
   *                Fpar(j,i,k)/parms * exp(1-Fpar(j,i,k)/parms)*
   *          min(SN(j,i,k),SP(j,i,k))
endif

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)