这块该如何修改呢

时光深处 · 发表于 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

[求助] 这块该如何修改呢

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章