olm.calcite.pwpRateTheory

olm.calcite.pwpRateTheory(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 theoretical relation for kappa4 from PWP.

Calculates PWP rate using theoretical relation for kappa4 from PWP (as described by equation 25 in Plummer, Wigley, and Parkhurst (1978) and in Dreybrodt [1988] equation 6.22b). In this case, a_H in the equation for kappa4 is the equilibrium surface value for the given carbonic acid concentration, as specified in the theory.

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)