Published November 8, 2022 | Version v1
Dataset Open

Model setup and output for 'Tidal conversion and dissipation at steep topography in a channel poleward of the critical latitude'

  • 1. University of Victoria


Data supplement to Hughes and Klymak 2019

Model input and output in a reduced form associated with the following paper:

Tidal conversion and dissipation at steep topography in a channel poleward of the critical latitude

Journal of Physical Oceanography. doi:10.1175/JPO-D-18-0132.1


As described in Table 1 of the associated paper, there are three main sets of simulations. The input files for these sets are contained in their respective directories (`vary_width`, `vary_forcing`, and `vary_freq`). The python script that creates all of the necessary files is ``. A fourth directory is titled `baroclinic_terms` and includes that simulation in which u' and p' are output at high temporal resolution.

A key point regarding the input files is that for the vary width and vary forcing cases, a single simulation involves multiple channels. This lets me compile a single executable `mitgcmuv` with a Nx × Ny grid of 600 × 1280, which I divide up into the necessary number of channels by putting vertical walls in appropriate places. For the vary width cases, the 'narrow' simulations are all channels from 0.2 to 32 km and the 'wide' simulations are all wider channels. Once the simulation has run, I use netcdf tools (`ncks`) to extract the individual channels using the scripts in the `extract_scripts` directory.

Most of the files in the `code` directories will be familiar to anyone that uses the MITgcm. An exception is the `energy_diagnostics_fill.F` (and `diagnostics_main_init.F` and `do_statevars_diags.F`, which have minor additions). The original, from ``, was modified slightly to suit this project.


The results directory contains five subdirectories to be described in turn.

Notes that in all cases, energy terms in the netCDF files do not include a factor of ρ. This was added in at the plotting stage.

All simulations used Checkpoint67b and were run on Graham:


The majority of the files are of the form `` where `FFF` is $1000 ω/f$ and `YY` is the channel width in kilometres. These files contain the tidally averaged, depth-integrated energy diagnostics for the seventh tidal cycle at all points (x, y) within the energy control volume.

There are also three files entitled ``, which contain fields of U, V, and T (which gives density with α = 0.0002) at two levels. These fields are used as examples for weakly and strongly responding channels.


These files are of the form `` where `UU` is the deep-water tidal current amplitude U_0 in cm/s. They contain the same energy terms as for the vary width simulations.


These files are of the form `` and contain fields of U, V, and T at two levels. Energy terms are not included because the vary frequency simulations were only run to get estimates of the along-ridge wavelength.


The single file within this directory contains u' and p' at a single x position every five minutes for four tidal cycles.


This directory, named for its obstacle and width, contains U, V, and T at every grid point for a snapshot in time and another file with the corresponding snapshots of all energy terms.


Files (440.5 MB)

Name Size Download all
11.5 kB Preview Download
2.1 kB Preview Download
7.9 kB Download
172.7 kB Preview Download
439.8 MB Preview Download
33.5 kB Preview Download
97.2 kB Preview Download
333.3 kB Preview Download

Additional details

Related works

Is supplement to
Journal article: 10.1175/JPO-D-18-0132.1 (DOI)