Stability

Vertical stability functions.

These work with ndarrays of profiles; use the axis keyword argument to specify the axis along which pressure varies. For example, the default, following the Matlab versions, is axis=0, meaning the pressure varies along the first dimension. Use axis=-1 if pressure varies along the last dimension–that is, along a row, as the column index increases, in the 2-D case.

Docstrings will be added later, either manually or via an automated mechanism.

gsw.stability.IPV_vs_fNsquared_ratio(SA, CT, p, p_ref=0, axis=0)[source]

Calculates the ratio of the vertical gradient of potential density to the vertical gradient of locally-referenced potential density. This is also the ratio of the planetary Isopycnal Potential Vorticity (IPV) to f times N^2, hence the name for this variable, IPV_vs_fNsquared_ratio (see Eqn. (3.20.17) of IOC et al. (2010)).

Parameters:
SAarray-like

Absolute Salinity, g/kg

tarray-like

In-situ temperature (ITS-90), degrees C

parray-like

Sea pressure (absolute pressure minus 10.1325 dbar), dbar

p_reffloat

Reference pressure, dbar

Returns:
IPV_vs_fNsquared_ratioarray

The ratio of the vertical gradient of potential density referenced to p_ref, to the vertical gradient of locally-referenced potential density, dimensionless.

p_midarray

Pressure at midpoints of p, dbar. The array shape matches IPV_vs_fNsquared_ratio.

gsw.stability.Nsquared(SA, CT, p, lat=None, axis=0)[source]

Calculate the square of the buoyancy frequency.

Parameters:
SAarray-like

Absolute Salinity, g/kg

CTarray-like

Conservative Temperature (ITS-90), degrees C

parray-like

Sea pressure (absolute pressure minus 10.1325 dbar), dbar

latarray-like, 1-D, optional

Latitude, degrees.

axisint, optional

The dimension along which pressure increases.

Returns:
N2array

Buoyancy frequency-squared at pressure midpoints, 1/s^2. The shape along the pressure axis dimension is one less than that of the inputs. (Frequency N is in radians per second.)

p_midarray

Pressure at midpoints of p, dbar. The array shape matches N2.

gsw.stability.Turner_Rsubrho(SA, CT, p, axis=0)[source]

Calculate the Turner Angle and the Stability Ratio.

Parameters:
SAarray-like

Absolute Salinity, g/kg

CTarray-like

Conservative Temperature (ITS-90), degrees C

parray-like

Sea pressure (absolute pressure minus 10.1325 dbar), dbar

axisint, optional

The dimension along which pressure increases.

Returns:
Tuarray

Turner Angle at pressure midpoints, degrees. The shape along the pressure axis dimension is one less than that of the inputs.

Rsubrhoarray

Stability Ratio, dimensionless. The shape matches Tu.

p_midarray

Pressure at midpoints of p, dbar. The array shape matches Tu.