Python Ocean Model 2.0 (pyOM2)
Introduction
pyOM2.0 (Python Ocean Model) is a numerical circulation ocean model which was written for educational purpose. It is meant to be a simple and easy to use numerical tool to configure and to integrate idealized and realistic numerical simulations of the ocean in Boussinesq approximation. Non-hydrostatic situations as well as large-scale oceanic flows can be considered, Cartesian or pseudo-spherical coordinate systems can be used. Several idealized experiments and examples are preconfigured and can be easily chosen and modified using two alternative configuration methods based on Fortran90 or Python. Prerequisites for the installation is a Fortran 90 compiler and the Lapack library, and for the Fortran front the NetCDF-library (since IO is realized mainly using the NetCDF format). For the Python front end, the numerical module \verb+numpy+ is required and several other modules can be used in addition, e.g. to provide
- a graphical user interface. Both version are based on identical Fortran90 code which is fully parallelized based on the MPI-library to enhance performance.
Resources
Prerequisites and Installation
Prerequisites
...
Installation
...
Sample Configurations
Kelvin-Helmholtz Instability
kelvin_helm1.py
Rayleigh–Bénard Convection
rayleigh.py
Eady's Baroclinic Instability
eady1.py / eady2.py
Eddy-driven zonal jets
jets1.py
Thermohaline Circulation
THC1.py
Southern Ocean Circulation
acc1.py
ENSO Response
enso1.py
Equatorial Waves
eq_waves1.py
Isopycnal Diffusion
isopyc_test1.py
An example of the python GUI for Eady's baroclinic instability case is shown below.
