olm.calcite.pwpRateFranci¶
- olm.calcite.pwpRateFranci(a_Ca=0.0, a_H2CO3s=0.0, a_H=0.0, a_HCO3=0.0, T_K=298.15, PCO2=0.0)[source]¶
Calculates PWP rate using relation for kappa4 used in Franci Gabrovsek’s code.
Calculates PWP rate using relation for kappa4 used in Franci Gabrovsek’s speleogenesis code (pers. commun.). This slight difference was discovered during testing of this code against Franci’s calculations. In this case, a_H in the equation for kappa4 is the bulk value and not the equilibrium surface value for the given carbonic acid concentration.
- Parameters:
- a_Cafloat
activity of calcium (mol/L)
- a_H2CO3sfloat
activity of carbonic acid (mol/L)
- a_Hfloat
activity of hydrogen (mol/L)
- a_HCO3float
activity of bicarbonate (mol/L)
- T_Kfloat
temperature degrees Kelvin
- PCO2float
partial pressure of CO2 (atm)
- Returns:
- Rfloat
calcite dissolution rate according to the PWP equation (mmol/cm^2/s)