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