olm.calcite.activityHFromPCO2

olm.calcite.activityHFromPCO2(PCO2, T_C=25.0, CaEq=None)

Calculates equilibrium activity of H+ given PCO2 using relaxed charge balance.

Calculates hydrogen activity at equilibrium given PCO2, temperature, and (optionally) equilibrium calcium concentration (mol/L). Assumes a relaxed charge balance (see 2.18a in Dreybrodt [1988]). If keyword CaEq is not given, then it is iteratively calculated using concCaEqFromPCO2().

Parameters:

PCO2 : float

partial pressure of CO2 (atm)

T_C : float, optional

temperature of solution in degrees Celsius (default = 25 C)

CaEq : float

Equilibrium calcium concentration (mol/L), optional

Returns:

aHeq : float

equilibrium activity of hydrogen ion (mol/L)

Notes

Assumes a H20-CO2-CaCO3 system.