Analysis code to: Pinner et al., 2024

This is the analysis code to the for open peer-review submitted publication
Pinner et al., 2024: Internal-wave-induced dissipation rates in the Weddell Sea Bottom Water gravity current,

The Weddell Sea Bottom Water gravity current transports dense water from the continental shelf to the deep sea and is crucial for the formation of new deep sea water. Build on vertical profiles and time series measured in the northwestern Weddell Sea, we apply 3 methods to distinguish turbulence caused by internal waves from turbulence by other sources. We find that in the upper part of the gravity current, internal waves are important for the mixing of less dense water down into the current.

Disclaimer

Important

• Although this code produces the results and figures to the accompanying paper, this repository still contains unused code snippets and unfinished documentation. The repository will be cleaned up further during the review process, while making sure that the output remains the same. This may include the renaming of files or functions to better communicate their purpose.
• Comments or corrections to the code can be given on GitHub as issues or discussions or on EGUsphere as a community reviewer.
• Note that figures created here can differ slightly from the submitted versions, as some post-processing (adjustements and labeling) were made with Inkscape.

Structure

The figures 1 to 6 are produced from data as follows:

flowchart TD
    bathymetry[(bathymetry)] ----> f01((f01: overview))
    CTD1[(CTD profiles)] --> id[[""eos80_legacy_gamma_n matlab toolbox""]]
    id --> gamma["neutral density γⁿ"]
    gamma --> f01

    f02((f02: spectrum))
    TS[(velocity time series)] ----> f02
    TS ----> f03((f03: flowfield))

    TS --> idemix[[wave energy/IDEMIX parameterization]]
    idemix --> epsidemix["ε_{IGW, IDEMIX}"]

    CTD2[(CTD profiles)]
    CTD2 --> fine[[finestructure]]
    fine --> epsfine["ε_{IGW, fine}"]
    f05(("f05: ε_{IGW} transect"))
    epsfine --> f05
    epsidemix --> f05
    epsidemix --> f04(("f04: ε transect"))
    CTD2 --> Thorpe[[Thorpe scales]]
    Thorpe --> epstotal["ε_{total}"]
    epstotal --> f04

    f06(("f06: ε profile"))
    epstotal --> f06
    epsidemix --> f06
    epsfine  --> f06

Loading

• Data (1st level) is referenced from the ./data directory (but not committed to repository)
• Methods (2nd level) are in the ./scripts directory
• Methods results (3rd level) are in the respective ./scripts/<method_name>/method_results directories
• Figures (4th level) are produced by the scripts in the ./results directory