-
Notifications
You must be signed in to change notification settings - Fork 6
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
new functions for P-adsorp from salinity/temperature #47
Comments
Hi Peisheng I think this is actually a whole new adsorption model, rather than a salinity correction to the existing model This needs to be added as a new PO4AdsorptionModel libaed-water/src/aed_phosphorus.F90 Line 314 in a27a543
The issue is also that the Pexch is changing dynamically and is not a constant at 0.78. |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
libaed-water/src/aed_phosphorus.F90
Line 308 in a27a543
Below is a new function for the Kpo4p prediction as function of salinity and temperature. This function comes out directly from Zhang and Huang 2011, and is much more accurate and reliable than the previous one.
Note this function calculate the Kpo4p directly, not a scale factor. So in the aed_phosphorus.f90 the line 308
Kpo4p = data%Kpo4p * Kpo4p_fT_fSal(data%theta_Kpo4, data%K_sal, salt, temp)
should be changed to
Kpo4p = Kpo4p_fT_fSal(data%Pexch, salt, temp)
The function only requirs one parameter, which is the exchangable phosphate content in the particles (Pexch).
The text was updated successfully, but these errors were encountered: