Heat transport across the Antarctic Slope Front controlled by cross-slope salinity gradients

Feb 2023 updates: Add code for EKE spectral analysis to MITgcm_ASF-heat-ver3/analysis/spectrum/ Add products of 5km and 10km runs to products_new-ver3 Add MITgcm source code, copied from http://mitgcm.org This release contains updates on analysis code and products. MITgcm_ASF-heat-ver3/<strong>newexp</strong>/: the Matlab scripts used to generate and run the MITgcm simulations MITgcm_ASF-heat-ver3/<strong>analysis</strong>/<strong>cross_slope</strong>/ and MITgcm_ASF-heat-ver2/<strong>analysis</strong>/<strong>plots</strong>/: the Matlab scripts used to analyze model output and make plots. MITgcm_ASF-heat-ver3/analysis/<strong>spectrum</strong>/: the<strong> </strong>Matlab<strong> </strong>scripts to calculate EKE spectra<strong> </strong> exps_configuration.zip: the configurations of the MITgcm simulations. <strong>products_new-ver3.zip</strong>: the products calculated from MITgcm diagnostics, including 7-year means of all the model outputs, overturning streamfunctions, neutral density, shoreward heat transport, kinetic energy, temporal decomposition, isopycnal thickness fluxes of the 5km and 10km runs, etc. ThicknessFlux_FreshShelf.zip: products of isopycnal thickness flux, used to calculate the decomposition of eddy/tidal heat advection/diffusion, for the "fresh-shelf" simulation. ThicknessFlux_ref.zip: as above, but for the reference simulation. ThicknessFlux_DenseShelf.zip: as above, but for the "dense-shelf" simulation. The source code of the Massachusetts Institute of Technology General Circulation Model (MITgcm) is available at: http://mitgcm.org. <strong>All the raw data of the model output are available at: https://doi.org/10.15144/S47P49.</strong> To reproduce MITgcm_ASF simulations: Start each simulation with a 20-year spin-up integration. Before running each simulation, you need to substitute <em>&amp;OBCS_PARM04</em> with <em>&amp;OBCS_PARM05 </em>in the file <em>input/<strong>d</strong>ata.obcs</em>, and substitute <em>&amp;EXF_NML_05 </em>with <em>&amp;EXF_NML_OBCS</em> <em> </em>in the file <em>input/data.exf</em>. For simulations with very fresh shelf waters (e.g., shelf salinity = 33 psu), you need to spin up the simulation with a very small time step (e.g., 60s) for ~ two months, and then use a larger time step. Initialize the production run from the corresponding spin-up run, using the Matlab script <em>initialize.m</em> in the folder<em> MITgcm_ASF-heat-ver2/newexp/. </em>When using the LAYERS package, you need to substitute<em> numperlist = 1</em> with <em>numperlist = 2 </em>in the file<em> code/DIAGNOSTICS_SIZE.h</em> before running the simulations. Notes on calculationg the overturning streamfunction and its mean/eddy/tidal decomposition using the MITgcm LAYERS package: avg_t: Calculate time averages. It has been modified since the vertical number of layers can be different from Nr. calc_Overturning_pt, usscar_plot_overturning_pt: calculate and plot eddy/mean/isopycnal overturning streamfunction using potential temperature layer fluxes. calc_Overturning_rho, usscar_plot_overturning_rho: calculate and plot eddy/mean/isopycnal overturning streamfunction using potential density layer fluxes. calc_Overturning_pt_Aocean, usscar_pt_overturning_rho_Aocean (<strong>recommended if your bathymetry is not flat</strong>): calculate and plot eddy/mean/isopycnal overturning streamfunction using <em>potential temperature</em> layer fluxes. For each latitude, use the total ocean area below a certain level to interpolate the streamfunction from pt space to z space. calc_Overturning_rho_Aocean, usscar_plot_overturning_rho_Aocean (<strong>recommended <strong>if your bathymetry is not flat</strong></strong>): calculate and plot eddy/mean/isopycnal overturning streamfunction using <em>potential density</em> layer fluxes. For each latitude, use the total ocean area below a certain level to interpolate the streamfunction from potential density space to z space. calc_decomposition_OT, plot_OT_rho_Aocean_TidalEddyMean: decompose the isopycnal overturning streamfunction into <strong>tidal</strong>/eddy/mean components, using potential density layer fluxes. Feel free to contact Yidongfang Si via <strong>[e-mail]</strong> if you have any questions.