这块该如何修改呢

fcode

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