esmtools.carbon.co2_sol

esmtools.carbon.co2_sol(t, s)[source]

Compute CO2 solubility per the equation used in CESM.

Note

See co2calc.F90 for the calculation of CO2 solubility in CESM.

Parameters:
  • t (xarray object) – SST (degC)
  • s (xarray object) – SSS (PSU)
Returns:

Value of solubility in mol/kg/atm

Return type:

ff (xarray object)

References

Weiss & Price (1980, Mar. Chem., 8, 347-359; Eq 13 with table 6 values)

Examples

>>> from esmtools.carbon import co2_sol
>>> import numpy as np
>>> import xarray as xr
>>> t = xr.DataArray(np.random.randint(10, 25, size=(100, 10, 10)),
        dims=['time', 'lat', 'lon'])
>>> s = xr.DataArray(np.random.randint(30, 35, size=(100, 10, 10)),
        dims=['time', 'lat', 'lon'])
>>> ff = co2_sol(t, s)