Constraining an eddy energy dissipation rate due to relative wind stress for use in energy budget-based eddy parameterisations

Code and data to reproduce results in the EGU Ocean Science journal paper, entitled 'Constraining an eddy energy dissipation rate due to relative wind stress for use in energy budget-based eddy parameterisations'. The data and corresponding scripts in this repository are: MITgcm simulation data for an anticyclone and cyclone under absolute and relative wind stress. ACE_eta_{absolute or relative}_10km_A25_v2_ext.nc ACE_tau_{absolute or relative}_10km_A25_v2_ext.nc ACE_temp_{absolute or relative}_10km_A25_v2_ext.nc CE_eta_{absolute or relative}_10km_A25_v2.nc CE_tau_{absolute or relative}_10km_A25_v2.nc CE_temp_{absolute or relative}_10km_A25_v2.nc Mean energetics computed using eddy_energy_mean.m. User may need to comment out lines of code relating to wvel. This also requires geostrophic_uv.m, rho_ref.mat, and ocean_vertical_grid.nc. {ACE or CE}_energy_{absolute or relative}_total_10km_A25_mean.nc Predicted eddy energy computed using predict_*.m scripts. n.b. viscous term left in absolute scripts, though we only need the first index of total eddy energy in abs cases. The abs cases are really just making the abs.mat files, see time_series_plots.py. Time-series figures (3 and 5) produced using time_series_plots.py. Initial figures (1 and 2) produced using initial_plots.py. Relative vorticity in Fig. 4 computed and plotted using eddy_evolution_plot.py. Dissipation rates computed in calc_dissipation_rate.m. Requires dynmodes.m and GSW toolbox to run. Also needs data Chelton2011_Le.csv and Chelton2011_Rd.csv. diss_rate*.nc Dissipation rate and climatology figures (6, 7, 8, and 9) made in diss_rate_plot.py Information and data needed to run the MITgcm can be found here.