Code

Short codes I have developed as part of my research. Please let me know if you decide to use one or have any questions.

Binning Data by ω500 Velocities

Function to bin data by ω500 velocities, as done by Bony et al. (2004) and many others.

w500_bin.py

Multi-Taper Spectral Analysis

This script includes functions to compute the power spectrum of a signal and to perform a coherence analysis of two time-series using the multi-taper method. These are basically just Python versions of the Matlab functions available from Peter Huybers' website. Requires the spectrum Python module.

multi_taper_analysis.py

Eddy Flux Co-Spectra

Function for calculating eddy flux co-spectra, following method of Hayashi (1971) (see also Randel and Held (1991)). Includes functions to calculate space-time cross-spectra and phase-speed cross-spectrum.

eddy_co_spectra.py

"1-dimensional" 2-layer QG Model

Python script for a 2-layer quasi-geostrophic model in which only the most unstable zonal wavenumber is retained, as described by Panetta and Held (1988). The domain is doubly-periodic and, out of the box, the model is forced by interacting with a zonally-uniform mean-flow, as well as by surface friction and radiative damping, but this can be easily changed. Non-linear terms are de-aliased using the Orszag (1971) 2/3 rule and note that all variables are non-dimensionalized.

Trunc_2l_QG.py

Growth rates and streamfunction