qgs is a a Python implementation of a set of idealized reduced-order models representing atmospheric mid-latitude variability. It consists of a two-layer quasi-geostrophic spectral model of the atmosphere on a beta-plane, coupled either to a simple land surface or to a shallow-water ocean. The model’s dynamical fields include the atmospheric and oceanic streamfunction and temperature fields, and the land temperature field