I have been examining the implicit scheme for phosphorus in the Fennel model and it seems to me that there
is no constraint that ensures that the ratio of phosphorus uptake/nitrogen uptake is equal to the model's specified P/N ratio
(R_P2N). If this constraint is not enforced then, it seems to me that phosphorus will not necessarily be conserved (summing over
phytoplankton and dissolved phosphate). I suggest that the factor cff6 (used in the PO4 update: Bio(i,k,iPO4_)=Bio(i,k,iPO4_)/(1+cff6) )
should be calculated from the previously computed cff4 and cff5 as shown below.
the NO3, NH4, and PO4 update equations are:
Bio(i,k,iNO3_)=Bio(i,k,iNO3_)/(1+cff4)
Bio(i,k,iNH4_)=Bio(i,k,iNH4_)/(1+cff5)
Bio(i,k,iPO4_)=Bio(i,k,iPO4_)/(1+cff6)
the uptakes of NO3 and NH4 are:
NO3_uptake=Bio(i,k,iNO3_)*cff4/(1+cff4)
NH4_uptake=Bio(i,k,iNH4_)*cff5/(1+cff5)
the total nitrogen uptake (TNU) is given by: TNU=NO3_uptake+NH4_uptake
the phosphorus uptake should be equal to R_P2N*TNU:
Bio(i,k,iPO4_)*cff6/(1+cff6)=R_P2N*TNU.
solving for cff6 gives:
cff6=R_P2N*TNU/(Bio(i,k,iPO4_) - R_P2N*TNU),
where TNU=Bio(i,k,iNO3_)*cff4/(1+cff4) + Bio(i,k,iNH4_)*cff5/(1+cff5)
Does this make sense? Am I incorrect in my conclusion that the cff6 calculation in the code does not ensure the correct P/N ratio
in the nutrient uptakes?
thanks,
dave