Phosphorus uptake in Fennel model's implicit scheme

[dullman]

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

[arnaudlaurent]

The formulations of nutrient uptake with the implicit scheme are detailed in the attached file below. You can see from the derivation of PO4_uptake that with the formulation Bio(i,k,iPO4_)=Bio(i,k,iPO4_)/(1+cff6), the phosphorus uptake is indeed equal to R_P2N*TNU.