waqlocal

waqlocal(lon, lat, q, omega=None, sigma=None, flip=True, skip=10)[source]

Return the local finite-amplitude wave activity. See [HN16] for details.

Parameters

  • lon, lat (ndarray) – The longitude and latitude coordinates.

  • q (ndarray) – The potential vorcitiy data. Should be sorted lon x lat x … (this will be changed in a future version).

  • skip (int, optional) – The interval of sorted q from which we select output Q values.

flipbool, optional

Whether to flip the input data along the latitude dimension. Use this if your zonal average gradient is negative poleward.

omegandarray, optional

The vertical wind. Shape must match q.

sigmandarray, optional

The data weights, useful for isentropic coordinates. Shape must match q.

References

HN16

Clare S. Y. Huang and Noboru Nakamura. Local Finite-Amplitude Wave Activity as a Diagnostic of Anomalous Weather Events. Journal of the Atmospheric Sciences, 73(1):211–229, 2016. doi:10.1175/JAS-D-15-0194.1.

Todo

Support using omega.