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