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)

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_Ca : float

activity of calcium (mol/L)

a_H2CO3s : float

activity of carbonic acid (mol/L)

a_H : float

activity of hydrogen (mol/L)

a_HCO3 : float

activity of bicarbonate (mol/L)

T_K : float

temperature degrees Kelvin

PCO2 : float

partial pressure of CO2 (atm)

Returns:

R : float

calcite dissolution rate according to the PWP equation (mmol/cm^2/s)