Scripts and datas for "Global Estimation of the Eddy Kinetic Energy Dissipation from a Diagnostic Energy Balance"

Input and output datasets used for a global reconstruction of the eddy kinetic energy (EKE) dissipation rate in relation to a submitted work :

R. Torres, R. Waldman, J. Mak and R. Séférian : Global Estimation of the Eddy Kinetic Energy Dissipation from a Diagnostic Energy Balance.

Inputs datas include a merge of 2 datasets from the World Ocean Atlas 2018 (WOA18, Garcia et al., 2019) and cover the 1995-2017 (95B7) period. Folder structure for the surface altimetry L4 datasets from the EU-Copernicus Marine Services (2021) is kept empty in order to limit the archive size. Datas can be download here.

Optional datasets include CMEMS MDT product (CMEMS/SEALEVEL_GLO_PHY_MDT_008_063/P20Y) downloaded here and ocean masks (misc/basins/doi_10.5281) from Martinez-Moreno et al. (2021).

In addition, simulation outputs from the NEMO-OMIP2 model runned with the GEOMETRIC parameterization are processed (mainly time-averaged) and stored in CNRM/runs/omip2_LR.Geom_Emin0-alpha01_1cyc-trd/post. These files are used in the uncertainties and errors quantification.

Outputs and published results are stored in each individual product post-processing folder while final EKE dissipation computation are located in the EKE_dissipation_rate folder since it results from a combination of multiple products.

IPython notebooks for computing and plotting global maps are also provided :

1-post_process_climato.ipynb : compute from the climatology (e.g. WOA18 datas) the EKE dissipation timescales (units in days) and the surface modes with rough topography (LaCasce and Groeskamp, 2020).

2-post_process_altimetry.ipynb : compute from altimetry (CMEMS) datasets the EKE at surface and eventually coarsen the grid from 0.25 to 1 degree in order to match the climatology grid.

3-compute_global_eke_dissipation.ipynb : combine both outputs from the two above scripts to compute the global EKE dissipation. The script also plots new maps.

0-plot_global_maps.ipynb : plot the global maps for climatology and altimetry products.

0-plot_lbekedis_ogcm.ipynb : plot and analyse EKE timescale errors from the NEMO-OMIP2 simulation outputs.

Note however that these scripts use the author python library XOCE availbale on GitHub: https://github.com/torresr-cnrm/xoce. All the scripts have been runned using the version 0.2 of XOCE. Feel free to contact ([e-mail]) for any help in installing and using this library.