profile_highlat_winter#

pycraf.atm.profile_highlat_winter(height)[source]#

High latitude winter height profiles according to ITU-R P.835-5.

Valid for geographic latitudes \(\vert \phi\vert > 45^\circ\).

Parameters:

heightQuantity: Height above ground [km]

Returns:

temperatureQuantity: Temperature [K]
pressureQuantity: Total pressure [hPa]
rho_waterQuantity: Water vapor density [g / m**3]
pressure_waterQuantity: Water vapor partial pressure [hPa]
ref_indexQuantity: Refractive index [dimless]
humidity_waterQuantity: Relative humidity if water vapor was in form of liquid water [%]
humidity_iceQuantity: Relative humidity if water vapor was in form of ice [%]

Notes

For convenience, derived quantities like water density/pressure and refraction indices are also returned.

The return value is actually a namedtuple, so it is possible to do the following:

>>> import numpy as np
>>> from pycraf import atm
>>> from astropy import units as u

>>> heights = np.linspace(0, 80, 9) * u.km
>>> aprof = atm.profile_highlat_winter(heights)
>>> for height, temp, press in zip(
...         heights, aprof.temperature, aprof.pressure
...         ):
...     print('{:2.0f}: {:5.1f} {:6.1f}'.format(height, temp, press))
 0 km: 257.4 K 1010.9 hPa
10 km: 217.5 K  243.9 hPa
20 km: 217.5 K   56.1 hPa
30 km: 217.5 K   12.9 hPa
40 km: 238.8 K    3.0 hPa
50 km: 260.0 K    0.7 hPa
60 km: 250.0 K    0.2 hPa
70 km: 233.3 K    0.0 hPa
80 km: 216.7 K    0.0 hPa