时光深处 发表于 2023-8-19 00:00:08

这块该如何修改呢

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:17

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
页: [1]
查看完整版本: 这块该如何修改呢