Barotropic-Python

A simple barotropic model written in Python using the spharm package for spherical harmonics.

Currently set up to use NCEP reanalysis for initial conditions

Written by Luke Madaus (5/2012) - University of Washington

Requires: PySPHARM -- Python interface to NCAR SPHEREPACK library: https://code.google.com/p/pyspharm/

netCDF4 -- Python interface to netCDF4 library matplotlib numpy

Based on the Held-Suarez Barotropic model, including hyperdiffusion. A brief description of their model may be found at: http://data1.gfdl.noaa.gov/~arl/pubrel/m/atm_dycores/src/atmos_spectral_barotropic/barotropic.pdf

The basic premise of the code is to take 500 hPa u and v winds from the NCEP reanalysis, extract the non-divergent component, compute vorticity, and advect along this vorticity using the barotropic vorticity equation. As this model uses "real" atmospheric data (which is not barotropic), the results are rarely stable beyond ~5 days of forecasting.