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 : <strong>R. Torres, R. Waldman, J. Mak and R. Séférian </strong>: <em>Global Estimation of the Eddy Kinetic Energy Dissipation from a Diagnostic Energy Balance</em>. 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 (<em>CMEMS/SEALEVEL_GLO_PHY_MDT_008_063/P20Y</em>) downloaded here and ocean masks (<em>misc/basins/doi_10.5281</em>) 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 <em>CNRM/runs/omip2_LR.Geom_Emin0-alpha01_1cyc-trd/post</em>. 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 <em>EKE_dissipation_rate</em> folder since it results from a combination of multiple products. IPython notebooks for computing and plotting global maps are also provided : <em>1-post_process_climato.ipynb</em> : 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). <em>2-post_process_altimetry.ipynb</em> : 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. <em>3-compute_global_eke_dissipation.ipynb</em> : combine both outputs from the two above scripts to compute the global EKE dissipation. The script also plots new maps. <em>0-plot_global_maps.ipynb</em> : plot the global maps for climatology and altimetry products. <em>0-plot_lbekedis_ogcm.ipynb</em> : 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.