olm.calcite.pwpRatePascal

olm.calcite.pwpRatePascal(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 found in PASCAL code.

Calculates PWP dissolution rate for calcite using the relation for kappa4 that is found in the PASCAL code in Dreybrodt (1988). This is also given in equation 3 from Buhmann and Dreybrodt (1985), The kinetics of calcite dissolution and precipitation in geologically relevant situations of karst areas: 1. Open system. They say that it is a fit to the experimental data of Plummer et al. for values of PCO2 < 0.05 atm.

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)