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 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
+ ]
@@ -0,0 +1,7 @@
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
+
@@ -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,5 @@
1
+ Wheel-Version: 1.0
2
+ Generator: setuptools (84.0.0)
3
+ Root-Is-Purelib: true
4
+ Tag: py3-none-any
5
+
@@ -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