Script in Matlab, R and Python that computes the Work done by buoyancy required to bring a parcel of water from different depths (zint) to a reference depth (defualt = -10 m). It is considered a resting ocean and that the plot in zint is in hydrostatic equilibrium. The vertical discretization of the vertical profile can be non-equidistant.
To view the demos, go to the demo folder, choose your preferred programming language, and run it.
Important
If you use this method, please include a reference to the following:
Moreles, E., Romero, E., Ramos-Musalem, K., and Tenorio-Fernandez, L.: The global ocean mixed layer depth derived from an energy approach based on buoyancy work, Ocean Sci., 21, 2019–2039, https://doi.org/10.5194/os-21-2019-2025, 2025.
Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it. (http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System.
This work is licensed under the GNU General Public License v3.0.
