This repository hosts Julia code for obtaining numerical solutions to the planetary geostrophic equations (with momentum fluxes parameterized by Fickian diffusion) for general topography. For a similar model using Rayleigh drag, see https://github.com/joernc/pgcm. The 1D and 2D models are described in Peterson & Callies (2022): Rapid spin up and spin down of flow along slopes, https://doi.org/10.1175/JPO-D-21-0173.1.
The PG models are in src/
(3D under development).
See examples/
for some examples of running them.
non_pg_models/
also has some code for a non-PG 1D model (both dimensional and non-dimensional).
Please don't hesitate to contact me if you have any questions!