Scripts for the numerical experiments in JPO-D-23-0059 "Vertical Exchange Induced by Mixed-Layer Instabilities"

Description

These four files are the scripts for the numerical experiments in JPO-D-23-0059 "Vertical Exchange Induced by Mixed-Layer Instabilities".

growth_rate.py is a Python script that runs Dedalus to solve an eigenvalue problem in the Section "Linear Stability Analysis" in the paper. It computes the growth rate of the mixed layer instability at different x-direction wavenumbers, given the Richardson numbers of the mixed layer and the thermocline and the depth ratio of these two layers.

vertical_structure.py is also a Python script that runs Dedalus to solve an eigenvalue problem in the Section "Linear Stability Analysis" in the paper. It computes the vertical structures of the eigenfunctions for u, v, w, p, and b at a given x-direction wavenumber, given the Richardson numbers of the mixed layer and the thermocline and the depth ratio of these two layers.

nonlinear_simulations.py is a Python script that runs Dedalus to solve an initial value problem in the Section "Nonlinear Simulations" in the paper. It computes the fields of u, v, w, p, and b at all times, given the Richardson numbers of the mixed layer and the thermocline.

dedalus.cfg is a Dedalus configuration file that should be placed in the same directory as nonlinear_simulations.py when running the latter.

Files

Additional details