Shallow-water modelling in Python

This repository contains two implementations of the shallow-water equations that are suitable to study a wide range of wave and ocean circulation phenomena, including non-linear effects.

They are a product of the Bornö summer school 2018, led by Markus Jochum and Carsten Eden.

A preview of the non-linear setup:

Features

Simple (linear) implementation

• Mass conserving on (Cartesian) Arakawa C-grid
• Mixed-time discretization
• Coriolis force on an f-plane
• Conditionally stable for Δt <= √2 / f

Fully non-linear implementation

All features of the simple implementation, plus...

• Adams-Bashforth time stepping scheme
• Lateral friction
• Varying Coriolis parameter (β-plane)
• Fully non-linear momentum and continuity equations
• Energy conserving scheme by Sadourny (1975)
• Rigid or periodic boundary conditions

Problems?

Feel free to open an issue or submit a PR!