Code to set up the simulations in Arnscheidt et al. (2021) "On the settling depth of meltwater escaping from underneath an Antarctic ice shelf". Each model runs from a single script, and the relevant parameter values can be straightforwardly modified within that script.

The simple line plume model can be run from 1d_line_plume.py with a working Python3 environment.

The large-eddy simulations run in Julia using Oceananigans.jl. Ideally, the more computationally intensive simulations should be run using Graphical Processing Units (GPUs). 2-D meltwater plume large-eddy simulations can be run from les_2d_constant.jl and les_2d_restoring.jl, which use a constant buoyancy source and a restoring buoyancy source respectively. 3-D large-eddy simulations can be run from les_3d.jl. These rely on a specific version of Oceananigans.jl and thus require the Project.toml file: they can be run using julia --project and then typing include("filename.jl").

pig_data contains temperature and salinity profiles from near the meltwater outflow from beneath Pine Island glacier in 2009 and 2014; these can be used as initial conditions in les_3d.jl.