profile_standard#

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

Standard height profiles according to ITU-R P.835-5, Annex 1.

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_standard(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: 288.1 K 1013.2 hPa
10 km: 223.1 K  264.4 hPa
20 km: 216.6 K   54.7 hPa
30 km: 226.6 K   11.7 hPa
40 km: 251.0 K    2.8 hPa
50 km: 270.6 K    0.8 hPa
60 km: 245.4 K    0.2 hPa
70 km: 217.4 K    0.0 hPa
80 km: 196.6 K    0.0 hPa