olm.calcite.activityHFromPCO2

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

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:
PCO2float

partial pressure of CO2 (atm)

T_Cfloat, optional

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

CaEqfloat

Equilibrium calcium concentration (mol/L), optional

Returns:
aHeqfloat

equilibrium activity of hydrogen ion (mol/L)

Notes

Assumes a H20-CO2-CaCO3 system.