oceanspy.utils.densmdjwf

oceanspy.utils.densmdjwf(s, t, p)[source]

Density of Sea Water using McDougall et al. 2003 (JAOT 20) polynomial (Gibbs Potential) [McJa03]. mdjwf.py: https://github.com/MITgcm/MITgcm/blob/master/utils/python/MITgcmutils/MITgcmutils/mdjwf.py

Parameters:

s: xarray.DatArray, array-like: salinity [psu (PSS-78)]
t: xarray.DatArray, array-like: potential temperature [degree C (IPTS-68)]
p: xarray.DatArray, array-like: pressure [dbar] (p may have dims 1x1, mx1, 1xn or mxn for S(mxn))

Returns:

rho: xarray.DatArray, array-like: density [kg/m^3]

References

[McJa03]

McDougall, T.J., D.R. Jackett, D.G. Wright, and R. Feistel, 2003: Accurate and Computationally Efficient Algorithms for Potential Temperature and Density of Seawater. J. Atmos. Oceanic Technol., 20, 730–741, https://doi.org/10.1175/1520-0426(2003)20<730:AACEAF>2.0.CO;2