olm.calcite.dissRateFromCaPCO2¶
- olm.calcite.dissRateFromCaPCO2(Ca, PCO2, T_C, rho=2.6, method=None, impure=True, per_tol=0.001, error=False, error_num=100, Ca_err=None, PCO2_err=None, molL=False, confidence=0.0, return_samples=False)[source]¶
Calculates the calcite/limestone dissolution rate from given calcium concentration and PCO2. Optionally uses Monte Carlo error propagation to calculate uncertainty in rates.
- Parameters:
- Cafloat, numpy.ndarray or pandas Series
Calcium concentration, default units are mg/L. Change to mol/L by setting keyword mol_L=true.
- PCO2float, numpy.ndarray, or pandas Series
The partial pressure of CO2 for the solution(s).
- T_Cfloat, numpy.ndarray, or pandas Series
The temperature of the water in degrees Celcius.
- rhofloat
Density of rock in g/cm^3. (default=2.6)
- methodstring
Determines method used to calculate dissolution rates. Set to either ‘PWP’ or ‘Palmer’. This keyword is required.
- impureboolean
Used when calculating Palmer rates. Determines whether to use the table values for impure calcite (True) or pure calcite (False). Impure calcite is more representative of typical limestone. (default = True)
- per_tolfloat
the fractional change in H concentration between iterations upon which the iteration is terminated (see solutionFromCaPCO2). default=0.001
- error: boolean
Set to true if you want to use Monte Carlo Error propagation to estimate error in dissolution rate. Requires values for Ca_err and PCO2_err. (default=False)
- error_numinteger
Size of random sample used in Monte Carlo Error propagation. default=100
- Ca_err: float, numpy.ndarray, or pandas Series
Percent error in calcite concentration(s) (1=100%)
- PCO2_err: float, numpy.ndarray, or pandas Series
Percent error in PCO2 values (1=100%)
- molLboolean
Are Ca units in mol/L. If so, set to true. Otherwise, units assumed are mg/L. (default=False, i.e. mg/L)
- confidencefloat
If non-zero then confidence intervals will be used in error estimation (e.g. 90 = 90% confidence). Default is 0.
- return_samplesboolean
If true (default is false), then return entire random samples within error arrays.
- Returns:
- Rfloat, numpy.ndarray, or pandas Series
calcite dissolution rate in mm/yr
- R_errfloat, numpy.ndarray, or pandas Series
error in dissolution rate (returned with R if keyword error=True)