This README.txt file was generated on 16 OCT 2024 by Sam Kelly ------------------- GENERAL INFORMATION ------------------- These Matlab scripts and data files replicate the figures in Kelly et al. (2024). See the paper captions for detailed figure descriptions. All input data files can be generated using the Matlab scripts in the "analysis_scripts" folder. ------------------------------------ SCRIPT DESCRIPTION OF plot_fig1.m ------------------------------------ This script plots figure 1. Function dependencies: (1) LAT2KM.m (2) LON2KM.m (3) FIXXLIM.m Input files: (1) data.mat (2) trans.mat (3) bathy.mat Output files: (1) fig1.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig2.m ------------------------------------ This script plots figure 2. Function dependencies: (1) soradna1.m (air-sea toolbox) (2) hfbulktc.m (air-sea toolbox) (3) blwhf.m (air-sea toolbox) Input files: (1) data.mat (2) wind.mat (3) sadcp.mat Output files: (1) fig2.pdf ------------------------------------ SCRIPT DESCRIPTION OF plot_fig3.m ------------------------------------ This script plots figure 3. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat (2) slab.mat (3) model.mat Output files: (1) fig3.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig4.m ------------------------------------ This script plots figure 4. Function dependencies: (1) cbrewer2.m Input files: (1) trans.mat Output files: (1) fig4.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig5.m ------------------------------------ This script plots figure 5. Function dependencies: none Input files: (1) trans.mat Output files: (1) fig5.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig6.m ------------------------------------ This script plots figure 6. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat Output files: (1) fig6.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig7.m ------------------------------------ This script plots figure 7. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat Output files: (1) fig7.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig8.m ------------------------------------ This script plots figure 8. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat Output files: (1) fig8.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig9.m ------------------------------------ This script plots figure 9. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat Output files: (1) fig9.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig10.m ------------------------------------ This script plots figure 10. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat Output files: (1) fig10.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig11.m ------------------------------------ This script plots figure 11. Function dependencies: none Input files: (1) data.mat (2) trans.mat (3) slab.mat Output files: (1) fig11.pdf ------------------------------------ SCRIPT DESCRIPTION OF plot_fig12.m ------------------------------------ This script plots figure 12. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat (2) model.mat Output files: (1) fig12.png ------------------------------------ SCRIPT DESCRIPTION OF plot_fig13.m ------------------------------------ This script plots figure 13. Function dependencies: (1) REDBLUE5.m Input files: (1) trans.mat (2) slab.mat (3) model.mat Output files: (1) fig12.png ------------------------------------ FUNCTION DESCRIPTION OF cbrewer2.m ------------------------------------ Colormap developed by Cynthia Brewer (http://colorbrewer2.org/). Try https://www.mathworks.com/matlabcentral/fileexchange/58350-cbrewer2 ------------------------------------ FUNCTION DESCRIPTION OF LAT2KM.m ------------------------------------ Input: None Output: The number of kilometers per degree of latitude. ------------------------------------ FUNCTION DESCRIPTION OF LON2KM.m ------------------------------------ Input: Latitude [deg] Output: The number of kilometers per degree of longitude. ------------------------------------ FUNCTION DESCRIPTION OF FIXXLIM.m ------------------------------------ Input: gca (Matlab handle, see Matlab documentation) right x-limit scale (i.e., the desired x/y aspect ratio of the figure) Output: The x-limits that have the desired aspect ratio (given the y-limits and the right x limit). ------------------------------------ FUNCTION DESCRIPTION OF REDBLUE5.m ------------------------------------ A custom Matlab color map that can is called like a standard color map (see Matlab documentation). ------------------------------------ TOOLBOX DESCRIPTION of air-sea ------------------------------------ Surface fluxes are computed using functions from the air-sea toolbox at https://github.com/sea-mat/air-sea BLWHF: estimates net long-wave heat flux using bulk formulas HFBULKTC: computes sensible and latent heat fluxes and other variables. SORADNA1: computes no-sky solar radiation and solar altitude. ------------------------------------ DATA FILE DESCRIPTION OF bathy.mat ------------------------------------ Great Lakes Bathymetry (https://www.ncei.noaa.gov/products/great-lakes-bathymetry) in Matlab format Structure: bathy fields: x (1 x 9841) [m] distance east y (1 x 4201) [m] distance north lon (1 x 9841) [deg] longitude lat (1 x 4201) [deg] latitude H (9841 x 4201) [m] bottom depth ------------------------------------ DATA FILE DESCRIPTION OF data.mat ------------------------------------ Observations from shipboard transects (data) and mooorings (moorings). Structure: data fields: lon (2794 x 1) [deg] longitude lat (2794 x 1) [deg] latitude H (2794 x 1) [m] bottom depth dist (2794 x 1) [m] distance from shore along transect z (2000 x 1) [m] depth time (2794 x 1) [days] days since 1-Jan-0000 (Matlab datenum) T (2000 x 2794) [C] water temperature rho (2000 x 2794) [C] water density u (2000 x 2794) [m/s] east water velocity v (2000 x 2794) [m/s] north water velocity eps (2000 x 2794) [W/kg] TKE dissipation rate tau (2794 x 1) [Pa] wind stress (real = east, imag = north) p_atm (2794 x 1) [Pa] atmospheric pressure N2 (2000 x 2794) [1/s^2] buoyancy frequency squared S2 (2000 x 2794)) [1/s^2] vertical shear squared Structure: moorings fields: name (1 x 3) mooring names lon (1 x 3) [deg] longitude lat (1 x 3) [deg] latitude H (1 x 3) [m] bottom depth z (1 x 200) [m] depth time (1 x 43825) [days] days since 1-Jan-0000 (Matlab datenum) T (43825 x 200 x 3) [C] water temperature u (43825 x 200 x 3) [m/s] water velocity (real = east, imag = north) Note: several unused fields in data.mat that are not documented here. ------------------------------------ DATA FILE DESCRIPTION OF sadcp.mat ------------------------------------ Shipboard currents (ADCP) and meteorological data Structure: met fields: name (1 x 13) weather station names code (1 x 13) weather station number lon (17830 x 1) [deg] longitude lat (17830 x 1) [deg] latitude time (17830 x 1) [days] days since 1-Jan-0000 (Matlab datenum) H (17830 x 1) [m] bottom depth U10 (17830 x 1) [m/s] 10-m wind velocity (real = east, imag = north) tau (17830 x 1) [Pa] wind stress (real = east, imag = north) p_atm (17830 x 1) [Pa] atmospheric pressure Ta (17830 x 1) [C] air temperature Tw (17830 x 1) [C] water temperature RG (17830 x 1) [%] relative humidity Note: several unused structures and fields in sadcp.mat are not documented here. ------------------------------------ DATA FILE DESCRIPTION OF slab.mat ------------------------------------ Slab model predictions Structure: slab fields: x (1 x 41) [m] distance from shore z (1 x 2000) [m] depth hr (97 x 1) [hr] time in hours f (1 x 1) [1/s] inertial frequency tau (97 x 1) [Pa] wind stress (real = east, imag = north) uI (1 x 1) [m/s] initial velocity eps_P (2000 x 97) [W/kg] TKE shear production eps_r (2000 x 97) [W/kg] dissipation by linear drag P (97 x 1) [W/m^2] TKE shear production D (97 x 1) [W/m^2] dissipation by linear drag WW (97 x 1) [W/m^2] wind work dEdt (97 x 1) [W/m^2] energy tendency u0 (2000 x 41 x 97) [m/s] observed velocity u (2000 x 41 x 97) [m/s] slab-model velocity F (2000 x 1) [] vertical profile of stress divergence (forcing) Z (2000 x 1) [] vertical profile of stress Note: several unused fields in slab.mat are not documented here. ------------------------------------ DATA FILE DESCRIPTION OF trans.mat ------------------------------------ Shipboard observations organized into transects Structure: trans fields: x (41 x 1) [m] distance from shore z (2000 x 1) [m] depth lon (41 x 1) [deg] longitude lat (41 x 1) [deg] latitude H (41 x 1) [m] bottom depth time (97 x 1) [days] days since 1-Jan-0000 (Matlab datenum) hr (97 x 1) [hr] hours since start of survey Tref (2000 x 41) [C] background temperature rho_ref (2000 x 41) [kg/m^3] background density N2_ref (2000 x 41) [1/s^2] background buoyancy frequency squared u (2000 x 41 x 97) [m/s] east velocity v (2000 x 41 x 97) [m/s] north velocity T (2000 x 41 x 97) [C] water temperature p (2000 x 41 x 97) [p] baroclinic pressure b (2000 x 41 x 97) [m/s^2] buoyancy rho (2000 x 41 x 97) [kg/m^3] density N2 (2000 x 41 x 97) [1/s^2] buoyancy frequency squared S2 (2000 x 41 x 97) [1/s^2] velocity shear squared eps (2000 x 41 x 97) [W/kg] TKE dissipation rate eps0 (41 x 97) [W/m^2] depth-integrated TKE dissipation rate tau (41 x 97) [Pa] wind stress (real = east, imag = north) c (41 x 10) [m/s] eigenspeed of first 10 vertical modes PHI (2000 x 41 x 10) [] vertical modes un (41 x 97 x 10) [m/s] east velocity modal amplitude vn (41 x 97 x 10) [m/s] north velocity modal amplitude pn (41 x 97 x 10) [Pa] pressure modal amplitude Note: several unused fields in trans.mat are not documented here. ------------------------------------ DATA FILE DESCRIPTION OF wind.mat ------------------------------------ Wind data from several NOAA weather stations Structure: wind fields: name (1 x 13) weather station names code (1 x 13) weather station number lon (1 x 13) [deg] longitude lat (1 x 13) [deg] latitude time (26281 x 1) [days] days since 1-Jan-0000 (Matlab datenum) U10 (26281 x 13) [m/s] 10-m wind velocity (real = east, imag = north) tau (26281 x 13) [Pa] wind stress (real = east, imag = north) p_atm (26281 x 13) [Pa] atmospheric pressure Ta (26281 x 13) [C] air temperature Tw (26281 x 13) [C] water temperature Note: several unused fields in trans.mat are not documented here.