p835 0.1.0__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- p835/__init__.py +87 -0
- p835/common/__init__.py +7 -0
- p835/common/constants.py +172 -0
- p835/common/exceptions.py +40 -0
- p835/conversions.py +106 -0
- p835/reference_atmosphere.py +328 -0
- p835/seasonal_atmosphere.py +499 -0
- p835-0.1.0.dist-info/METADATA +70 -0
- p835-0.1.0.dist-info/RECORD +13 -0
- p835-0.1.0.dist-info/WHEEL +5 -0
- p835-0.1.0.dist-info/licenses/LICENSE +21 -0
- p835-0.1.0.dist-info/licenses/NOTICE +26 -0
- p835-0.1.0.dist-info/top_level.txt +1 -0
p835/__init__.py
ADDED
|
@@ -0,0 +1,87 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
p835 — ITU-R P.835-7 (Reference standard atmospheres).
|
|
9
|
+
|
|
10
|
+
Использование:
|
|
11
|
+
from p835 import global_temperature, global_pressure
|
|
12
|
+
|
|
13
|
+
T = global_temperature(10.0) # K
|
|
14
|
+
P = global_pressure(10.0) # гПа
|
|
15
|
+
"""
|
|
16
|
+
|
|
17
|
+
from .reference_atmosphere import (
|
|
18
|
+
global_temperature,
|
|
19
|
+
global_pressure,
|
|
20
|
+
global_water_vapour_density,
|
|
21
|
+
global_wet_pressure,
|
|
22
|
+
MODEL_NAME,
|
|
23
|
+
)
|
|
24
|
+
from .seasonal_atmosphere import (
|
|
25
|
+
AtmosphereRegion,
|
|
26
|
+
Season,
|
|
27
|
+
interpolated_temperature,
|
|
28
|
+
interpolated_pressure,
|
|
29
|
+
interpolated_water_vapour,
|
|
30
|
+
get_region_for_latitude,
|
|
31
|
+
regional_temperature,
|
|
32
|
+
regional_pressure,
|
|
33
|
+
regional_water_vapour,
|
|
34
|
+
low_latitude_annual_temperature,
|
|
35
|
+
low_latitude_annual_pressure,
|
|
36
|
+
low_latitude_annual_water_vapour,
|
|
37
|
+
mid_latitude_summer_temperature,
|
|
38
|
+
mid_latitude_summer_pressure,
|
|
39
|
+
mid_latitude_summer_water_vapour,
|
|
40
|
+
mid_latitude_winter_temperature,
|
|
41
|
+
mid_latitude_winter_pressure,
|
|
42
|
+
mid_latitude_winter_water_vapour,
|
|
43
|
+
high_latitude_summer_temperature,
|
|
44
|
+
high_latitude_summer_pressure,
|
|
45
|
+
high_latitude_summer_water_vapour,
|
|
46
|
+
high_latitude_winter_temperature,
|
|
47
|
+
high_latitude_winter_pressure,
|
|
48
|
+
high_latitude_winter_water_vapour,
|
|
49
|
+
)
|
|
50
|
+
from .conversions import (
|
|
51
|
+
convert_to_geopotential_height,
|
|
52
|
+
convert_to_geometric_height,
|
|
53
|
+
water_vapour_density_to_pressure,
|
|
54
|
+
water_vapour_pressure_to_density,
|
|
55
|
+
)
|
|
56
|
+
|
|
57
|
+
__version__ = "0.1.0"
|
|
58
|
+
|
|
59
|
+
__all__ = [
|
|
60
|
+
"global_temperature", "global_pressure",
|
|
61
|
+
"global_water_vapour_density", "global_wet_pressure",
|
|
62
|
+
"MODEL_NAME",
|
|
63
|
+
"AtmosphereRegion", "Season",
|
|
64
|
+
"interpolated_temperature", "interpolated_pressure",
|
|
65
|
+
"interpolated_water_vapour", "get_region_for_latitude",
|
|
66
|
+
"regional_temperature", "regional_pressure",
|
|
67
|
+
"regional_water_vapour",
|
|
68
|
+
"low_latitude_annual_temperature",
|
|
69
|
+
"low_latitude_annual_pressure",
|
|
70
|
+
"low_latitude_annual_water_vapour",
|
|
71
|
+
"mid_latitude_summer_temperature",
|
|
72
|
+
"mid_latitude_summer_pressure",
|
|
73
|
+
"mid_latitude_summer_water_vapour",
|
|
74
|
+
"mid_latitude_winter_temperature",
|
|
75
|
+
"mid_latitude_winter_pressure",
|
|
76
|
+
"mid_latitude_winter_water_vapour",
|
|
77
|
+
"high_latitude_summer_temperature",
|
|
78
|
+
"high_latitude_summer_pressure",
|
|
79
|
+
"high_latitude_summer_water_vapour",
|
|
80
|
+
"high_latitude_winter_temperature",
|
|
81
|
+
"high_latitude_winter_pressure",
|
|
82
|
+
"high_latitude_winter_water_vapour",
|
|
83
|
+
"convert_to_geopotential_height",
|
|
84
|
+
"convert_to_geometric_height",
|
|
85
|
+
"water_vapour_density_to_pressure",
|
|
86
|
+
"water_vapour_pressure_to_density",
|
|
87
|
+
]
|
p835/common/__init__.py
ADDED
p835/common/constants.py
ADDED
|
@@ -0,0 +1,172 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
core/itu/common/constants.py — Физические и математические константы,
|
|
9
|
+
используемые во всех моделях ITU-R.
|
|
10
|
+
|
|
11
|
+
Источник: Рекомендации ITU-R P-серии, U.S. Standard Atmosphere 1976.
|
|
12
|
+
"""
|
|
13
|
+
|
|
14
|
+
import math
|
|
15
|
+
|
|
16
|
+
# =====================================================================
|
|
17
|
+
# ГЕОМЕТРИЧЕСКИЕ И ГРАВИТАЦИОННЫЕ КОНСТАНТЫ
|
|
18
|
+
# =====================================================================
|
|
19
|
+
|
|
20
|
+
# Радиус Земли (средний), км — ITU-R P.452, P.526, P.530, P.1812
|
|
21
|
+
EARTH_RADIUS_KM = 6371.0
|
|
22
|
+
|
|
23
|
+
# Радиус Земли по U.S. Standard Atmosphere 1976, км — P.835 (1a/1b)
|
|
24
|
+
# ВАЖНО: в P.835 используется именно это значение (6 356.766 км),
|
|
25
|
+
# а не 6371. Не путать!
|
|
26
|
+
EARTH_RADIUS_P835_KM = 6356.766
|
|
27
|
+
|
|
28
|
+
# Ускорение свободного падения на поверхности Земли, м/с²
|
|
29
|
+
G0 = 9.80665
|
|
30
|
+
|
|
31
|
+
# =====================================================================
|
|
32
|
+
# АТМОСФЕРНЫЕ КОНСТАНТЫ (стандартные, на уровне моря)
|
|
33
|
+
# =====================================================================
|
|
34
|
+
|
|
35
|
+
# Стандартные условия U.S. Standard Atmosphere 1976 / ITU-R P.835
|
|
36
|
+
T0_K = 288.15 # температура, K
|
|
37
|
+
P0_HPA = 1013.25 # давление, гПа
|
|
38
|
+
RHO0_WATER_G_M3 = 7.5 # плотность водяного пара, г/м³
|
|
39
|
+
|
|
40
|
+
# Универсальная газовая постоянная, Дж/(моль·К)
|
|
41
|
+
R_UNIVERSAL = 8.314462618
|
|
42
|
+
|
|
43
|
+
# Газовая постоянная сухого воздуха, Дж/(кг·К)
|
|
44
|
+
R_DRY = 287.053
|
|
45
|
+
|
|
46
|
+
# Газовая постоянная водяного пара, Дж/(кг·К)
|
|
47
|
+
R_H2O = 461.495
|
|
48
|
+
|
|
49
|
+
# Молярные массы (кг/моль)
|
|
50
|
+
M_DRY_AIR = 28.9644e-3
|
|
51
|
+
M_H2O = 18.0153e-3
|
|
52
|
+
|
|
53
|
+
# Молекулярные массы (относительные, для P.676)
|
|
54
|
+
# M_d = 28.966, M_w = 18.01528
|
|
55
|
+
M_DRY_AIR_P676 = 28.966
|
|
56
|
+
M_H2O_P676 = 18.01528
|
|
57
|
+
|
|
58
|
+
# =====================================================================
|
|
59
|
+
# МАТЕМАТИЧЕСКИЕ КОНСТАНТЫ
|
|
60
|
+
# =====================================================================
|
|
61
|
+
|
|
62
|
+
# Скорость света в вакууме, м/с
|
|
63
|
+
SPEED_OF_LIGHT = 299792458.0
|
|
64
|
+
|
|
65
|
+
# Скорость света, приближённая, км/с — для удобства формул (c ≈ 300 000)
|
|
66
|
+
SPEED_OF_LIGHT_KM_S = 300000.0
|
|
67
|
+
|
|
68
|
+
# Постоянная Больцмана, Дж/К
|
|
69
|
+
K_BOLTZMANN = 1.380649e-23
|
|
70
|
+
|
|
71
|
+
# Температура стандартного шума, K (ITU-R P.372)
|
|
72
|
+
T_NOISE_K = 290.0
|
|
73
|
+
|
|
74
|
+
# Постоянная Эйлера-Маскерони (не путать с числом e!)
|
|
75
|
+
EULER_GAMMA = 0.5772156649015329
|
|
76
|
+
|
|
77
|
+
# Число e
|
|
78
|
+
E = math.e
|
|
79
|
+
|
|
80
|
+
# π
|
|
81
|
+
PI = math.pi
|
|
82
|
+
|
|
83
|
+
# =====================================================================
|
|
84
|
+
# ПРЕДЕЛЫ ПРИМЕНИМОСТИ МОДЕЛЕЙ ITU-R
|
|
85
|
+
# (для валидации входных параметров)
|
|
86
|
+
# =====================================================================
|
|
87
|
+
|
|
88
|
+
# P.835: диапазон высот
|
|
89
|
+
P835_Z_MIN_KM = 0.0
|
|
90
|
+
P835_Z_MAX_KM = 100.0
|
|
91
|
+
P835_Z_REGIME1_MAX_KM = 86.0 # переход на regime 2
|
|
92
|
+
P835_H_REGIME1_MAX_KM = 84.852 # геопотенциальная граница
|
|
93
|
+
|
|
94
|
+
# P.676: диапазон частот
|
|
95
|
+
P676_F_MIN_GHZ = 1.0
|
|
96
|
+
P676_F_MAX_GHZ = 1000.0
|
|
97
|
+
|
|
98
|
+
# P.528: диапазон частот, расстояний, высот
|
|
99
|
+
P528_F_MIN_MHZ = 100.0
|
|
100
|
+
P528_F_MAX_MHZ = 30000.0
|
|
101
|
+
P528_H_MIN_M = 1.5
|
|
102
|
+
P528_H_MAX_M = 20000.0
|
|
103
|
+
P528_D_MAX_KM = 1800.0
|
|
104
|
+
|
|
105
|
+
# P.530: диапазон частот, расстояний
|
|
106
|
+
P530_F_MIN_GHZ = 0.45
|
|
107
|
+
P530_F_MAX_GHZ = 45.0
|
|
108
|
+
P530_D_MIN_KM = 7.5
|
|
109
|
+
P530_D_MAX_KM = 300.0
|
|
110
|
+
|
|
111
|
+
# P.452: диапазон частот, расстояний
|
|
112
|
+
P452_F_MIN_GHZ = 0.1
|
|
113
|
+
P452_F_MAX_GHZ = 50.0
|
|
114
|
+
P452_D_MAX_KM = 10000.0
|
|
115
|
+
|
|
116
|
+
# P.1812: диапазон частот, расстояний
|
|
117
|
+
P1812_F_MIN_MHZ = 30.0
|
|
118
|
+
P1812_F_MAX_MHZ = 6000.0
|
|
119
|
+
P1812_D_MIN_KM = 0.25
|
|
120
|
+
P1812_D_MAX_KM = 3000.0
|
|
121
|
+
|
|
122
|
+
# P.2108: диапазон частот
|
|
123
|
+
P2108_F_MIN_GHZ = 0.03
|
|
124
|
+
P2108_F_MAX_GHZ = 100.0
|
|
125
|
+
|
|
126
|
+
# P.533: диапазон частот
|
|
127
|
+
P533_F_MIN_MHZ = 2.0
|
|
128
|
+
P533_F_MAX_MHZ = 30.0
|
|
129
|
+
|
|
130
|
+
# =====================================================================
|
|
131
|
+
# ЗНАЧЕНИЯ ПО УМОЛЧАНИЮ
|
|
132
|
+
# =====================================================================
|
|
133
|
+
|
|
134
|
+
DEFAULT_TEMPERATURE_C = 15.0
|
|
135
|
+
DEFAULT_PRESSURE_HPA = 1013.25
|
|
136
|
+
DEFAULT_HUMIDITY_PERCENT = 50.0
|
|
137
|
+
DEFAULT_WATER_VAPOUR_DENSITY_G_M3 = 7.5
|
|
138
|
+
DEFAULT_POLARIZATION = 'V'
|
|
139
|
+
|
|
140
|
+
# __all__ — публичный API модуля
|
|
141
|
+
__all__ = [
|
|
142
|
+
# Геометрия и гравитация
|
|
143
|
+
'EARTH_RADIUS_KM',
|
|
144
|
+
'EARTH_RADIUS_P835_KM',
|
|
145
|
+
'G0',
|
|
146
|
+
# Атмосфера
|
|
147
|
+
'T0_K', 'P0_HPA', 'RHO0_WATER_G_M3',
|
|
148
|
+
'R_UNIVERSAL', 'R_DRY', 'R_H2O',
|
|
149
|
+
'M_DRY_AIR', 'M_H2O',
|
|
150
|
+
'M_DRY_AIR_P676', 'M_H2O_P676',
|
|
151
|
+
# Математика
|
|
152
|
+
'SPEED_OF_LIGHT', 'SPEED_OF_LIGHT_KM_S',
|
|
153
|
+
'K_BOLTZMANN', 'T_NOISE_K',
|
|
154
|
+
'EULER_GAMMA', 'E', 'PI',
|
|
155
|
+
# Пределы
|
|
156
|
+
'P835_Z_MIN_KM', 'P835_Z_MAX_KM',
|
|
157
|
+
'P835_Z_REGIME1_MAX_KM', 'P835_H_REGIME1_MAX_KM',
|
|
158
|
+
'P676_F_MIN_GHZ', 'P676_F_MAX_GHZ',
|
|
159
|
+
'P528_F_MIN_MHZ', 'P528_F_MAX_MHZ',
|
|
160
|
+
'P528_H_MIN_M', 'P528_H_MAX_M', 'P528_D_MAX_KM',
|
|
161
|
+
'P530_F_MIN_GHZ', 'P530_F_MAX_GHZ',
|
|
162
|
+
'P530_D_MIN_KM', 'P530_D_MAX_KM',
|
|
163
|
+
'P452_F_MIN_GHZ', 'P452_F_MAX_GHZ', 'P452_D_MAX_KM',
|
|
164
|
+
'P1812_F_MIN_MHZ', 'P1812_F_MAX_MHZ',
|
|
165
|
+
'P1812_D_MIN_KM', 'P1812_D_MAX_KM',
|
|
166
|
+
'P2108_F_MIN_GHZ', 'P2108_F_MAX_GHZ',
|
|
167
|
+
'P533_F_MIN_MHZ', 'P533_F_MAX_MHZ',
|
|
168
|
+
# Умолчания
|
|
169
|
+
'DEFAULT_TEMPERATURE_C', 'DEFAULT_PRESSURE_HPA',
|
|
170
|
+
'DEFAULT_HUMIDITY_PERCENT', 'DEFAULT_WATER_VAPOUR_DENSITY_G_M3',
|
|
171
|
+
'DEFAULT_POLARIZATION',
|
|
172
|
+
]
|
|
@@ -0,0 +1,40 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
core/itu/common/exceptions.py — Исключения для моделей ITU-R.
|
|
9
|
+
"""
|
|
10
|
+
|
|
11
|
+
|
|
12
|
+
class ITUError(Exception):
|
|
13
|
+
"""Базовый класс для всех ошибок моделей ITU-R."""
|
|
14
|
+
pass
|
|
15
|
+
|
|
16
|
+
|
|
17
|
+
class ITUValidationError(ITUError):
|
|
18
|
+
"""
|
|
19
|
+
Ошибка валидации входных параметров.
|
|
20
|
+
|
|
21
|
+
Поднимается, когда параметр критически вне диапазона применимости
|
|
22
|
+
модели, и клиент запросил строгий режим (strict=True).
|
|
23
|
+
"""
|
|
24
|
+
|
|
25
|
+
def __init__(self, model: str, param: str, value, expected: str):
|
|
26
|
+
self.model = model
|
|
27
|
+
self.param = param
|
|
28
|
+
self.value = value
|
|
29
|
+
self.expected = expected
|
|
30
|
+
super().__init__(
|
|
31
|
+
f"[{model}] Параметр '{param}' = {value} вне диапазона: {expected}"
|
|
32
|
+
)
|
|
33
|
+
|
|
34
|
+
|
|
35
|
+
class ITUComputationError(ITUError):
|
|
36
|
+
"""Ошибка в процессе вычисления (NaN, деление на ноль и т.п.)."""
|
|
37
|
+
pass
|
|
38
|
+
|
|
39
|
+
|
|
40
|
+
__all__ = ['ITUError', 'ITUValidationError', 'ITUComputationError']
|
p835/conversions.py
ADDED
|
@@ -0,0 +1,106 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
core/itu/p835/conversions.py — Преобразования высот и давлений
|
|
9
|
+
по ITU-R P.835-7 (Annex 1).
|
|
10
|
+
|
|
11
|
+
Формулы:
|
|
12
|
+
(1a): H = 6356.766·Z / (6356.766 + Z)
|
|
13
|
+
(1b): Z = 6356.766·H / (6356.766 − H)
|
|
14
|
+
(7): e = ρ·T / 216.7
|
|
15
|
+
"""
|
|
16
|
+
|
|
17
|
+
from .common.constants import EARTH_RADIUS_P835_KM
|
|
18
|
+
|
|
19
|
+
|
|
20
|
+
def convert_to_geopotential_height(z_km: float) -> float:
|
|
21
|
+
"""
|
|
22
|
+
Геометрическая высота Z (км) → геопотенциальная H (км).
|
|
23
|
+
|
|
24
|
+
P.835-7, ур. (1a):
|
|
25
|
+
H = 6356.766·Z / (6356.766 + Z)
|
|
26
|
+
|
|
27
|
+
Args:
|
|
28
|
+
z_km: геометрическая высота, км (0..100).
|
|
29
|
+
|
|
30
|
+
Returns:
|
|
31
|
+
Геопотенциальная высота, км.
|
|
32
|
+
"""
|
|
33
|
+
r = EARTH_RADIUS_P835_KM
|
|
34
|
+
return (r * z_km) / (r + z_km)
|
|
35
|
+
|
|
36
|
+
|
|
37
|
+
def convert_to_geometric_height(h_km: float) -> float:
|
|
38
|
+
"""
|
|
39
|
+
Геопотенциальная высота H (км) → геометрическая Z (км).
|
|
40
|
+
|
|
41
|
+
P.835-7, ур. (1b):
|
|
42
|
+
Z = 6356.766·H / (6356.766 − H)
|
|
43
|
+
|
|
44
|
+
Args:
|
|
45
|
+
h_km: геопотенциальная высота, км (0..84.852).
|
|
46
|
+
|
|
47
|
+
Returns:
|
|
48
|
+
Геометрическая высота, км.
|
|
49
|
+
|
|
50
|
+
Raises:
|
|
51
|
+
ValueError: если h_km >= 6356.766 (вырожденный случай).
|
|
52
|
+
"""
|
|
53
|
+
r = EARTH_RADIUS_P835_KM
|
|
54
|
+
if h_km >= r:
|
|
55
|
+
raise ValueError(
|
|
56
|
+
f"convert_to_geometric_height: H={h_km} км ≥ {r} км"
|
|
57
|
+
)
|
|
58
|
+
return (r * h_km) / (r - h_km)
|
|
59
|
+
|
|
60
|
+
|
|
61
|
+
def water_vapour_density_to_pressure(
|
|
62
|
+
rho_g_m3: float, temperature_k: float
|
|
63
|
+
) -> float:
|
|
64
|
+
"""
|
|
65
|
+
Плотность водяного пара ρ (г/м³) → парциальное давление e (гПа).
|
|
66
|
+
|
|
67
|
+
P.835-7, ур. (7):
|
|
68
|
+
e = ρ·T / 216.7
|
|
69
|
+
|
|
70
|
+
Args:
|
|
71
|
+
rho_g_m3: плотность водяного пара, г/м³.
|
|
72
|
+
temperature_k: температура, K.
|
|
73
|
+
|
|
74
|
+
Returns:
|
|
75
|
+
Парциальное давление водяного пара, гПа.
|
|
76
|
+
"""
|
|
77
|
+
return rho_g_m3 * temperature_k / 216.7
|
|
78
|
+
|
|
79
|
+
|
|
80
|
+
def water_vapour_pressure_to_density(
|
|
81
|
+
e_hpa: float, temperature_k: float
|
|
82
|
+
) -> float:
|
|
83
|
+
"""
|
|
84
|
+
Парциальное давление e (гПа) → плотность ρ (г/м³).
|
|
85
|
+
|
|
86
|
+
Обратная к (7):
|
|
87
|
+
ρ = e·216.7 / T
|
|
88
|
+
|
|
89
|
+
Args:
|
|
90
|
+
e_hpa: парциальное давление, гПа.
|
|
91
|
+
temperature_k: температура, K.
|
|
92
|
+
|
|
93
|
+
Returns:
|
|
94
|
+
Плотность водяного пара, г/м³.
|
|
95
|
+
"""
|
|
96
|
+
if temperature_k <= 0:
|
|
97
|
+
return 0.0
|
|
98
|
+
return e_hpa * 216.7 / temperature_k
|
|
99
|
+
|
|
100
|
+
|
|
101
|
+
__all__ = [
|
|
102
|
+
'convert_to_geopotential_height',
|
|
103
|
+
'convert_to_geometric_height',
|
|
104
|
+
'water_vapour_density_to_pressure',
|
|
105
|
+
'water_vapour_pressure_to_density',
|
|
106
|
+
]
|
|
@@ -0,0 +1,328 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
core/itu/p835/reference_atmosphere.py — ITU-R reference atmosphere
|
|
9
|
+
по ITU-R P.835-7, Annex 1.
|
|
10
|
+
|
|
11
|
+
Реализует:
|
|
12
|
+
* GlobalTemperature(Z) — температура, K
|
|
13
|
+
* GlobalPressure(Z) — давление, гПа
|
|
14
|
+
* GlobalWaterVapourDensity(Z, rho_0) — плотность H₂O, г/м³
|
|
15
|
+
* GlobalWetPressure(Z, rho_0) — парциальное давление H₂O, гПа
|
|
16
|
+
|
|
17
|
+
Диапазон высот: Z ∈ [0, 100] км.
|
|
18
|
+
|
|
19
|
+
Два режима (P.835-7, Annex 1 §1.1):
|
|
20
|
+
* Regime 1: Z ≤ 86 км → вычисляется через геопотенциальную высоту H.
|
|
21
|
+
* Regime 2: 86 < Z ≤ 100 км → вычисляется напрямую по Z.
|
|
22
|
+
|
|
23
|
+
ВАЖНО: в документе P.835-7 (08/2024) формула (4b) для T(Z) в regime 2
|
|
24
|
+
содержит опечатку OCR: T = 263.1905 − 76.3232·[1 − ((Z−91)/19.9429)^½]².
|
|
25
|
+
Правильная формула — с квадратом (без корня) внутри скобок:
|
|
26
|
+
T = 263.1905 − 76.3232·√(1 − ((Z−91)/19.9429)²)
|
|
27
|
+
Это классическая U.S. Standard Atmosphere 1976.
|
|
28
|
+
"""
|
|
29
|
+
|
|
30
|
+
import math
|
|
31
|
+
|
|
32
|
+
from .common.constants import (
|
|
33
|
+
P835_Z_MIN_KM,
|
|
34
|
+
P835_Z_MAX_KM,
|
|
35
|
+
P835_Z_REGIME1_MAX_KM,
|
|
36
|
+
P835_H_REGIME1_MAX_KM,
|
|
37
|
+
RHO0_WATER_G_M3,
|
|
38
|
+
)
|
|
39
|
+
from .common.exceptions import ITUValidationError
|
|
40
|
+
from .conversions import (
|
|
41
|
+
convert_to_geopotential_height,
|
|
42
|
+
water_vapour_density_to_pressure,
|
|
43
|
+
)
|
|
44
|
+
|
|
45
|
+
|
|
46
|
+
MODEL_NAME = "ITU-R P.835-7"
|
|
47
|
+
|
|
48
|
+
|
|
49
|
+
# =====================================================================
|
|
50
|
+
# ТЕМПЕРАТУРА
|
|
51
|
+
# =====================================================================
|
|
52
|
+
|
|
53
|
+
def _global_temperature_regime1(h_prime_km: float) -> float:
|
|
54
|
+
"""
|
|
55
|
+
Температура T(H) для geopotential height H ≤ 84.852 км.
|
|
56
|
+
P.835-7, ур. (2a–2g).
|
|
57
|
+
|
|
58
|
+
Args:
|
|
59
|
+
h_prime_km: геопотенциальная высота, км.
|
|
60
|
+
|
|
61
|
+
Returns:
|
|
62
|
+
Температура, K.
|
|
63
|
+
"""
|
|
64
|
+
if h_prime_km < 0.0:
|
|
65
|
+
raise ITUValidationError(
|
|
66
|
+
MODEL_NAME, "H'", h_prime_km, ">= 0"
|
|
67
|
+
)
|
|
68
|
+
|
|
69
|
+
if h_prime_km <= 11.0:
|
|
70
|
+
return 288.15 - 6.5 * h_prime_km
|
|
71
|
+
elif h_prime_km <= 20.0:
|
|
72
|
+
return 216.65
|
|
73
|
+
elif h_prime_km <= 32.0:
|
|
74
|
+
return 216.65 + (h_prime_km - 20.0)
|
|
75
|
+
elif h_prime_km <= 47.0:
|
|
76
|
+
return 228.65 + 2.8 * (h_prime_km - 32.0)
|
|
77
|
+
elif h_prime_km <= 51.0:
|
|
78
|
+
return 270.65
|
|
79
|
+
elif h_prime_km <= 71.0:
|
|
80
|
+
return 270.65 - 2.8 * (h_prime_km - 51.0)
|
|
81
|
+
elif h_prime_km <= P835_H_REGIME1_MAX_KM:
|
|
82
|
+
return 214.65 - 2.0 * (h_prime_km - 71.0)
|
|
83
|
+
else:
|
|
84
|
+
raise ITUValidationError(
|
|
85
|
+
MODEL_NAME, "H'", h_prime_km,
|
|
86
|
+
f"<= {P835_H_REGIME1_MAX_KM}"
|
|
87
|
+
)
|
|
88
|
+
|
|
89
|
+
|
|
90
|
+
def _global_temperature_regime2(z_km: float) -> float:
|
|
91
|
+
"""
|
|
92
|
+
Температура T(Z) для geometric height 86 ≤ Z ≤ 100 км.
|
|
93
|
+
P.835-7, ур. (4a, 4b).
|
|
94
|
+
|
|
95
|
+
Ур. (4a): T = 186.8673 для 86 < Z ≤ 91
|
|
96
|
+
Ур. (4b): T = 263.1905 − 76.3232·√(1 − ((Z−91)/19.9429)²) для 91 < Z ≤ 100
|
|
97
|
+
|
|
98
|
+
Args:
|
|
99
|
+
z_km: геометрическая высота, км.
|
|
100
|
+
|
|
101
|
+
Returns:
|
|
102
|
+
Температура, K.
|
|
103
|
+
"""
|
|
104
|
+
if z_km < P835_Z_REGIME1_MAX_KM:
|
|
105
|
+
raise ITUValidationError(
|
|
106
|
+
MODEL_NAME, "Z", z_km, f">= {P835_Z_REGIME1_MAX_KM}"
|
|
107
|
+
)
|
|
108
|
+
if z_km > P835_Z_MAX_KM:
|
|
109
|
+
raise ITUValidationError(
|
|
110
|
+
MODEL_NAME, "Z", z_km, f"<= {P835_Z_MAX_KM}"
|
|
111
|
+
)
|
|
112
|
+
|
|
113
|
+
if z_km <= 91.0:
|
|
114
|
+
return 186.8673
|
|
115
|
+
|
|
116
|
+
ratio = (z_km - 91.0) / 19.9429
|
|
117
|
+
arg = 1.0 - ratio * ratio
|
|
118
|
+
if arg < 0.0:
|
|
119
|
+
arg = 0.0
|
|
120
|
+
return 263.1905 - 76.3232 * math.sqrt(arg)
|
|
121
|
+
|
|
122
|
+
|
|
123
|
+
def global_temperature(z_km: float) -> float:
|
|
124
|
+
"""
|
|
125
|
+
Температура T(Z) на высоте Z, K.
|
|
126
|
+
|
|
127
|
+
Args:
|
|
128
|
+
z_km: геометрическая высота, км (0 ≤ Z ≤ 100).
|
|
129
|
+
|
|
130
|
+
Returns:
|
|
131
|
+
Температура, K.
|
|
132
|
+
|
|
133
|
+
Raises:
|
|
134
|
+
ITUValidationError: если Z вне [0, 100].
|
|
135
|
+
"""
|
|
136
|
+
if z_km < P835_Z_MIN_KM or z_km > P835_Z_MAX_KM:
|
|
137
|
+
raise ITUValidationError(
|
|
138
|
+
MODEL_NAME, "Z", z_km,
|
|
139
|
+
f"[{P835_Z_MIN_KM}, {P835_Z_MAX_KM}]"
|
|
140
|
+
)
|
|
141
|
+
|
|
142
|
+
if z_km < P835_Z_REGIME1_MAX_KM:
|
|
143
|
+
h_prime = convert_to_geopotential_height(z_km)
|
|
144
|
+
return _global_temperature_regime1(h_prime)
|
|
145
|
+
else:
|
|
146
|
+
return _global_temperature_regime2(z_km)
|
|
147
|
+
|
|
148
|
+
|
|
149
|
+
# =====================================================================
|
|
150
|
+
# ДАВЛЕНИЕ
|
|
151
|
+
# =====================================================================
|
|
152
|
+
|
|
153
|
+
def _global_pressure_regime1(h_prime_km: float) -> float:
|
|
154
|
+
"""
|
|
155
|
+
Давление P(H) для geopotential height H ≤ 84.852 км.
|
|
156
|
+
P.835-7, ур. (3a–3g).
|
|
157
|
+
|
|
158
|
+
Args:
|
|
159
|
+
h_prime_km: геопотенциальная высота, км.
|
|
160
|
+
|
|
161
|
+
Returns:
|
|
162
|
+
Давление, гПа.
|
|
163
|
+
"""
|
|
164
|
+
if h_prime_km < 0.0:
|
|
165
|
+
raise ITUValidationError(
|
|
166
|
+
MODEL_NAME, "H'", h_prime_km, ">= 0"
|
|
167
|
+
)
|
|
168
|
+
|
|
169
|
+
if h_prime_km <= 11.0:
|
|
170
|
+
return 1013.25 * (288.15 / (288.15 - 6.5 * h_prime_km)) ** (-34.1632 / 6.5)
|
|
171
|
+
elif h_prime_km <= 20.0:
|
|
172
|
+
return 226.3226 * math.exp(-34.1632 * (h_prime_km - 11.0) / 216.65)
|
|
173
|
+
elif h_prime_km <= 32.0:
|
|
174
|
+
return 54.74980 * (216.65 / (216.65 + (h_prime_km - 20.0))) ** 34.1632
|
|
175
|
+
elif h_prime_km <= 47.0:
|
|
176
|
+
return 8.680422 * (228.65 / (228.65 + 2.8 * (h_prime_km - 32.0))) ** (34.1632 / 2.8)
|
|
177
|
+
elif h_prime_km <= 51.0:
|
|
178
|
+
return 1.109106 * math.exp(-34.1632 * (h_prime_km - 47.0) / 270.65)
|
|
179
|
+
elif h_prime_km <= 71.0:
|
|
180
|
+
return 0.6694167 * (270.65 / (270.65 - 2.8 * (h_prime_km - 51.0))) ** (-34.1632 / 2.8)
|
|
181
|
+
elif h_prime_km <= P835_H_REGIME1_MAX_KM:
|
|
182
|
+
return 0.03956649 * (214.65 / (214.65 - 2.0 * (h_prime_km - 71.0))) ** (-34.1632 / 2.0)
|
|
183
|
+
else:
|
|
184
|
+
raise ITUValidationError(
|
|
185
|
+
MODEL_NAME, "H'", h_prime_km,
|
|
186
|
+
f"<= {P835_H_REGIME1_MAX_KM}"
|
|
187
|
+
)
|
|
188
|
+
|
|
189
|
+
|
|
190
|
+
def _global_pressure_regime2(z_km: float) -> float:
|
|
191
|
+
"""
|
|
192
|
+
Давление P(Z) для geometric height 86 ≤ Z ≤ 100 км.
|
|
193
|
+
P.835-7, ур. (5).
|
|
194
|
+
|
|
195
|
+
Args:
|
|
196
|
+
z_km: геометрическая высота, км.
|
|
197
|
+
|
|
198
|
+
Returns:
|
|
199
|
+
Давление, гПа.
|
|
200
|
+
"""
|
|
201
|
+
if z_km < P835_Z_REGIME1_MAX_KM:
|
|
202
|
+
raise ITUValidationError(
|
|
203
|
+
MODEL_NAME, "Z", z_km, f">= {P835_Z_REGIME1_MAX_KM}"
|
|
204
|
+
)
|
|
205
|
+
if z_km > P835_Z_MAX_KM:
|
|
206
|
+
raise ITUValidationError(
|
|
207
|
+
MODEL_NAME, "Z", z_km, f"<= {P835_Z_MAX_KM}"
|
|
208
|
+
)
|
|
209
|
+
|
|
210
|
+
a0 = 95.571899
|
|
211
|
+
a1 = -4.011801
|
|
212
|
+
a2 = 6.424731e-2
|
|
213
|
+
a3 = -4.789660e-4
|
|
214
|
+
a4 = 1.340543e-6
|
|
215
|
+
|
|
216
|
+
z = z_km
|
|
217
|
+
return math.exp(
|
|
218
|
+
a0 + a1 * z + a2 * z * z + a3 * z * z * z + a4 * z * z * z * z
|
|
219
|
+
)
|
|
220
|
+
|
|
221
|
+
|
|
222
|
+
def global_pressure(z_km: float) -> float:
|
|
223
|
+
"""
|
|
224
|
+
Давление P(Z) на высоте Z, гПа.
|
|
225
|
+
|
|
226
|
+
Args:
|
|
227
|
+
z_km: геометрическая высота, км (0 ≤ Z ≤ 100).
|
|
228
|
+
|
|
229
|
+
Returns:
|
|
230
|
+
Давление, гПа.
|
|
231
|
+
|
|
232
|
+
Raises:
|
|
233
|
+
ITUValidationError: если Z вне [0, 100].
|
|
234
|
+
"""
|
|
235
|
+
if z_km < P835_Z_MIN_KM or z_km > P835_Z_MAX_KM:
|
|
236
|
+
raise ITUValidationError(
|
|
237
|
+
MODEL_NAME, "Z", z_km,
|
|
238
|
+
f"[{P835_Z_MIN_KM}, {P835_Z_MAX_KM}]"
|
|
239
|
+
)
|
|
240
|
+
|
|
241
|
+
if z_km < P835_Z_REGIME1_MAX_KM:
|
|
242
|
+
h_prime = convert_to_geopotential_height(z_km)
|
|
243
|
+
return _global_pressure_regime1(h_prime)
|
|
244
|
+
else:
|
|
245
|
+
return _global_pressure_regime2(z_km)
|
|
246
|
+
|
|
247
|
+
|
|
248
|
+
# =====================================================================
|
|
249
|
+
# ПЛОТНОСТЬ ВОДЯНОГО ПАРА
|
|
250
|
+
# =====================================================================
|
|
251
|
+
|
|
252
|
+
def global_water_vapour_density(
|
|
253
|
+
z_km: float, rho_0_g_m3: float = RHO0_WATER_G_M3
|
|
254
|
+
) -> float:
|
|
255
|
+
"""
|
|
256
|
+
Плотность водяного пара ρ(Z), г/м³.
|
|
257
|
+
|
|
258
|
+
P.835-7, ур. (6), (8):
|
|
259
|
+
ρ(Z) = 7.5·exp(−Z/2) (6)
|
|
260
|
+
ρ(Z) = 2e−6 · P(Z)·216.7 / T(Z) (выше точки смешивания) (8)
|
|
261
|
+
|
|
262
|
+
Точка смешивания: e(Z)/P(Z) = 2e−6, где e(Z) — парциальное
|
|
263
|
+
давление H₂O из ур. (7). Выше этой точки используется ур. (8).
|
|
264
|
+
|
|
265
|
+
Args:
|
|
266
|
+
z_km: геометрическая высота, км (0 ≤ Z ≤ 100).
|
|
267
|
+
rho_0_g_m3: плотность H₂O у поверхности, г/м³ (по умолчанию 7.5).
|
|
268
|
+
|
|
269
|
+
Returns:
|
|
270
|
+
Плотность водяного пара, г/м³.
|
|
271
|
+
|
|
272
|
+
Raises:
|
|
273
|
+
ITUValidationError: если Z вне [0, 100].
|
|
274
|
+
"""
|
|
275
|
+
if z_km < P835_Z_MIN_KM or z_km > P835_Z_MAX_KM:
|
|
276
|
+
raise ITUValidationError(
|
|
277
|
+
MODEL_NAME, "Z", z_km,
|
|
278
|
+
f"[{P835_Z_MIN_KM}, {P835_Z_MAX_KM}]"
|
|
279
|
+
)
|
|
280
|
+
|
|
281
|
+
# Базовая экспонента (ур. 6)
|
|
282
|
+
rho_exp = rho_0_g_m3 * math.exp(-z_km / 2.0)
|
|
283
|
+
|
|
284
|
+
# Парциальное давление (ур. 7)
|
|
285
|
+
T = global_temperature(z_km)
|
|
286
|
+
e_hpa = water_vapour_density_to_pressure(rho_exp, T)
|
|
287
|
+
P_hpa = global_pressure(z_km)
|
|
288
|
+
|
|
289
|
+
# Проверяем mixing ratio
|
|
290
|
+
if P_hpa <= 0.0:
|
|
291
|
+
return 0.0
|
|
292
|
+
mixing_ratio = e_hpa / P_hpa
|
|
293
|
+
|
|
294
|
+
# Если mixing ratio ниже 2e-6 — используем ур. (8)
|
|
295
|
+
if mixing_ratio < 2.0e-6:
|
|
296
|
+
return 2.0e-6 * P_hpa * 216.7 / T
|
|
297
|
+
|
|
298
|
+
return rho_exp
|
|
299
|
+
|
|
300
|
+
|
|
301
|
+
def global_wet_pressure(
|
|
302
|
+
z_km: float, rho_0_g_m3: float = RHO0_WATER_G_M3
|
|
303
|
+
) -> float:
|
|
304
|
+
"""
|
|
305
|
+
Парциальное давление водяного пара e(Z), гПа.
|
|
306
|
+
|
|
307
|
+
P.835-7, ур. (7):
|
|
308
|
+
e(Z) = ρ(Z)·T(Z) / 216.7
|
|
309
|
+
|
|
310
|
+
Args:
|
|
311
|
+
z_km: геометрическая высота, км.
|
|
312
|
+
rho_0_g_m3: плотность H₂O у поверхности, г/м³.
|
|
313
|
+
|
|
314
|
+
Returns:
|
|
315
|
+
Парциальное давление, гПа.
|
|
316
|
+
"""
|
|
317
|
+
rho = global_water_vapour_density(z_km, rho_0_g_m3)
|
|
318
|
+
T = global_temperature(z_km)
|
|
319
|
+
return water_vapour_density_to_pressure(rho, T)
|
|
320
|
+
|
|
321
|
+
|
|
322
|
+
__all__ = [
|
|
323
|
+
'global_temperature',
|
|
324
|
+
'global_pressure',
|
|
325
|
+
'global_water_vapour_density',
|
|
326
|
+
'global_wet_pressure',
|
|
327
|
+
'MODEL_NAME',
|
|
328
|
+
]
|
|
@@ -0,0 +1,499 @@
|
|
|
1
|
+
# Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
2
|
+
# SPDX-License-Identifier: MIT
|
|
3
|
+
#
|
|
4
|
+
# Python implementation of ITU-R P.835-7.
|
|
5
|
+
# See README.md and NOTICE for attribution.
|
|
6
|
+
|
|
7
|
+
"""
|
|
8
|
+
core/itu/p835/seasonal_atmosphere.py — Сезонные эталонные атмосферы
|
|
9
|
+
по ITU-R P.835-7, Annex 2.
|
|
10
|
+
|
|
11
|
+
Реализует пять профилей:
|
|
12
|
+
1. Low latitude annual (15°N, все сезоны)
|
|
13
|
+
2. Mid latitude summer (45°N, лето)
|
|
14
|
+
3. Mid latitude winter (45°N, зима)
|
|
15
|
+
4. High latitude summer (60°N, лето)
|
|
16
|
+
5. High latitude winter (60°N, зима)
|
|
17
|
+
|
|
18
|
+
Интерполяция по широте:
|
|
19
|
+
0°–15°N: low annual
|
|
20
|
+
15°–45°N (лето/зима): линейная между low и mid
|
|
21
|
+
45°–60°N (лето/зима): линейная между mid и high
|
|
22
|
+
≥ 60°N: high
|
|
23
|
+
|
|
24
|
+
Все формулы дают T(K) и P(hPa) как функции Z(km), где Z —
|
|
25
|
+
геометрическая высота.
|
|
26
|
+
"""
|
|
27
|
+
|
|
28
|
+
import math
|
|
29
|
+
from enum import Enum
|
|
30
|
+
|
|
31
|
+
|
|
32
|
+
class AtmosphereRegion(Enum):
|
|
33
|
+
"""Тип региональной атмосферы."""
|
|
34
|
+
LOW_LATITUDE_ANNUAL = "low_annual"
|
|
35
|
+
MID_LATITUDE_SUMMER = "mid_summer"
|
|
36
|
+
MID_LATITUDE_WINTER = "mid_winter"
|
|
37
|
+
HIGH_LATITUDE_SUMMER = "high_summer"
|
|
38
|
+
HIGH_LATITUDE_WINTER = "high_winter"
|
|
39
|
+
|
|
40
|
+
|
|
41
|
+
class Season(Enum):
|
|
42
|
+
"""Сезон."""
|
|
43
|
+
SUMMER = "summer"
|
|
44
|
+
WINTER = "winter"
|
|
45
|
+
|
|
46
|
+
|
|
47
|
+
# =====================================================================
|
|
48
|
+
# ВСПОМОГАТЕЛЬНЫЕ ФУНКЦИИ
|
|
49
|
+
# =====================================================================
|
|
50
|
+
|
|
51
|
+
def _pressure_exponential_from_10_and_72(
|
|
52
|
+
z_km: float,
|
|
53
|
+
P10_hpa: float,
|
|
54
|
+
P72_hpa: float,
|
|
55
|
+
decay_10_72: float,
|
|
56
|
+
decay_72_100: float,
|
|
57
|
+
) -> float:
|
|
58
|
+
"""
|
|
59
|
+
Универсальная формула для P(Z) в верхних слоях.
|
|
60
|
+
|
|
61
|
+
P.835-7, ур. (10b/10c), (13b/13c), (16b/16c), (19b/19c), (22b/22c).
|
|
62
|
+
|
|
63
|
+
Args:
|
|
64
|
+
z_km: высота.
|
|
65
|
+
P10_hpa: давление на 10 км.
|
|
66
|
+
P72_hpa: давление на 72 км.
|
|
67
|
+
decay_10_72: коэффициент затухания на 10–72 км (обычно 0.147).
|
|
68
|
+
decay_72_100: коэффициент на 72–100 км (обычно 0.150–0.165).
|
|
69
|
+
|
|
70
|
+
Returns:
|
|
71
|
+
Давление, гПа.
|
|
72
|
+
"""
|
|
73
|
+
if z_km <= 10.0:
|
|
74
|
+
raise ValueError("Функция только для Z > 10 км")
|
|
75
|
+
if z_km <= 72.0:
|
|
76
|
+
return P10_hpa * math.exp(-decay_10_72 * (z_km - 10.0))
|
|
77
|
+
return P72_hpa * math.exp(-decay_72_100 * (z_km - 72.0))
|
|
78
|
+
|
|
79
|
+
|
|
80
|
+
# =====================================================================
|
|
81
|
+
# 1. LOW-LATITUDE ANNUAL (15°N, все сезоны)
|
|
82
|
+
# =====================================================================
|
|
83
|
+
|
|
84
|
+
def low_latitude_annual_temperature(z_km: float) -> float:
|
|
85
|
+
"""
|
|
86
|
+
T(Z) для low-latitude annual (15°N).
|
|
87
|
+
P.835-7, ур. (9a–9e).
|
|
88
|
+
"""
|
|
89
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
90
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
91
|
+
|
|
92
|
+
if z_km < 17.0:
|
|
93
|
+
return 300.4222 - 6.3533 * z_km + 0.005886 * z_km * z_km
|
|
94
|
+
elif z_km < 47.0:
|
|
95
|
+
return 194.0 + 2.533 * (z_km - 17.0)
|
|
96
|
+
elif z_km < 52.0:
|
|
97
|
+
return 270.0
|
|
98
|
+
elif z_km < 80.0:
|
|
99
|
+
return 270.0 - 3.0714 * (z_km - 52.0)
|
|
100
|
+
else:
|
|
101
|
+
return 184.0
|
|
102
|
+
|
|
103
|
+
|
|
104
|
+
def low_latitude_annual_pressure(z_km: float) -> float:
|
|
105
|
+
"""
|
|
106
|
+
P(Z) для low-latitude annual (15°N).
|
|
107
|
+
P.835-7, ур. (10a–10c).
|
|
108
|
+
"""
|
|
109
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
110
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
111
|
+
|
|
112
|
+
if z_km <= 10.0:
|
|
113
|
+
return 1012.0306 - 109.0338 * z_km + 3.6316 * z_km * z_km
|
|
114
|
+
P10 = 1012.0306 - 109.0338 * 10.0 + 3.6316 * 100.0
|
|
115
|
+
if z_km <= 72.0:
|
|
116
|
+
return P10 * math.exp(-0.147 * (z_km - 10.0))
|
|
117
|
+
P72 = P10 * math.exp(-0.147 * 62.0)
|
|
118
|
+
return P72 * math.exp(-0.165 * (z_km - 72.0))
|
|
119
|
+
|
|
120
|
+
|
|
121
|
+
def low_latitude_annual_water_vapour(z_km: float) -> float:
|
|
122
|
+
"""
|
|
123
|
+
ρ(Z) для low-latitude annual (15°N).
|
|
124
|
+
P.835-7, ур. (11a, 11b).
|
|
125
|
+
"""
|
|
126
|
+
if z_km < 0.0:
|
|
127
|
+
raise ValueError(f"Z={z_km} < 0")
|
|
128
|
+
if z_km > 15.0:
|
|
129
|
+
return 0.0
|
|
130
|
+
z = z_km
|
|
131
|
+
exponent = (-0.2313 * z - 0.1122 * z * z
|
|
132
|
+
+ 0.01351 * z * z * z
|
|
133
|
+
- 0.0005923 * z * z * z * z)
|
|
134
|
+
return 19.6542 * math.exp(exponent)
|
|
135
|
+
|
|
136
|
+
|
|
137
|
+
# =====================================================================
|
|
138
|
+
# 2. MID-LATITUDE SUMMER (45°N, лето)
|
|
139
|
+
# =====================================================================
|
|
140
|
+
|
|
141
|
+
def mid_latitude_summer_temperature(z_km: float) -> float:
|
|
142
|
+
"""T(Z) для mid-latitude summer. P.835-7, ур. (12a–12f)."""
|
|
143
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
144
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
145
|
+
|
|
146
|
+
if z_km < 13.0:
|
|
147
|
+
return 294.9838 - 5.2159 * z_km - 0.07109 * z_km * z_km
|
|
148
|
+
elif z_km < 17.0:
|
|
149
|
+
return 215.15
|
|
150
|
+
elif z_km < 47.0:
|
|
151
|
+
return 215.15 * math.exp(0.008128 * (z_km - 17.0))
|
|
152
|
+
elif z_km < 53.0:
|
|
153
|
+
return 275.0
|
|
154
|
+
elif z_km < 80.0:
|
|
155
|
+
return 275.0 + 111.57755 * (1.0 - math.exp(0.0237 * (z_km - 53.0)))
|
|
156
|
+
else:
|
|
157
|
+
return 175.0
|
|
158
|
+
|
|
159
|
+
|
|
160
|
+
def mid_latitude_summer_pressure(z_km: float) -> float:
|
|
161
|
+
"""P(Z) для mid-latitude summer. P.835-7, ур. (13a–13c)."""
|
|
162
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
163
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
164
|
+
|
|
165
|
+
if z_km <= 10.0:
|
|
166
|
+
return 1012.8186 - 111.5569 * z_km + 3.8646 * z_km * z_km
|
|
167
|
+
P10 = 1012.8186 - 111.5569 * 10.0 + 3.8646 * 100.0
|
|
168
|
+
if z_km <= 72.0:
|
|
169
|
+
return P10 * math.exp(-0.147 * (z_km - 10.0))
|
|
170
|
+
P72 = P10 * math.exp(-0.147 * 62.0)
|
|
171
|
+
return P72 * math.exp(-0.165 * (z_km - 72.0))
|
|
172
|
+
|
|
173
|
+
|
|
174
|
+
def mid_latitude_summer_water_vapour(z_km: float) -> float:
|
|
175
|
+
"""ρ(Z) для mid-latitude summer. P.835-7, ур. (14a, 14b)."""
|
|
176
|
+
if z_km < 0.0:
|
|
177
|
+
raise ValueError(f"Z={z_km} < 0")
|
|
178
|
+
if z_km > 15.0:
|
|
179
|
+
return 0.0
|
|
180
|
+
z = z_km
|
|
181
|
+
exponent = -0.4174 * z - 0.02290 * z * z + 0.001007 * z * z * z
|
|
182
|
+
return 14.3542 * math.exp(exponent)
|
|
183
|
+
|
|
184
|
+
|
|
185
|
+
# =====================================================================
|
|
186
|
+
# 3. MID-LATITUDE WINTER (45°N, зима)
|
|
187
|
+
# =====================================================================
|
|
188
|
+
|
|
189
|
+
def mid_latitude_winter_temperature(z_km: float) -> float:
|
|
190
|
+
"""T(Z) для mid-latitude winter. P.835-7, ур. (15a–15f)."""
|
|
191
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
192
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
193
|
+
|
|
194
|
+
if z_km < 10.0:
|
|
195
|
+
return 272.7241 - 3.6217 * z_km - 0.1759 * z_km * z_km
|
|
196
|
+
elif z_km < 33.0:
|
|
197
|
+
return 218.0
|
|
198
|
+
elif z_km < 47.0:
|
|
199
|
+
return 218.0 + 3.3571 * (z_km - 33.0)
|
|
200
|
+
elif z_km < 53.0:
|
|
201
|
+
return 265.0
|
|
202
|
+
elif z_km < 80.0:
|
|
203
|
+
return 265.0 - 2.0370 * (z_km - 53.0)
|
|
204
|
+
else:
|
|
205
|
+
return 210.0
|
|
206
|
+
|
|
207
|
+
|
|
208
|
+
def mid_latitude_winter_pressure(z_km: float) -> float:
|
|
209
|
+
"""P(Z) для mid-latitude winter. P.835-7, ур. (16a–16c)."""
|
|
210
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
211
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
212
|
+
|
|
213
|
+
if z_km <= 10.0:
|
|
214
|
+
return 1018.8627 - 124.2954 * z_km + 4.8307 * z_km * z_km
|
|
215
|
+
P10 = 1018.8627 - 124.2954 * 10.0 + 4.8307 * 100.0
|
|
216
|
+
if z_km <= 72.0:
|
|
217
|
+
return P10 * math.exp(-0.147 * (z_km - 10.0))
|
|
218
|
+
P72 = P10 * math.exp(-0.147 * 62.0)
|
|
219
|
+
return P72 * math.exp(-0.155 * (z_km - 72.0))
|
|
220
|
+
|
|
221
|
+
|
|
222
|
+
def mid_latitude_winter_water_vapour(z_km: float) -> float:
|
|
223
|
+
"""ρ(Z) для mid-latitude winter. P.835-7, ур. (17a, 17b)."""
|
|
224
|
+
if z_km < 0.0:
|
|
225
|
+
raise ValueError(f"Z={z_km} < 0")
|
|
226
|
+
if z_km > 10.0:
|
|
227
|
+
return 0.0
|
|
228
|
+
z = z_km
|
|
229
|
+
exponent = -0.2697 * z - 0.03604 * z * z + 0.0004489 * z * z * z
|
|
230
|
+
return 3.4742 * math.exp(exponent)
|
|
231
|
+
|
|
232
|
+
|
|
233
|
+
# =====================================================================
|
|
234
|
+
# 4. HIGH-LATITUDE SUMMER (60°N, лето)
|
|
235
|
+
# =====================================================================
|
|
236
|
+
|
|
237
|
+
def high_latitude_summer_temperature(z_km: float) -> float:
|
|
238
|
+
"""T(Z) для high-latitude summer. P.835-7, ур. (18a–18f)."""
|
|
239
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
240
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
241
|
+
|
|
242
|
+
if z_km < 10.0:
|
|
243
|
+
return 286.8374 - 4.7805 * z_km - 0.1402 * z_km * z_km
|
|
244
|
+
elif z_km < 23.0:
|
|
245
|
+
return 225.0
|
|
246
|
+
elif z_km < 48.0:
|
|
247
|
+
return 225.0 * math.exp(0.008317 * (z_km - 23.0))
|
|
248
|
+
elif z_km < 53.0:
|
|
249
|
+
return 277.0
|
|
250
|
+
elif z_km < 79.0:
|
|
251
|
+
return 277.0 - 4.0769 * (z_km - 53.0)
|
|
252
|
+
else:
|
|
253
|
+
return 171.0
|
|
254
|
+
|
|
255
|
+
|
|
256
|
+
def high_latitude_summer_pressure(z_km: float) -> float:
|
|
257
|
+
"""P(Z) для high-latitude summer. P.835-7, ур. (19a–19c)."""
|
|
258
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
259
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
260
|
+
|
|
261
|
+
if z_km <= 10.0:
|
|
262
|
+
return 1008.0278 - 113.2494 * z_km + 3.9408 * z_km * z_km
|
|
263
|
+
P10 = 1008.0278 - 113.2494 * 10.0 + 3.9408 * 100.0
|
|
264
|
+
if z_km <= 72.0:
|
|
265
|
+
return P10 * math.exp(-0.140 * (z_km - 10.0))
|
|
266
|
+
P72 = P10 * math.exp(-0.140 * 62.0)
|
|
267
|
+
return P72 * math.exp(-0.165 * (z_km - 72.0))
|
|
268
|
+
|
|
269
|
+
|
|
270
|
+
def high_latitude_summer_water_vapour(z_km: float) -> float:
|
|
271
|
+
"""ρ(Z) для high-latitude summer. P.835-7, ур. (20a, 20b)."""
|
|
272
|
+
if z_km < 0.0:
|
|
273
|
+
raise ValueError(f"Z={z_km} < 0")
|
|
274
|
+
if z_km > 15.0:
|
|
275
|
+
return 0.0
|
|
276
|
+
z = z_km
|
|
277
|
+
exponent = -0.3614 * z - 0.005402 * z * z - 0.001955 * z * z * z
|
|
278
|
+
return 8.988 * math.exp(exponent)
|
|
279
|
+
|
|
280
|
+
|
|
281
|
+
# =====================================================================
|
|
282
|
+
# 5. HIGH-LATITUDE WINTER (60°N, зима)
|
|
283
|
+
# =====================================================================
|
|
284
|
+
|
|
285
|
+
def high_latitude_winter_temperature(z_km: float) -> float:
|
|
286
|
+
"""T(Z) для high-latitude winter. P.835-7, ур. (21a–21e)."""
|
|
287
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
288
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
289
|
+
|
|
290
|
+
if z_km < 8.5:
|
|
291
|
+
return (257.4345 + 2.3474 * z_km
|
|
292
|
+
- 1.5479 * z_km * z_km
|
|
293
|
+
+ 0.08473 * z_km * z_km * z_km)
|
|
294
|
+
elif z_km < 30.0:
|
|
295
|
+
return 217.5
|
|
296
|
+
elif z_km < 50.0:
|
|
297
|
+
return 217.5 + 2.125 * (z_km - 30.0)
|
|
298
|
+
elif z_km < 54.0:
|
|
299
|
+
return 260.0
|
|
300
|
+
else:
|
|
301
|
+
return 260.0 - 1.667 * (z_km - 54.0)
|
|
302
|
+
|
|
303
|
+
|
|
304
|
+
def high_latitude_winter_pressure(z_km: float) -> float:
|
|
305
|
+
"""P(Z) для high-latitude winter. P.835-7, ур. (22a–22c)."""
|
|
306
|
+
if z_km < 0.0 or z_km > 100.0:
|
|
307
|
+
raise ValueError(f"Z={z_km} вне [0, 100]")
|
|
308
|
+
|
|
309
|
+
if z_km <= 10.0:
|
|
310
|
+
return 1010.8828 - 122.2411 * z_km + 4.554 * z_km * z_km
|
|
311
|
+
P10 = 1010.8828 - 122.2411 * 10.0 + 4.554 * 100.0
|
|
312
|
+
if z_km <= 72.0:
|
|
313
|
+
return P10 * math.exp(-0.147 * (z_km - 10.0))
|
|
314
|
+
P72 = P10 * math.exp(-0.147 * 62.0)
|
|
315
|
+
return P72 * math.exp(-0.150 * (z_km - 72.0))
|
|
316
|
+
|
|
317
|
+
|
|
318
|
+
def high_latitude_winter_water_vapour(z_km: float) -> float:
|
|
319
|
+
"""ρ(Z) для high-latitude winter. P.835-7, ур. (23a, 23b)."""
|
|
320
|
+
if z_km < 0.0:
|
|
321
|
+
raise ValueError(f"Z={z_km} < 0")
|
|
322
|
+
if z_km > 10.0:
|
|
323
|
+
return 0.0
|
|
324
|
+
z = z_km
|
|
325
|
+
exponent = 0.07481 * z - 0.0981 * z * z + 0.00281 * z * z * z
|
|
326
|
+
return 1.2319 * math.exp(exponent)
|
|
327
|
+
|
|
328
|
+
|
|
329
|
+
# =====================================================================
|
|
330
|
+
# УНИВЕРСАЛЬНЫЙ ДОСТУП ПО РЕГИОНУ
|
|
331
|
+
# =====================================================================
|
|
332
|
+
|
|
333
|
+
_REGION_TEMPERATURE = {
|
|
334
|
+
AtmosphereRegion.LOW_LATITUDE_ANNUAL: low_latitude_annual_temperature,
|
|
335
|
+
AtmosphereRegion.MID_LATITUDE_SUMMER: mid_latitude_summer_temperature,
|
|
336
|
+
AtmosphereRegion.MID_LATITUDE_WINTER: mid_latitude_winter_temperature,
|
|
337
|
+
AtmosphereRegion.HIGH_LATITUDE_SUMMER: high_latitude_summer_temperature,
|
|
338
|
+
AtmosphereRegion.HIGH_LATITUDE_WINTER: high_latitude_winter_temperature,
|
|
339
|
+
}
|
|
340
|
+
|
|
341
|
+
_REGION_PRESSURE = {
|
|
342
|
+
AtmosphereRegion.LOW_LATITUDE_ANNUAL: low_latitude_annual_pressure,
|
|
343
|
+
AtmosphereRegion.MID_LATITUDE_SUMMER: mid_latitude_summer_pressure,
|
|
344
|
+
AtmosphereRegion.MID_LATITUDE_WINTER: mid_latitude_winter_pressure,
|
|
345
|
+
AtmosphereRegion.HIGH_LATITUDE_SUMMER: high_latitude_summer_pressure,
|
|
346
|
+
AtmosphereRegion.HIGH_LATITUDE_WINTER: high_latitude_winter_pressure,
|
|
347
|
+
}
|
|
348
|
+
|
|
349
|
+
_REGION_WATER_VAPOUR = {
|
|
350
|
+
AtmosphereRegion.LOW_LATITUDE_ANNUAL: low_latitude_annual_water_vapour,
|
|
351
|
+
AtmosphereRegion.MID_LATITUDE_SUMMER: mid_latitude_summer_water_vapour,
|
|
352
|
+
AtmosphereRegion.MID_LATITUDE_WINTER: mid_latitude_winter_water_vapour,
|
|
353
|
+
AtmosphereRegion.HIGH_LATITUDE_SUMMER: high_latitude_summer_water_vapour,
|
|
354
|
+
AtmosphereRegion.HIGH_LATITUDE_WINTER: high_latitude_winter_water_vapour,
|
|
355
|
+
}
|
|
356
|
+
|
|
357
|
+
|
|
358
|
+
def regional_temperature(z_km: float, region: AtmosphereRegion) -> float:
|
|
359
|
+
"""T(Z) для заданного региона."""
|
|
360
|
+
return _REGION_TEMPERATURE[region](z_km)
|
|
361
|
+
|
|
362
|
+
|
|
363
|
+
def regional_pressure(z_km: float, region: AtmosphereRegion) -> float:
|
|
364
|
+
"""P(Z) для заданного региона."""
|
|
365
|
+
return _REGION_PRESSURE[region](z_km)
|
|
366
|
+
|
|
367
|
+
|
|
368
|
+
def regional_water_vapour(z_km: float, region: AtmosphereRegion) -> float:
|
|
369
|
+
"""ρ(Z) для заданного региона."""
|
|
370
|
+
return _REGION_WATER_VAPOUR[region](z_km)
|
|
371
|
+
|
|
372
|
+
|
|
373
|
+
def get_region_for_latitude(
|
|
374
|
+
lat_deg: float, season: Season
|
|
375
|
+
) -> AtmosphereRegion:
|
|
376
|
+
"""
|
|
377
|
+
Возвращает региональную атмосферу по широте и сезону.
|
|
378
|
+
|
|
379
|
+
Логика Annex 2:
|
|
380
|
+
0°–15°N: low annual (все сезоны)
|
|
381
|
+
≥ 60°N: high (по сезону)
|
|
382
|
+
между 15° и 45° или 45° и 60° — интерполяция (см. interp)
|
|
383
|
+
Для отрицательных широт (южное полушарие) — по модулю.
|
|
384
|
+
|
|
385
|
+
Args:
|
|
386
|
+
lat_deg: широта в градусах.
|
|
387
|
+
season: Season.SUMMER или Season.WINTER.
|
|
388
|
+
|
|
389
|
+
Returns:
|
|
390
|
+
AtmosphereRegion (без интерполяции). Для интерполяции
|
|
391
|
+
используйте interpolated_temperature / _pressure / _water_vapour.
|
|
392
|
+
"""
|
|
393
|
+
abs_lat = abs(lat_deg)
|
|
394
|
+
if abs_lat <= 15.0:
|
|
395
|
+
return AtmosphereRegion.LOW_LATITUDE_ANNUAL
|
|
396
|
+
if abs_lat >= 60.0:
|
|
397
|
+
return (AtmosphereRegion.HIGH_LATITUDE_SUMMER if season == Season.SUMMER
|
|
398
|
+
else AtmosphereRegion.HIGH_LATITUDE_WINTER)
|
|
399
|
+
if abs_lat <= 45.0:
|
|
400
|
+
return (AtmosphereRegion.MID_LATITUDE_SUMMER if season == Season.SUMMER
|
|
401
|
+
else AtmosphereRegion.MID_LATITUDE_WINTER)
|
|
402
|
+
return (AtmosphereRegion.HIGH_LATITUDE_SUMMER if season == Season.SUMMER
|
|
403
|
+
else AtmosphereRegion.HIGH_LATITUDE_WINTER)
|
|
404
|
+
|
|
405
|
+
|
|
406
|
+
def _interp_between_regions(
|
|
407
|
+
z_km: float,
|
|
408
|
+
lat_deg: float,
|
|
409
|
+
season: Season,
|
|
410
|
+
accessor,
|
|
411
|
+
) -> float:
|
|
412
|
+
"""
|
|
413
|
+
Внутренняя: интерполяция между двумя регионами по широте.
|
|
414
|
+
|
|
415
|
+
Args:
|
|
416
|
+
z_km: высота.
|
|
417
|
+
lat_deg: широта.
|
|
418
|
+
season: сезон.
|
|
419
|
+
accessor: функция, возвращающая T(Z, region) для нужного параметра.
|
|
420
|
+
"""
|
|
421
|
+
abs_lat = abs(lat_deg)
|
|
422
|
+
|
|
423
|
+
if abs_lat <= 15.0:
|
|
424
|
+
return accessor(z_km, AtmosphereRegion.LOW_LATITUDE_ANNUAL)
|
|
425
|
+
|
|
426
|
+
low_val = accessor(z_km, AtmosphereRegion.LOW_LATITUDE_ANNUAL)
|
|
427
|
+
|
|
428
|
+
if abs_lat <= 45.0:
|
|
429
|
+
mid_region = (AtmosphereRegion.MID_LATITUDE_SUMMER if season == Season.SUMMER
|
|
430
|
+
else AtmosphereRegion.MID_LATITUDE_WINTER)
|
|
431
|
+
mid_val = accessor(z_km, mid_region)
|
|
432
|
+
t = (abs_lat - 15.0) / (45.0 - 15.0)
|
|
433
|
+
return low_val * (1.0 - t) + mid_val * t
|
|
434
|
+
|
|
435
|
+
if abs_lat < 60.0:
|
|
436
|
+
mid_region = (AtmosphereRegion.MID_LATITUDE_SUMMER if season == Season.SUMMER
|
|
437
|
+
else AtmosphereRegion.MID_LATITUDE_WINTER)
|
|
438
|
+
high_region = (AtmosphereRegion.HIGH_LATITUDE_SUMMER if season == Season.SUMMER
|
|
439
|
+
else AtmosphereRegion.HIGH_LATITUDE_WINTER)
|
|
440
|
+
mid_val = accessor(z_km, mid_region)
|
|
441
|
+
high_val = accessor(z_km, high_region)
|
|
442
|
+
t = (abs_lat - 45.0) / (60.0 - 45.0)
|
|
443
|
+
return mid_val * (1.0 - t) + high_val * t
|
|
444
|
+
|
|
445
|
+
return accessor(z_km, (
|
|
446
|
+
AtmosphereRegion.HIGH_LATITUDE_SUMMER if season == Season.SUMMER
|
|
447
|
+
else AtmosphereRegion.HIGH_LATITUDE_WINTER
|
|
448
|
+
))
|
|
449
|
+
|
|
450
|
+
|
|
451
|
+
def interpolated_temperature(
|
|
452
|
+
z_km: float, lat_deg: float, season: Season
|
|
453
|
+
) -> float:
|
|
454
|
+
"""T(Z) с линейной интерполяцией по широте."""
|
|
455
|
+
return _interp_between_regions(z_km, lat_deg, season, regional_temperature)
|
|
456
|
+
|
|
457
|
+
|
|
458
|
+
def interpolated_pressure(
|
|
459
|
+
z_km: float, lat_deg: float, season: Season
|
|
460
|
+
) -> float:
|
|
461
|
+
"""P(Z) с линейной интерполяцией по широте."""
|
|
462
|
+
return _interp_between_regions(z_km, lat_deg, season, regional_pressure)
|
|
463
|
+
|
|
464
|
+
|
|
465
|
+
def interpolated_water_vapour(
|
|
466
|
+
z_km: float, lat_deg: float, season: Season
|
|
467
|
+
) -> float:
|
|
468
|
+
"""ρ(Z) с линейной интерполяцией по широте."""
|
|
469
|
+
return _interp_between_regions(z_km, lat_deg, season, regional_water_vapour)
|
|
470
|
+
|
|
471
|
+
|
|
472
|
+
__all__ = [
|
|
473
|
+
'AtmosphereRegion', 'Season',
|
|
474
|
+
# Low
|
|
475
|
+
'low_latitude_annual_temperature',
|
|
476
|
+
'low_latitude_annual_pressure',
|
|
477
|
+
'low_latitude_annual_water_vapour',
|
|
478
|
+
# Mid summer
|
|
479
|
+
'mid_latitude_summer_temperature',
|
|
480
|
+
'mid_latitude_summer_pressure',
|
|
481
|
+
'mid_latitude_summer_water_vapour',
|
|
482
|
+
# Mid winter
|
|
483
|
+
'mid_latitude_winter_temperature',
|
|
484
|
+
'mid_latitude_winter_pressure',
|
|
485
|
+
'mid_latitude_winter_water_vapour',
|
|
486
|
+
# High summer
|
|
487
|
+
'high_latitude_summer_temperature',
|
|
488
|
+
'high_latitude_summer_pressure',
|
|
489
|
+
'high_latitude_summer_water_vapour',
|
|
490
|
+
# High winter
|
|
491
|
+
'high_latitude_winter_temperature',
|
|
492
|
+
'high_latitude_winter_pressure',
|
|
493
|
+
'high_latitude_winter_water_vapour',
|
|
494
|
+
# Региональный доступ
|
|
495
|
+
'regional_temperature', 'regional_pressure', 'regional_water_vapour',
|
|
496
|
+
'get_region_for_latitude',
|
|
497
|
+
'interpolated_temperature', 'interpolated_pressure',
|
|
498
|
+
'interpolated_water_vapour',
|
|
499
|
+
]
|
|
@@ -0,0 +1,70 @@
|
|
|
1
|
+
Metadata-Version: 2.4
|
|
2
|
+
Name: p835
|
|
3
|
+
Version: 0.1.0
|
|
4
|
+
Summary: ITU-R P.835-7 reference standard atmospheres
|
|
5
|
+
Author-email: Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
6
|
+
License: MIT
|
|
7
|
+
Classifier: Programming Language :: Python :: 3
|
|
8
|
+
Classifier: Operating System :: OS Independent
|
|
9
|
+
Requires-Python: >=3.9
|
|
10
|
+
Description-Content-Type: text/markdown
|
|
11
|
+
License-File: LICENSE
|
|
12
|
+
License-File: NOTICE
|
|
13
|
+
Provides-Extra: test
|
|
14
|
+
Requires-Dist: pytest>=7; extra == "test"
|
|
15
|
+
Dynamic: license-file
|
|
16
|
+
|
|
17
|
+
# p835
|
|
18
|
+
|
|
19
|
+
Python-реализация **ITU-R P.835-7** — «Reference standard
|
|
20
|
+
atmospheres».
|
|
21
|
+
|
|
22
|
+
## Установка
|
|
23
|
+
|
|
24
|
+
```bash
|
|
25
|
+
pip install -e .
|
|
26
|
+
```
|
|
27
|
+
|
|
28
|
+
## Использование
|
|
29
|
+
|
|
30
|
+
```python
|
|
31
|
+
from p835 import global_temperature, global_pressure
|
|
32
|
+
|
|
33
|
+
T = global_temperature(10.0) # K
|
|
34
|
+
P = global_pressure(10.0) # гПа
|
|
35
|
+
```
|
|
36
|
+
|
|
37
|
+
## Тесты
|
|
38
|
+
|
|
39
|
+
```bash
|
|
40
|
+
pytest tests/ -v
|
|
41
|
+
```
|
|
42
|
+
|
|
43
|
+
## Attribution
|
|
44
|
+
|
|
45
|
+
ATTRIBUTION AND SOURCES
|
|
46
|
+
========================
|
|
47
|
+
|
|
48
|
+
This software (PyP835) is a Python implementation of the mathematical
|
|
49
|
+
model described in:
|
|
50
|
+
|
|
51
|
+
ITU-R P.835-7 — "Reference standard atmospheres"
|
|
52
|
+
International Telecommunication Union, Radiocommunication Sector (ITU-R)
|
|
53
|
+
Geneva, Switzerland
|
|
54
|
+
|
|
55
|
+
The ITU-R Recommendations are publicly available at:
|
|
56
|
+
https://www.itu.int/rec/R-REC-P.835-7/en
|
|
57
|
+
|
|
58
|
+
The original reference implementation in C++ was published by the
|
|
59
|
+
U.S. National Telecommunications and Information Administration (NTIA)
|
|
60
|
+
and is available at:
|
|
61
|
+
https://github.com/NTIA/p835-7
|
|
62
|
+
|
|
63
|
+
The Python port was created by:
|
|
64
|
+
Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>, 2026
|
|
65
|
+
|
|
66
|
+
This software is provided under the MIT License (see LICENSE).
|
|
67
|
+
|
|
68
|
+
The author of this port is not affiliated with ITU or NTIA. All
|
|
69
|
+
formulas and algorithms are taken from the publicly available
|
|
70
|
+
ITU-R Recommendations.
|
|
@@ -0,0 +1,13 @@
|
|
|
1
|
+
p835/__init__.py,sha256=2YEnoyrUcWaG5uSJXekUECijfxgQvstbGXGGfiuyFtc,2697
|
|
2
|
+
p835/conversions.py,sha256=ESW7U_N6joZONEx5U8_zHW7IXjVWuf1UzE3Nv8blQII,2935
|
|
3
|
+
p835/reference_atmosphere.py,sha256=8wuEy4FEGzBDP9t_EZ-vo_8fKPP6BwvjdXuEX3kTL3w,10243
|
|
4
|
+
p835/seasonal_atmosphere.py,sha256=3Ul5eECIse_H5DoJy0ZvZu-bKJ7ignH4zq8xdzNmRcY,17571
|
|
5
|
+
p835/common/__init__.py,sha256=AY7RrZJI-hHT_X8mXldyJ2Wz9fuvwPaVLSB-jbub9LQ,194
|
|
6
|
+
p835/common/constants.py,sha256=7lGZRkKlLaeXN5GtOCPsUqYgdP0QMsPEt89pCfthR7c,5977
|
|
7
|
+
p835/common/exceptions.py,sha256=pRV0tWU_eirLy_DlqNGGWXoEe0PlLx0fgu6sUKeWd1o,1308
|
|
8
|
+
p835-0.1.0.dist-info/licenses/LICENSE,sha256=OnEyEitlqkdmtth_6FPN5BuLyQTqPrVwKXuiyjf2QZo,1108
|
|
9
|
+
p835-0.1.0.dist-info/licenses/NOTICE,sha256=5nq1y470n6BsCo9QAKqPYrcSPPpxt680mX3Tz3ZvEEE,899
|
|
10
|
+
p835-0.1.0.dist-info/METADATA,sha256=YCAIRUycPZcnKSKTFanpjjSDdvufWYZpMjEepmTL2Gg,1751
|
|
11
|
+
p835-0.1.0.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
|
|
12
|
+
p835-0.1.0.dist-info/top_level.txt,sha256=eMb386RSuQ1JZ3spKQDPtboxJQXfXTl-fBzKj6PjM2I,5
|
|
13
|
+
p835-0.1.0.dist-info/RECORD,,
|
|
@@ -0,0 +1,21 @@
|
|
|
1
|
+
MIT License
|
|
2
|
+
|
|
3
|
+
Copyright (c) 2026 Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>
|
|
4
|
+
|
|
5
|
+
Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
6
|
+
of this software and associated documentation files (the "Software"), to deal
|
|
7
|
+
in the Software without restriction, including without limitation the rights
|
|
8
|
+
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
9
|
+
copies of the Software, and to permit persons to whom the Software is
|
|
10
|
+
furnished to do so, subject to the following conditions:
|
|
11
|
+
|
|
12
|
+
The above copyright notice and this permission notice shall be included in all
|
|
13
|
+
copies or substantial portions of the Software.
|
|
14
|
+
|
|
15
|
+
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
16
|
+
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
17
|
+
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
18
|
+
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
19
|
+
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
20
|
+
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
|
|
21
|
+
SOFTWARE.
|
|
@@ -0,0 +1,26 @@
|
|
|
1
|
+
ATTRIBUTION AND SOURCES
|
|
2
|
+
========================
|
|
3
|
+
|
|
4
|
+
This software (PyP835) is a Python implementation of the mathematical
|
|
5
|
+
model described in:
|
|
6
|
+
|
|
7
|
+
ITU-R P.835-7 — "Reference standard atmospheres"
|
|
8
|
+
International Telecommunication Union, Radiocommunication Sector (ITU-R)
|
|
9
|
+
Geneva, Switzerland
|
|
10
|
+
|
|
11
|
+
The ITU-R Recommendations are publicly available at:
|
|
12
|
+
https://www.itu.int/rec/R-REC-P.835-7/en
|
|
13
|
+
|
|
14
|
+
The original reference implementation in C++ was published by the
|
|
15
|
+
U.S. National Telecommunications and Information Administration (NTIA)
|
|
16
|
+
and is available at:
|
|
17
|
+
https://github.com/NTIA/p835-7
|
|
18
|
+
|
|
19
|
+
The Python port was created by:
|
|
20
|
+
Evgeny Gennadievich Prutkov <sviazist.pr@gmail.com>, 2026
|
|
21
|
+
|
|
22
|
+
This software is provided under the MIT License (see LICENSE).
|
|
23
|
+
|
|
24
|
+
The author of this port is not affiliated with ITU or NTIA. All
|
|
25
|
+
formulas and algorithms are taken from the publicly available
|
|
26
|
+
ITU-R Recommendations.
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
p835
|