olm.calcite.concCaEqFromPCO2

olm.calcite.concCaEqFromPCO2(PCO2, T_C=25.0)

Calculates the equilibrium concentration of calcium using PCO2 and temp.

Iteratively solves for the equilibrium concentration of calcium from PCO2 and temp. First guesses that activity coefficients are 1, and then iterates to solution using scipy’s brentq function.

Parameters:

PCO2float, numpy.ndarray, or pandas Series

partial pressure of CO2 (atm)

T_Cfloat, numpy.ndarray, or pandas Series (optional)

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

Returns:

CaEqfloat, numpy.ndarray, or pandas Series

Equilibrium concentration(s) of calcium in mol/L

Notes

Assumes a H20-CO2-CaCO3 system. If a numpy array or pandas Series object are passed in as PCO2 arguments, then equilibrium concentrations will be found iteratively in a for-loop and returned as the same data type given in the argument.