Balanced simulations of Eady turbulence and strain-induced frontogenesis

This is the repository for code used in the paper

• Dù, R.S., Smith, K.S. , Bühler, O., 2024. Next-order balanced model captures submesoscale physics and statistics

It contains simulation and analysis of simulations of Eady turbulence under the QG⁺¹ model and of strain-induced frontogenesis under the QG⁺¹ and semigeostrophic model.

Instructions for Eady turbulence under QG⁺¹

• The file Eady_QGpl.py simulates Eady turbulence under the QG⁺¹ model according to the parameters in the paper.
• Use the notebook plot_snap_both.ipynb in the plot folder to generate Fig. 1-6 in the paper.
• Use the notebook plot_jointpdf.ipynb to generate the latter two figures in Fig. 7-8 in the paper.
• Use the notebook plot_btend.ipynb to generate Fig. 9 in the paper.

Instructions for strain-induced frontogenesis simulations

• Use the notebook IC_gen.ipynb to generate the initial conditions for the simulations. It generates the first figure in Fig. 11 in the paper.
• The file 2DFrontQGpl.py and 2DFrontSemiG.py simulate strain induced frontogenesis under the QG⁺¹ and semigeostrophic model.
• Use the notebook Front_plotboth.ipynb in the plot folder to generate the latter two figures in Fig. 12 in the paper.
• Use the notebook maxby_growth.ipynb to generate Fig. 13 in the paper.