earthkit.meteo.vertical.array.pressure_at_height_levels

earthkit.meteo.vertical.array.pressure_at_height_levels(height, t, q, sp, A, B, alpha_top='ifs')

Calculate the pressure at a height above the surface from model full-levels.

Deprecated in version 1.0.0 See pressure_at_height_levels() is deprecated for details.

Parameters:

  • height (number) – height above the surface for which the pressure needs to be computed (m)

  • t (ndarray) – temperature at model full-levels (K). First dimension must correspond to the model full-levels.

  • q (ndarray) – specific humidity at model full-levels (kg/kg). First dimension must correspond to the model full-levels.

  • sp (ndarray) – surface pressure (Pa)

  • A (ndarray) – A-coefficients defining the model levels

  • B (ndarray) – B-coefficients defining the model levels

  • alpha_top (str, optional) – Option passed to pressure_at_model_levels(). The possible values are: “ifs” (default) or “arpege”.

Returns:

pressure at the given height level (Pa)

Return type:

number or ndarray

See also

pressure_at_model_levels, relative_geopotential_thickness

Notes

t and q must contain the same model levels in ascending order with respect to the model level number. The model level range must be contiguous and must include the bottom-most level, but not all the levels must be present. E.g. if the vertical coordinate system has 137 model levels using only a subset of levels between e.g. 137-96 is allowed.

A and B must be defined on the model half-levels corresponding to the model full-levels in t and q. So the number of levels in A and B must be one more than the number of levels in t and q.

The pressure at height level is calculated by finding the model level above and below the specified height and interpolating the pressure with linear interpolation.