This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Next revision Both sides next revision | ||
start:hype_model_description:hype_np_soil [2018/08/10 16:33] cpers [Balance SP – PartP] |
start:hype_model_description:hype_np_soil [2018/08/10 16:37] cpers [Balance SP – PartP] |
||
---|---|---|---|
Line 334: | Line 334: | ||
ENDDO | ENDDO | ||
- | The coefficients //Kfr//, and //Nfr// are soil dependent, bulkdensity is the dry density of soil (= 1300 kg/m3) and soilthick is the soil layer thickness. //fxn// is function value for the concentration //xn//, //fprimx// is the function value of the derivative of //fxn//, and <m>xn_{1}</m> is the concentration of the next iteration. The resulting SP equilibrium concentration (//xn=xn_1//) is used to calculate the PP eqiuilibrium concentration (<m> PPequi{_}conc</m>). | + | The coefficients //Kfr//, and //Nfr// are soil dependent, bulkdensity is the dry density of soil (= 1300 kg/m3) and soilthick is the soil layer thickness. //fxn// is function value for the concentration //xn//, //fprimx// is the function value of the derivative of //fxn//, and <m>xn_1</m> is the concentration of the next iteration. The resulting SP equilibrium concentration (<m>xn=xn_1</m>) is used to calculate the PP eqiuilibrium concentration (<m>PPequi_{conc}</m>). |
- | A new SP concentration (//SP_conc//) is calculated from the previous concentration (//SP_conc_old//) and the equilibrium concentration, taking into account the balance is not reached during the time step. | + | A new SP concentration (<m>SP_conc</m>) is calculated from the previous concentration (<m>SP_conc_old</m>) and the equilibrium concentration, taking into account the balance is not reached during the time step. |
- | <m> PPequi{_}conc = Kfr * ({xn_{1}}^{Nfr}) </m> | + | <m> PPequi_{conc} = Kfr * ({xn_1}^{Nfr}) </m> |
- | <m> Adsdes = (PPequi{_}conc - concsolid)*(1-e^{-Kadsdes}) </m> | + | <m> Adsdes = (PPequi_{conc} - concsolid)*(1-e^{-Kadsdes}) </m> |
- | <m> SP{_}conc = SP{_}conc{_}old - {(Adsdes * bulkdensity * LayerThick)/ vol} </m> | + | <m> SP_conc = SP_conc_old - {(Adsdes * bulkdensity * LayerThick)/ vol} </m> |
//Adsdes// is the amount of adsorbed/desorbed P, //Kadsdes// is a soil dependent rate constant, //equi_concsolid// is the solid phase equilibrium concentration, and //concsolid// is the initial concentration the solid phase. | //Adsdes// is the amount of adsorbed/desorbed P, //Kadsdes// is a soil dependent rate constant, //equi_concsolid// is the solid phase equilibrium concentration, and //concsolid// is the initial concentration the solid phase. |