brunt_vaisala_frequency#
- metpy.calc.brunt_vaisala_frequency(height, potential_temperature, vertical_dim=0)#
Calculate the Brunt-Vaisala frequency.
This function will calculate the Brunt-Vaisala frequency as follows:
\[N = \left( \frac{g}{\theta} \frac{d\theta}{dz} \right)^\frac{1}{2}\]This formula based off of Equations 3.75 and 3.77 in [Hobbs2006].
This function is a wrapper for
brunt_vaisala_frequency_squared
that filters out negative (unstable) quantities and takes the square root.- Parameters
height (
xarray.DataArray
orpint.Quantity
) – Atmospheric (geopotential) heightpotential_temperature (
xarray.DataArray
orpint.Quantity
) – Atmospheric potential temperaturevertical_dim (int, optional) – The axis corresponding to vertical in the potential temperature array, defaults to 0, unless height and
potential_temperature
given asxarray.DataArray
, in which case it is automatically determined from the coordinate metadata.
- Returns
pint.Quantity
orxarray.DataArray
– Brunt-Vaisala frequency. Given aspint.Quantity
, unless both height andpotential_temperature
arguments are given asxarray.DataArray
, in which case will bexarray.DataArray
.
Changed in version 1.0: Renamed
heights
,axis
parameters toheight
,vertical_dim