francispoulin / QGModel

Python code that solves nonlinear QG equations

QGModel

This is Python code that solves the nonlinear QG equations for PV anomalies on a doubly-periodic domain. It consists of two parts:

1. a fully-dealiased pseudo-spectral code that steps forward the PV conservation equations and
2. implementations of specific model types that provide an inversion relation.

At the moment, seven model types have been implemented:

1. two-dimensional dynamics (TwoDim),
2. surface QG dynamics (Surface),
3. multi-layer dynamics (Layered),
4. Eady dynamics (Eady),
5. floating Eady dynamics (FloatingEady),
6. two-Eady dynamics (TwoEady),
7. two-Eady dynamics with buoyancy jump (TwoEadyJump).

See run.py for an example of how a model is initialized and run.

This model makes use of PyFFTW, a Python wrapper of FFTW. Follow the installation instructions for PyFFTW if you do not have it installed already.