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

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