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_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)