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.

  • flip (bool, optional) – Whether to flip the input data along the latitude dimension. Use this if your zonal average gradient is negative poleward.

  • omega (ndarray, optional) – The vertical wind. Shape must match q.

  • sigma (ndarray, 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.