profile_midlat_winter#

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

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

Valid for geographic latitudes \(22^\circ < \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_midlat_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: 272.7 K 1018.9 hPa
10 km: 218.0 K  259.0 hPa
20 km: 218.0 K   59.5 hPa
30 km: 218.0 K   13.7 hPa
40 km: 241.5 K    3.1 hPa
50 km: 265.0 K    0.7 hPa
60 km: 250.7 K    0.2 hPa
70 km: 230.4 K    0.0 hPa
80 km: 210.0 K    0.0 hPa