konrad.utils.get_quadratic_pgrid

konrad.utils.get_quadratic_pgrid(surface_pressure=100000.0, top_pressure=1, num=128)[source]

Create matching pressure levels and half-levels.

The half-levels range from surface_pressure to 1 Pa.

The pressure for every level i ([0, N]) is given by:

\[ln(p/p_\mathrm{t}) = -\frac{\ln(p_\mathrm{s}/p_\mathrm{t})}{2} \left(\frac{i^2}{N^2} + \frac{i}{N}\right) + \ln(p_\mathrm{s}/p_\mathrm{t})\]
Parameters

  • surface_pressure (float) – Pressure of the lowest first grid point [Pa].

  • top_pressure (float) – Pressure at the highest last grid point [Pa].

  • num (int) – Number of grid points.

Returns

ndarray – Pressure grid [Pa].