earthkit.meteo.vertical.array.pressure_at_height_levels ======================================================= .. py:function:: 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. :param height: height above the surface for which the pressure needs to be computed (m) :type height: :class:`number` :param t: temperature at model full-levels (K). First dimension must correspond to the model full-levels. :type t: :class:`ndarray` :param q: specific humidity at model full-levels (kg/kg). First dimension must correspond to the model full-levels. :type q: :class:`ndarray` :param sp: surface pressure (Pa) :type sp: :class:`ndarray` :param A: A-coefficients defining the model levels :type A: :class:`ndarray` :param B: B-coefficients defining the model levels :type B: :class:`ndarray` :param alpha_top: Option passed to :func:`pressure_at_model_levels`. The possible values are: "ifs" (default) or "arpege". :type alpha_top: :class:`str`, *optional* :returns: pressure at the given height level (Pa) :rtype: :class:`number` or :class:`ndarray` .. rubric:: 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. .. seealso:: :obj:`pressure_at_model_levels`, :obj:`relative_geopotential_thickness`