Source code for sasktran2.climatology.atomic_oxygen
from __future__ import annotations
import calendar
import datetime as _dt
from typing import Any
import numpy as np
import numpy.typing as npt
MONTHS = (
"January",
"February",
"March",
"April",
"May",
"June",
"July",
"August",
"September",
"October",
"November",
"December",
)
ALTITUDES_M = np.array(
[
130000.0,
125000.0,
120000.0,
115000.0,
110000.0,
105000.0,
100000.0,
95000.0,
90000.0,
85000.0,
80000.0,
75000.0,
70000.0,
65000.0,
60000.0,
55000.0,
50000.0,
45000.0,
40000.0,
],
dtype=float,
)
LATITUDES_DEGREES = np.array(
[-80, -70, -60, -50, -40, -30, -20, -10, 0, 10, 20, 30, 40, 50, 60, 70, 80],
dtype=float,
)
# Atomic oxygen number density in cm^-3. The table is ordered as
# [month, altitude, latitude], using MONTHS, ALTITUDES_M, and LATITUDES_DEGREES.
# Source: tabulated zonal mean reference data from NASA NTRS report 19910018329,
# Tables 1 and 2, originally published as Adv. Space Res. 17(11), 1996.
_TABLE_CM3 = np.array(
[
[
[
2.500000e10,
2.800000e10,
3.100000e10,
3.600000e10,
3.600000e10,
3.700000e10,
3.700000e10,
3.700000e10,
3.200000e10,
3.100000e10,
3.100000e10,
3.100000e10,
3.000000e10,
2.900000e10,
2.700000e10,
2.600000e10,
2.400000e10,
],
[
3.600000e10,
3.900000e10,
4.400000e10,
5.100000e10,
5.100000e10,
5.100000e10,
5.100000e10,
5.000000e10,
4.300000e10,
4.300000e10,
4.300000e10,
4.300000e10,
4.200000e10,
4.100000e10,
3.900000e10,
3.600000e10,
3.400000e10,
],
[
5.400000e10,
5.900000e10,
6.600000e10,
7.600000e10,
7.600000e10,
7.700000e10,
7.500000e10,
7.300000e10,
6.300000e10,
6.200000e10,
6.300000e10,
6.400000e10,
6.400000e10,
6.200000e10,
5.900000e10,
5.500000e10,
5.200000e10,
],
[
8.900000e10,
9.800000e10,
1.100000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
1.000000e11,
1.000000e11,
1.000000e11,
1.100000e11,
1.100000e11,
1.000000e11,
9.700000e10,
9.000000e10,
8.500000e10,
],
[
1.200000e11,
1.400000e11,
1.600000e11,
1.000000e11,
1.000000e11,
2.000000e11,
2.100000e11,
2.200000e11,
1.900000e11,
1.900000e11,
1.800000e11,
1.700000e11,
1.600000e11,
1.500000e11,
1.400000e11,
1.300000e11,
1.200000e11,
],
[
1.700000e11,
1.900000e11,
2.200000e11,
2.800000e11,
2.800000e11,
3.000000e11,
3.200000e11,
3.400000e11,
2.900000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.300000e11,
2.100000e11,
1.000000e11,
1.700000e11,
1.600000e11,
],
[
2.500000e11,
2.600000e11,
2.900000e11,
3.700000e11,
3.700000e11,
4.000000e11,
4.200000e11,
4.200000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.400000e11,
3.100000e11,
2.800000e11,
2.500000e11,
2.400000e11,
2.300000e11,
],
[
3.000000e11,
3.000000e11,
3.200000e11,
3.800000e11,
3.500000e11,
4.000000e11,
4.000000e11,
3.900000e11,
3.300000e11,
3.300000e11,
3.400000e11,
5.400000e11,
3.200000e11,
3.000000e11,
2.800000e11,
2.700000e11,
2.800000e11,
],
[
2.400000e11,
2.300000e11,
2.200000e11,
2.400000e11,
2.400000e11,
2.300000e11,
2.200000e11,
2.100000e11,
1.000000e11,
1.500000e11,
1.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.300000e11,
],
[
9.100000e10,
8.100000e10,
7.300000e10,
6.500000e10,
6.500000e10,
5.800000e10,
5.200000e10,
4.600000e10,
3.900000e10,
3.900000e10,
4.400000e10,
4.900000e10,
5.500000e10,
6.000000e10,
6.500000e10,
7.400000e10,
8.600000e10,
],
[
4.300000e09,
4.500000e09,
5.000000e09,
7.500000e09,
1.100000e10,
1.500000e10,
1.900000e10,
2.100000e10,
3.800000e10,
2.300000e10,
2.400000e10,
2.400000e10,
2.400000e10,
2.900000e10,
2.500000e10,
0.000000e00,
0.000000e00,
],
[
4.900000e09,
4.400000e09,
4.100000e09,
4.200000e09,
4.300000e09,
4.600000e09,
5.200000e09,
5.200000e09,
5.200000e09,
5.700000e09,
7.100000e09,
9.700000e09,
1.300000e10,
1.500000e10,
1.400000e10,
0.000000e00,
0.000000e00,
],
[
6.300000e09,
5.900000e09,
5.600000e09,
5.200000e09,
4.600000e09,
4.100000e09,
4.200000e09,
4.500000e09,
4.700000e09,
4.500000e09,
4.700000e09,
5.400000e09,
6.300000e09,
9.400000e09,
1.200000e10,
0.000000e00,
0.000000e00,
],
[
7.200000e09,
7.000000e09,
6.800000e09,
6.200000e09,
5.800000e09,
5.700000e09,
5.900000e09,
6.100000e09,
6.000000e09,
5.800000e09,
5.800000e09,
6.000000e09,
7.600000e09,
7.600000e09,
4.900000e09,
0.000000e00,
0.000000e00,
],
[
6.400000e09,
6.600000e09,
6.800000e09,
6.800000e09,
6.800000e09,
6.900000e09,
7.200000e09,
7.300000e09,
7.200000e09,
6.900000e09,
6.600000e09,
6.400000e09,
6.000000e09,
7.200000e09,
6.400000e09,
0.000000e00,
0.000000e00,
],
[
4.800000e09,
5.200000e09,
5.600000e09,
5.900000e09,
6.100000e09,
6.300000e09,
6.500000e09,
6.700000e09,
6.700000e09,
6.500000e09,
6.200000e09,
5.800000e09,
5.500000e09,
5.400000e09,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
3.000000e09,
2.800000e09,
3.700000e09,
4.000000e09,
4.300000e09,
4.600000e09,
4.700000e09,
4.800000e09,
4.800000e09,
4.700000e09,
4.400000e09,
4.100000e09,
3.700000e09,
2.800000e09,
1.900000e09,
0.000000e00,
0.000000e00,
],
[
1.200000e09,
1.400000e09,
1.600000e09,
1.800000e09,
2.000000e09,
2.200000e09,
2.300000e09,
2.300000e09,
2.300000e09,
2.100000e09,
2.000000e09,
1.700000e09,
1.300000e09,
8.900000e08,
4.500000e08,
0.000000e00,
0.000000e00,
],
[
2.900000e08,
3.900000e08,
4.600000e08,
5.200000e08,
5.800000e08,
6.100000e08,
6.300000e08,
6.400000e08,
6.200000e08,
5.600000e08,
4.900000e08,
4.000000e08,
2.900000e08,
2.000000e08,
1.100000e08,
0.000000e00,
0.000000e00,
],
],
[
[
2.700000e10,
3.000000e10,
3.300000e10,
3.600000e10,
3.800000e10,
3.900000e10,
3.800000e10,
3.700000e10,
3.400000e10,
3.300000e10,
3.300000e10,
3.300000e10,
3.300000e10,
3.100000e10,
3.000000e10,
2.800000e10,
2.800000e10,
],
[
3.800000e10,
4.200000e10,
4.600000e10,
5.100000e10,
5.300000e10,
5.300000e10,
5.200000e10,
5.100000e10,
4.600000e10,
4.500000e10,
4.475000e10,
4.450000e10,
4.425000e10,
4.400000e10,
4.150000e10,
3.900000e10,
3.700000e10,
],
[
5.800000e10,
6.300000e10,
6.900000e10,
7.500000e10,
7.800000e10,
7.800000e10,
7.600000e10,
7.400000e10,
6.600000e10,
6.600000e10,
6.600000e10,
6.700000e10,
6.700000e10,
6.500000e10,
6.200000e10,
4.800000e10,
5.600000e10,
],
[
9.300000e10,
1.000000e11,
1.100000e11,
1.200000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.000000e11,
9.400000e10,
4.000000e10,
],
[
1.300000e11,
1.400000e11,
1.600000e11,
1.500000e11,
1.000000e11,
2.000000e11,
2.100000e11,
2.200000e11,
2.000000e11,
1.900000e11,
1.900000e11,
1.500000e11,
1.700000e11,
1.600000e11,
1.500000e11,
1.300000e11,
1.300000e11,
],
[
1.000000e11,
2.000000e11,
2.300000e11,
2.500000e11,
2.500000e11,
3.000000e11,
3.200000e11,
3.300000e11,
3.000000e11,
3.000000e11,
2.800000e11,
2.600000e11,
2.400000e11,
2.200000e11,
2.000000e11,
1.900000e11,
1.800000e11,
],
[
2.600000e11,
2.700000e11,
3.000000e11,
3.300000e11,
3.700000e11,
3.000000e11,
4.100000e11,
4.000000e11,
3.700000e11,
3.600000e11,
3.600000e11,
5.500000e11,
3.200000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.500000e11,
],
[
3.000000e11,
3.100000e11,
3.200000e11,
3.400000e11,
3.700000e11,
3.500000e11,
3.500000e11,
3.700000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.200000e11,
3.000000e11,
2.900000e11,
2.900000e11,
2.900000e11,
],
[
2.200000e11,
2.200000e11,
2.100000e11,
2.200000e11,
2.200000e11,
2.200000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.800000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.200000e11,
],
[
7.200000e10,
6.800000e10,
6.400000e10,
6.200000e10,
5.800000e10,
5.300000e10,
4.800000e10,
4.300000e10,
3.900000e10,
3.900000e10,
4.300000e10,
4.700000e10,
5.200000e10,
5.500000e10,
5.900000e10,
6.500000e10,
7.100000e10,
],
[
5.000000e09,
5.200000e09,
6.400000e09,
8.900000e09,
1.300000e10,
1.800000e10,
2.200000e10,
2.200000e10,
2.300000e10,
2.800000e10,
2.800000e10,
2.800000e10,
2.700000e10,
2.600000e10,
2.400000e10,
1.700000e10,
0.000000e00,
],
[
5.000000e09,
3.500000e09,
4.600000e09,
4.600000e09,
4.500000e09,
5.200000e09,
5.200000e09,
5.400000e09,
5.700000e09,
6.000000e09,
7.300000e09,
9.500000e09,
1.400000e10,
1.800000e10,
1.800000e10,
1.500000e10,
0.000000e00,
],
[
6.100000e09,
4.600000e09,
5.800000e09,
5.300000e09,
4.900000e09,
4.500000e09,
4.600000e09,
4.600000e09,
4.500000e09,
4.600000e09,
5.000000e09,
5.500000e09,
6.200000e09,
7.800000e09,
1.100000e10,
1.200000e10,
0.000000e00,
],
[
7.400000e09,
5.600000e09,
6.600000e09,
6.200000e09,
5.900000e09,
5.900000e09,
6.100000e09,
6.000000e09,
5.900000e09,
5.500000e09,
5.500000e09,
5.700000e09,
6.100000e09,
7.600000e09,
8.200000e09,
6.500000e09,
0.000000e00,
],
[
6.900000e09,
5.800000e09,
6.700000e09,
6.600000e09,
6.800000e09,
6.900000e09,
7.000000e09,
7.000000e09,
6.800000e09,
6.600000e09,
6.200000e09,
6.000000e09,
5.700000e09,
5.800000e09,
6.200000e09,
5.600000e09,
0.000000e00,
],
[
5.200000e09,
4.800000e09,
5.700000e09,
5.800000e09,
6.000000e09,
6.200000e09,
6.200000e09,
6.500000e09,
6.500000e09,
6.500000e09,
6.100000e09,
5.800000e09,
5.600000e09,
5.400000e09,
5.200000e09,
2.700000e09,
0.000000e00,
],
[
3.200000e09,
3.300000e09,
3.800000e09,
4.100000e09,
4.400000e09,
4.600000e09,
4.700000e09,
4.800000e09,
4.800000e09,
4.700000e09,
4.500000e09,
4.300000e09,
4.100000e09,
3.700000e09,
2.800000e09,
7.600000e08,
0.000000e00,
],
[
1.100000e09,
1.400000e09,
1.600000e09,
1.800000e09,
2.000000e09,
2.200000e09,
2.300000e09,
2.400000e09,
2.400000e09,
2.200000e09,
2.100000e09,
1.900000e09,
1.600000e09,
1.200000e09,
7.600000e08,
1.900000e08,
0.000000e00,
],
[
2.700000e08,
3.600000e08,
4.300000e08,
4.800000e08,
5.400000e08,
5.900000e08,
6.400000e08,
6.800000e08,
6.700000e08,
6.200000e08,
5.500000e08,
4.600000e08,
3.700000e08,
2.700000e08,
1.600000e08,
7.200000e07,
0.000000e00,
],
],
[
[
3.100000e10,
3.300000e10,
3.700000e10,
4.000000e10,
4.100000e10,
4.100000e10,
4.000000e10,
3.900000e10,
3.800000e10,
3.900000e10,
3.900000e10,
4.000000e10,
3.900000e10,
3.800000e10,
3.700000e10,
3.500000e10,
3.500000e10,
],
[
4.300000e10,
4.600000e10,
5.100000e10,
5.500000e10,
5.700000e10,
5.600000e10,
5.500000e10,
5.300000e10,
5.200000e10,
5.200000e10,
5.300000e10,
5.400000e10,
5.400000e10,
5.300000e10,
5.100000e10,
4.900000e10,
4.700000e10,
],
[
6.300000e10,
6.800000e10,
7.400000e10,
8.000000e10,
8.200000e10,
8.200000e10,
7.900000e10,
7.600000e10,
7.400000e10,
7.500000e10,
7.600000e10,
7.800000e10,
7.800000e10,
7.700000e10,
7.400000e10,
7.100000e10,
6.900000e10,
],
[
1.000000e11,
1.100000e11,
1.200000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.150000e11,
1.100000e11,
],
[
1.400000e11,
1.600000e11,
1.700000e11,
1.900000e11,
2.000000e11,
2.100000e11,
2.100000e11,
2.200000e11,
2.200000e11,
2.200000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.500000e11,
1.700000e11,
1.600000e11,
1.600000e11,
],
[
2.100000e11,
2.200000e11,
2.400000e11,
2.700000e11,
2.900000e11,
3.100000e11,
3.200000e11,
3.300000e11,
3.300000e11,
3.300000e11,
3.200000e11,
3.000000e11,
2.800000e11,
2.600000e11,
2.500000e11,
2.300000e11,
2.300000e11,
],
[
2.800000e11,
2.900000e11,
3.200000e11,
3.500000e11,
3.800000e11,
4.000000e11,
4.000000e11,
4.000000e11,
4.000000e11,
4.000000e11,
3.900000e11,
3.800000e11,
3.600000e11,
3.400000e11,
3.200000e11,
3.100000e11,
3.000000e11,
],
[
3.000000e11,
3.100000e11,
3.200000e11,
3.500000e11,
3.700000e11,
3.700000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.300000e11,
3.200000e11,
3.100000e11,
3.100000e11,
],
[
1.000000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.000000e11,
2.000000e11,
2.000000e11,
1.900000e11,
1.000000e11,
2.000000e11,
2.000000e11,
1.900000e11,
1.000000e11,
1.900000e11,
1.500000e11,
],
[
5.200000e10,
5.300000e10,
5.300000e10,
5.300000e10,
5.100000e10,
4.900000e10,
4.600000e10,
4.200000e10,
4.000000e10,
4.200000e10,
4.500000e10,
4.600000e10,
4.700000e10,
4.700000e10,
4.700000e10,
4.600000e10,
4.600000e10,
],
[
9.800000e09,
1.500000e10,
1.600000e10,
1.900000e10,
2.300000e10,
2.500000e10,
2.300000e10,
2.100000e10,
2.000000e10,
2.200000e10,
2.600000e10,
2.900000e10,
3.200000e10,
3.100000e10,
2.800000e10,
2.300000e10,
1.500000e10,
],
[
6.900000e09,
6.200000e09,
5.800000e09,
5.700000e09,
5.900000e09,
5.800000e09,
5.600000e09,
5.400000e09,
5.800000e09,
5.900000e09,
6.700000e09,
8.600000e09,
1.200000e10,
1.700000e10,
2.100000e10,
2.000000e10,
1.900000e10,
],
[
6.500000e09,
6.300000e09,
5.900000e09,
5.200000e09,
4.900000e09,
5.000000e09,
4.900000e09,
4.700000e09,
4.500000e09,
4.500000e09,
5.100000e09,
5.800000e09,
6.000000e09,
7.200000e09,
9.400000e09,
1.200000e10,
1.100000e10,
],
[
7.200000e09,
6.700000e09,
6.300000e09,
5.900000e09,
5.800000e09,
6.100000e09,
6.200000e09,
6.100000e09,
5.800000e09,
5.600000e09,
5.700000e09,
5.800000e09,
6.200000e09,
6.600000e09,
7.000000e09,
7.400000e09,
6.700000e09,
],
[
7.000000e09,
6.700000e09,
6.600000e09,
6.400000e09,
6.400000e09,
6.600000e09,
6.700000e09,
6.800000e09,
6.800000e09,
6.600000e09,
6.200000e09,
6.100000e09,
6.100000e09,
6.100000e09,
5.900000e09,
5.900000e09,
4.200000e09,
],
[
5.600000e09,
5.200000e09,
5.700000e09,
5.800000e09,
6.000000e09,
6.000000e09,
6.100000e09,
6.400000e09,
6.500000e09,
6.400000e09,
6.200000e09,
6.000000e09,
5.800000e09,
5.700000e09,
5.500000e09,
5.300000e09,
2.200000e09,
],
[
3.500000e09,
4.100000e09,
4.200000e09,
4.400000e09,
4.500000e09,
4.600000e09,
4.600000e09,
4.700000e09,
4.800000e09,
4.700000e09,
4.600000e09,
4.500000e09,
4.400000e09,
4.200000e09,
3.900000e09,
3.000000e09,
7.000000e08,
],
[
1.000000e09,
1.400000e09,
1.700000e09,
1.900000e09,
2.000000e09,
2.200000e09,
2.400000e09,
2.300000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.100000e09,
2.000000e09,
1.700000e09,
1.300000e09,
8.100000e08,
8.100000e08,
],
[
2.000000e08,
3.000000e08,
3.800000e08,
4.500000e08,
5.100000e08,
5.600000e08,
6.300000e08,
6.700000e08,
6.800000e08,
6.500000e08,
6.000000e08,
5.300000e08,
4.800000e08,
4.000000e08,
1.400000e08,
1.700000e08,
7.400000e07,
],
],
[
[
3.400000e10,
3.600000e10,
3.900000e10,
4.100000e10,
4.300000e10,
4.200000e10,
4.100000e10,
3.900000e10,
3.800000e10,
3.800000e10,
3.800000e10,
3.800000e10,
3.700000e10,
3.600000e10,
3.400000e10,
3.200000e10,
3.100000e10,
],
[
4.600000e10,
4.900000e10,
5.300000e10,
5.600000e10,
5.800000e10,
5.700000e10,
5.500000e10,
5.300000e10,
5.200000e10,
5.100000e10,
5.200000e10,
5.200000e10,
5.200000e10,
5.000000e10,
4.700000e10,
4.500000e10,
4.300000e10,
],
[
6.700000e10,
7.100000e10,
7.700000e10,
8.200000e10,
8.400000e10,
8.200000e10,
7.900000e10,
7.600000e10,
7.400000e10,
7.400000e10,
7.500000e10,
7.600000e10,
7.600000e10,
7.400000e10,
7.000000e10,
6.600000e10,
6.400000e10,
],
[
1.000000e11,
1.100000e11,
1.200000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.000000e11,
],
[
1.550000e11,
1.600000e11,
1.800000e11,
2.900000e11,
2.000000e11,
2.100000e11,
2.200000e11,
2.200000e11,
2.200000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.800000e11,
1.700000e11,
1.600000e11,
1.500000e11,
1.400000e11,
],
[
2.200000e11,
2.300000e11,
2.500000e11,
2.800000e11,
3.000000e11,
3.200000e11,
3.300000e11,
3.300000e11,
3.300000e11,
3.200000e11,
3.100000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.300000e11,
2.100000e11,
2.000000e11,
],
[
2.900000e11,
3.000000e11,
3.300000e11,
3.600000e11,
3.800000e11,
4.000000e11,
4.000000e11,
4.000000e11,
3.900000e11,
3.900000e11,
3.900000e11,
3.700000e11,
3.500000e11,
3.200000e11,
3.000000e11,
2.800000e11,
2.800000e11,
],
[
2.800000e11,
3.000000e11,
3.200000e11,
3.400000e11,
3.600000e11,
3.700000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.600000e11,
3.500000e11,
3.400000e11,
3.200000e11,
3.100000e11,
3.000000e11,
3.100000e11,
],
[
1.600000e11,
1.700000e11,
1.800000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
1.000000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.100000e11,
],
[
3.700000e10,
3.900000e10,
4.200000e10,
4.500000e10,
4.600000e10,
4.600000e10,
4.500000e10,
4.200000e10,
3.900000e10,
4.100000e10,
4.400000e10,
4.700000e10,
5.000000e10,
5.300000e10,
5.500000e10,
5.800000e10,
6.100000e10,
],
[
0.000000e00,
0.000000e00,
3.000000e10,
3.200000e10,
3.100000e10,
2.700000e10,
2.300000e10,
2.500000e10,
2.500000e10,
2.500000e10,
2.600000e10,
3.000000e10,
3.100000e10,
2.900000e10,
2.800000e10,
2.500000e10,
1.800000e10,
],
[
0.000000e00,
0.000000e00,
1.600000e10,
1.100000e10,
7.800000e09,
6.200000e09,
5.700000e09,
5.500000e09,
5.500000e09,
5.400000e09,
6.000000e09,
7.500000e09,
9.200000e09,
1.000000e10,
1.200000e10,
1.300000e10,
1.200000e10,
],
[
0.000000e00,
0.000000e00,
4.200000e09,
5.800000e09,
5.300000e09,
5.700000e09,
5.000000e09,
4.600000e09,
4.000000e09,
6.700000e09,
6.300000e09,
5.900000e09,
6.200000e09,
6.050000e09,
5.900000e09,
7.800000e09,
8.700000e09,
],
[
0.000000e00,
0.000000e00,
7.400000e09,
6.600000e09,
6.300000e09,
6.500000e09,
6.300000e09,
5.900000e09,
5.700000e09,
5.600000e09,
5.900000e09,
6.000000e09,
6.200000e09,
6.400000e09,
1.500000e10,
6.700000e09,
1.900000e10,
],
[
0.000000e00,
0.000000e00,
7.300000e09,
6.900000e09,
6.600000e09,
6.700000e09,
6.800000e09,
6.700000e09,
6.700000e09,
6.600000e09,
6.400000e09,
6.400000e09,
6.500000e09,
6.500000e09,
6.400000e09,
6.200000e09,
6.000000e09,
],
[
0.000000e00,
0.000000e00,
6.600000e09,
6.400000e09,
6.100000e09,
6.100000e09,
6.200000e09,
6.400000e09,
6.500000e09,
6.400000e09,
6.100000e09,
6.100000e09,
6.100000e09,
6.000000e09,
5.800000e09,
5.600000e09,
5.400000e09,
],
[
3.000000e09,
5.200000e09,
4.500000e09,
4.500000e09,
4.700000e09,
4.600000e09,
4.600000e09,
4.700000e09,
4.000000e09,
8.700000e09,
4.300000e09,
4.450000e09,
4.600000e09,
4.300000e09,
4.100000e09,
3.700000e09,
2.800000e09,
],
[
6.100000e08,
1.400000e09,
1.600000e09,
1.800000e09,
2.000000e09,
2.100000e09,
2.200000e09,
2.200000e09,
2.200000e09,
2.300000e09,
2.200000e09,
2.200000e09,
2.200000e09,
2.000000e09,
5.100000e09,
1.200000e09,
2.400000e09,
],
[
1.100000e08,
2.400000e08,
3.100000e08,
3.700000e08,
4.300000e08,
4.000000e08,
5.600000e08,
6.200000e08,
6.400000e08,
6.500000e08,
6.200000e08,
5.900000e08,
5.800000e08,
5.200000e08,
4.600000e08,
2.800000e08,
1.800000e08,
],
],
[
[
3.300000e10,
3.500000e10,
3.800000e10,
4.000000e10,
4.100000e10,
4.000000e10,
3.900000e10,
3.700000e10,
3.600000e10,
3.500000e10,
3.500000e10,
3.500000e10,
3.400000e10,
3.300000e10,
3.100000e10,
2.900000e10,
2.800000e10,
],
[
4.600000e10,
4.800000e10,
5.200000e10,
5.400000e10,
5.600000e10,
5.500000e10,
5.200000e10,
5.000000e10,
4.800000e10,
4.800000e10,
4.800000e10,
4.800000e10,
4.800000e10,
4.600000e10,
4.300000e10,
4.100000e10,
3.900000e10,
],
[
6.600000e10,
7.000000e10,
7.400000e10,
7.800000e10,
8.000000e10,
7.800000e10,
7.500000e10,
7.200000e10,
7.000000e10,
6.900000e10,
7.000000e10,
7.100000e10,
7.100000e10,
6.900000e10,
6.500000e10,
6.100000e10,
5.800000e10,
],
[
1.000000e11,
1.100000e11,
1.100000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
9.800000e10,
9.300000e10,
],
[
1.500000e11,
1.600000e11,
1.800000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.800000e11,
1.800000e11,
1.600000e11,
1.500000e11,
1.400000e11,
1.300000e11,
],
[
2.200000e11,
2.300000e11,
2.500000e11,
2.800000e11,
3.000000e11,
3.100000e11,
3.200000e11,
3.200000e11,
3.200000e11,
3.100000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.300000e11,
2.100000e11,
1.000000e11,
1.800000e11,
],
[
2.700000e11,
2.900000e11,
3.200000e11,
3.500000e11,
3.700000e11,
3.900000e11,
3.900000e11,
3.900000e11,
3.800000e11,
3.800000e11,
3.800000e11,
3.600000e11,
3.300000e11,
3.000000e11,
2.700000e11,
2.600000e11,
2.500000e11,
],
[
2.500000e11,
2.700000e11,
3.000000e11,
3.300000e11,
3.500000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.500000e11,
3.500000e11,
3.400000e11,
3.250000e11,
3.100000e11,
2.900000e11,
2.900000e11,
3.000000e11,
],
[
1.300000e11,
1.400000e11,
1.600000e11,
1.700000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
1.900000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.300000e11,
],
[
2.800000e10,
3.100000e10,
3.400000e10,
3.900000e10,
4.300000e10,
4.500000e10,
4.600000e10,
4.300000e10,
4.000000e10,
4.000000e10,
4.500000e10,
4.900000e10,
5.400000e10,
5.800000e10,
6.300000e10,
7.000000e10,
7.900000e10,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
3.700000e10,
3.100000e10,
2.800000e10,
6.900000e10,
3.200000e10,
3.300000e10,
3.000000e10,
2.700000e10,
2.300000e10,
1.900000e10,
1.600000e10,
1.200000e10,
9.900000e09,
7.400000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
1.400000e10,
9.100000e09,
6.800000e09,
6.100000e09,
3.800000e09,
5.800000e09,
5.700000e09,
5.700000e09,
6.200000e09,
5.800000e09,
5.900000e09,
5.400000e09,
5.100000e09,
5.300000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
7.800000e09,
5.800000e09,
5.400000e09,
4.800000e09,
4.300000e09,
4.500000e09,
4.700000e09,
5.200000e09,
5.200000e09,
5.500000e09,
5.800000e09,
6.000000e09,
6.400000e09,
6.900000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
7.200000e09,
7.000000e09,
6.500000e09,
6.300000e09,
5.900000e09,
5.700000e09,
6.000000e09,
6.200000e09,
6.100000e09,
6.100000e09,
6.400000e09,
6.800000e09,
6.800000e09,
6.700000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
5.900000e09,
7.000000e09,
6.800000e09,
7.000000e09,
7.000000e09,
6.900000e09,
6.900000e09,
6.800000e09,
6.700000e09,
6.500000e09,
6.700000e09,
6.700000e09,
6.600000e09,
6.500000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
3.200000e09,
6.400000e09,
6.000000e09,
6.200000e09,
6.500000e09,
6.500000e09,
6.300000e09,
6.200000e09,
6.200000e09,
6.200000e09,
6.100000e09,
5.900000e09,
5.600000e09,
5.200000e09,
],
[
0.000000e00,
3.400000e09,
5.400000e09,
4.400000e09,
4.500000e09,
4.200000e09,
3.900000e09,
8.500000e09,
4.600000e09,
4.600000e09,
4.700000e09,
4.700000e09,
4.500000e09,
4.300000e09,
4.000000e09,
3.600000e09,
3.100000e09,
],
[
0.000000e00,
6.800000e08,
1.400000e09,
1.500000e09,
1.800000e09,
2.300000e09,
2.400000e09,
2.000000e09,
2.100000e09,
2.200000e09,
2.200000e09,
2.200000e09,
2.200000e09,
2.000000e09,
1.800000e09,
1.300000e09,
1.100000e09,
],
[
0.000000e00,
1.300000e08,
2.400000e08,
2.900000e08,
3.500000e08,
4.600000e08,
5.400000e08,
5.500000e08,
5.900000e08,
6.200000e08,
6.100000e08,
6.000000e08,
6.000000e08,
5.500000e08,
4.700000e08,
3.700000e08,
2.700000e08,
],
],
[
[
3.200000e10,
3.300000e10,
3.500000e10,
3.700000e10,
3.800000e10,
3.700000e10,
3.600000e10,
3.400000e10,
3.300000e10,
3.200000e10,
3.200000e10,
3.200000e10,
3.100000e10,
3.000000e10,
2.800000e10,
2.600000e10,
2.500000e10,
],
[
4.400000e10,
4.600000e10,
4.900000e10,
5.100000e10,
5.200000e10,
5.100000e10,
4.900000e10,
4.600000e10,
4.400000e10,
4.400000e10,
4.400000e10,
4.400000e10,
4.400000e10,
4.200000e10,
4.000000e10,
3.700000e10,
3.500000e10,
],
[
6.300000e10,
6.600000e10,
7.000000e10,
7.400000e10,
7.500000e10,
7.300000e10,
7.000000e10,
6.600000e10,
6.400000e10,
6.400000e10,
6.500000e10,
6.600000e10,
6.600000e10,
6.400000e10,
6.000000e10,
5.600000e10,
5.300000e10,
],
[
9.800000e10,
1.000000e11,
1.100000e11,
1.100000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
9.900000e10,
9.200000e10,
8.700000e10,
],
[
1.500000e11,
1.600000e11,
1.700000e11,
1.800000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
1.000000e11,
1.500000e11,
1.700000e11,
1.700000e11,
1.600000e11,
1.400000e11,
1.300000e11,
1.200000e11,
],
[
2.100000e11,
2.200000e11,
2.500000e11,
2.700000e11,
2.800000e11,
3.000000e11,
3.000000e11,
3.000000e11,
3.000000e11,
2.900000e11,
2.800000e11,
2.600000e11,
2.800000e11,
2.200000e11,
2.000000e11,
1.800000e11,
1.700000e11,
],
[
2.600000e11,
2.800000e11,
3.100000e11,
3.400000e11,
3.600000e11,
3.700000e11,
3.800000e11,
3.700000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.400000e11,
3.200000e11,
2.900000e11,
2.600000e11,
2.400000e11,
2.300000e11,
],
[
2.300000e11,
2.500000e11,
2.800000e11,
3.100000e11,
3.400000e11,
3.500000e11,
3.500000e11,
3.500000e11,
3.400000e11,
3.800000e11,
3.400000e11,
3.400000e11,
3.300000e11,
3.000000e11,
2.800000e11,
2.800000e11,
2.900000e11,
],
[
1.200000e11,
1.300000e11,
1.400000e11,
1.600000e11,
1.500000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.100000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.400000e11,
],
[
2.400000e10,
2.700000e10,
3.000000e10,
3.600000e10,
4.100000e10,
4.500000e10,
4.600000e10,
4.400000e10,
4.000000e10,
4.000000e10,
4.500000e10,
5.100000e10,
5.700000e10,
6.100000e10,
6.600000e10,
7.600000e10,
8.900000e10,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
3.100000e10,
2.600000e10,
2.400000e10,
2.500000e10,
2.500000e10,
2.300000e10,
2.000000e10,
1.700000e10,
1.300000e10,
9.600000e09,
6.200000e09,
4.800000e09,
4.400000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
9.400000e09,
6.800000e09,
5.700000e09,
5.100000e09,
5.100000e09,
5.000000e09,
5.100000e09,
5.100000e09,
4.900000e09,
4.700000e09,
4.500000e09,
4.300000e09,
4.300000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
6.400000e09,
5.700000e09,
4.700000e09,
4.200000e09,
4.200000e09,
4.400000e09,
4.500000e09,
4.900000e09,
4.900000e09,
5.100000e09,
5.500000e09,
5.900000e09,
6.200000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
7.700000e09,
6.500000e09,
6.100000e09,
6.000000e09,
6.000000e09,
6.200000e09,
6.400000e09,
6.100000e09,
6.200000e09,
6.700000e09,
7.200000e09,
7.300000e09,
7.200000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
7.400000e09,
7.000000e09,
7.000000e09,
7.200000e09,
7.300000e09,
7.400000e09,
7.400000e09,
7.100000e09,
7.000000e09,
7.200000e09,
7.200000e09,
6.900000e09,
6.500000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
6.500000e09,
5.900000e09,
6.100000e09,
6.400000e09,
6.500000e09,
6.400000e09,
6.300000e09,
6.200000e09,
6.100000e09,
6.000000e09,
5.600000e09,
5.300000e09,
4.800000e09,
],
[
0.000000e00,
0.000000e00,
4.400000e09,
5.400000e09,
4.400000e09,
4.300000e09,
4.400000e09,
4.500000e09,
4.600000e09,
4.700000e09,
4.700000e09,
4.600000e09,
4.300000e09,
4.100000e09,
3.700000e09,
3.300000e09,
2.900000e09,
],
[
0.000000e00,
0.000000e00,
9.900000e08,
1.500000e09,
1.600000e09,
1.700000e09,
1.800000e09,
1.900000e09,
2.100000e09,
2.100000e09,
2.200000e09,
2.200000e09,
2.100000e09,
1.900000e09,
1.700000e09,
1.400000e09,
1.100000e09,
],
[
0.000000e00,
0.000000e00,
2.100000e08,
2.600000e08,
3.000000e08,
3.700000e08,
4.400000e08,
5.000000e08,
5.500000e08,
5.800000e08,
5.800000e08,
5.800000e08,
5.600000e08,
5.200000e08,
4.600000e08,
3.800000e08,
3.000000e08,
],
],
[
[
3.100000e10,
3.300000e10,
3.500000e10,
3.600000e10,
3.700000e10,
3.600000e10,
3.500000e10,
3.300000e10,
3.200000e10,
3.100000e10,
3.100000e10,
3.100000e10,
3.000000e10,
2.900000e10,
2.700000e10,
2.600000e10,
2.400000e10,
],
[
4.300000e10,
4.500000e10,
4.800000e10,
5.000000e10,
5.100000e10,
5.000000e10,
4.700000e10,
4.500000e10,
4.500000e10,
4.300000e10,
4.300000e10,
4.300000e10,
4.200000e10,
4.100000e10,
3.900000e10,
3.600000e10,
3.400000e10,
],
[
6.200000e10,
6.500000e10,
6.900000e10,
7.200000e10,
7.300000e10,
7.100000e10,
6.800000e10,
6.500000e10,
6.300000e10,
6.200000e10,
6.300000e10,
6.400000e10,
6.400000e10,
6.200000e10,
5.900000e10,
5.500000e10,
5.200000e10,
],
[
9.700000e10,
1.000000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.000000e11,
1.000000e11,
1.000000e11,
1.100000e11,
1.100000e11,
1.000000e11,
9.700000e10,
9.000000e10,
8.500000e10,
],
[
1.500000e11,
1.500000e11,
1.700000e11,
1.900000e11,
1.900000e11,
1.450000e11,
1.000000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.800000e11,
1.700000e11,
1.600000e11,
1.500000e11,
1.400000e11,
1.300000e11,
1.200000e11,
],
[
2.100000e11,
2.200000e11,
2.400000e11,
2.600000e11,
2.800000e11,
2.900000e11,
3.000000e11,
3.000000e11,
2.900000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.300000e11,
2.100000e11,
1.900000e11,
1.700000e11,
1.600000e11,
],
[
2.600000e11,
2.800000e11,
3.100000e11,
3.300000e11,
3.800000e11,
3.700000e11,
3.700000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.400000e11,
3.100000e11,
2.500000e11,
2.500000e11,
2.400000e11,
2.300000e11,
],
[
2.300000e11,
2.500000e11,
2.800000e11,
3.100000e11,
3.300000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.300000e11,
3.300000e11,
3.400000e11,
3.400000e11,
3.200000e11,
3.000000e11,
2.800000e11,
2.700000e11,
2.800000e11,
],
[
1.200000e11,
1.350000e11,
1.500000e11,
1.650000e11,
1.800000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.900000e11,
1.800000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.100000e11,
2.300000e11,
],
[
2.500000e10,
2.700000e10,
3.100000e10,
3.600000e10,
4.100000e10,
4.400000e10,
4.500000e10,
4.300000e10,
3.900000e10,
3.900000e10,
4.400000e10,
4.900000e10,
5.500000e10,
6.000000e10,
6.500000e10,
7.400000e10,
8.600000e10,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
2.500000e10,
2.500000e10,
2.500000e10,
2.300000e10,
2.100000e10,
2.000000e10,
2.000000e10,
1.800000e10,
1.700000e10,
1.300000e10,
8.900000e09,
5.600000e09,
4.600000e09,
4.300000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
1.200000e10,
8.800000e09,
6.600000e09,
5.500000e09,
4.900000e09,
5.200000e09,
5.300000e09,
5.300000e09,
5.100000e09,
4.800000e09,
4.800000e09,
4.200000e09,
4.100000e09,
4.700000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
7.600000e09,
5.800000e09,
5.400000e09,
4.600000e09,
4.400000e09,
4.400000e09,
4.300000e09,
4.100000e09,
4.300000e09,
4.900000e09,
5.200000e09,
5.300000e09,
5.800000e09,
6.200000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
7.900000e09,
7.100000e09,
6.400000e09,
6.100000e09,
6.000000e09,
6.100000e09,
6.300000e09,
6.300000e09,
6.200000e09,
6.300000e09,
6.600000e09,
7.100000e09,
7.500000e09,
7.400000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
6.900000e09,
6.700000e09,
6.800000e09,
7.100000e09,
7.300000e09,
7.500000e09,
7.500000e09,
7.300000e09,
7.100000e09,
7.200000e09,
7.200000e09,
7.100000e09,
6.600000e09,
],
[
0.000000e00,
0.000000e00,
0.000000e00,
0.000000e00,
6.100000e09,
5.900000e09,
6.100000e09,
6.400000e09,
6.600000e09,
6.600000e09,
6.400000e09,
6.300000e09,
6.000000e09,
5.900000e09,
5.600000e09,
5.200000e09,
4.700000e09,
],
[
0.000000e00,
0.000000e00,
3.800000e09,
4.800000e09,
4.100000e09,
4.300000e09,
4.400000e09,
4.600000e09,
4.700000e09,
4.700000e09,
4.700000e09,
8.600000e09,
4.400000e09,
4.000000e09,
3.600000e09,
3.300000e09,
2.900000e09,
],
[
0.000000e00,
0.000000e00,
8.800000e08,
1.400000e09,
1.500000e09,
1.600000e09,
1.900000e09,
2.000000e09,
2.100000e09,
2.200000e09,
2.200000e09,
2.100000e09,
2.000000e09,
1.800000e09,
1.600000e09,
1.400000e09,
1.100000e09,
],
[
0.000000e00,
0.000000e00,
1.800000e08,
2.600000e08,
2.900000e08,
3.600000e08,
3.300000e08,
4.800000e08,
5.500000e08,
5.700000e08,
5.700000e08,
5.500000e08,
5.300000e08,
4.900000e08,
4.300000e08,
3.600000e08,
2.900000e08,
],
],
[
[
3.200000e10,
3.400000e10,
3.600000e10,
3.800000e10,
3.900000e10,
3.800000e10,
3.700000e10,
3.500000e10,
3.400000e10,
3.300000e10,
3.300000e10,
3.300000e10,
3.300000e10,
3.100000e10,
2.950000e10,
2.800000e10,
2.700000e10,
],
[
4.400000e10,
4.600000e10,
5.000000e10,
5.200000e10,
5.300000e10,
5.200000e10,
5.000000e10,
4.700000e10,
4.600000e10,
4.500000e10,
4.500000e10,
4.600000e10,
4.500000e10,
4.400000e10,
4.100000e10,
3.900000e10,
3.700000e10,
],
[
6.400000e10,
6.700000e10,
7.200000e10,
7.600000e10,
7.700000e10,
7.500000e10,
7.200000e10,
6.800000e10,
6.600000e10,
6.600000e10,
6.600000e10,
6.700000e10,
6.700000e10,
6.500000e10,
6.200000e10,
5.800000e10,
5.600000e10,
],
[
1.000000e11,
1.000000e11,
1.100000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.100000e11,
1.000000e11,
9.400000e10,
9.000000e10,
],
[
1.500000e11,
1.600000e11,
1.700000e11,
1.500000e11,
1.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
1.900000e11,
1.900000e11,
1.800000e11,
1.700000e11,
1.600000e11,
1.500000e11,
1.300000e11,
1.300000e11,
],
[
2.100000e11,
2.300000e11,
2.500000e11,
2.700000e11,
2.900000e11,
3.000000e11,
3.000000e11,
3.100000e11,
3.000000e11,
3.000000e11,
2.800000e11,
2.600000e11,
2.400000e11,
2.200000e11,
2.000000e11,
1.900000e11,
1.800000e11,
],
[
2.700000e11,
2.900000e11,
3.200000e11,
3.400000e11,
3.600000e11,
3.800000e11,
3.800000e11,
3.800000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.200000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.500000e11,
],
[
2.600000e11,
2.800000e11,
3.000000e11,
3.300000e11,
3.400000e11,
3.500000e11,
3.500000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.400000e11,
3.200000e11,
3.000000e11,
2.900000e11,
2.900000e11,
2.900000e11,
],
[
1.400000e11,
1.500000e11,
1.600000e11,
1.500000e11,
1.900000e11,
1.900000e11,
1.000000e11,
1.900000e11,
1.900000e11,
1.800000e11,
1.900000e11,
1.900000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.200000e11,
],
[
3.100000e10,
3.300000e10,
3.700000e10,
4.100000e10,
4.400000e10,
4.500000e10,
4.500000e10,
4.200000e10,
3.900000e10,
3.900000e10,
4.300000e10,
4.700000e10,
5.200000e10,
5.500000e10,
5.900000e10,
6.500000e10,
7.100000e10,
],
[
0.000000e00,
0.000000e00,
2.900000e10,
2.800000e10,
2.800000e10,
2.600000e10,
2.400000e10,
2.100000e10,
2.000000e10,
3.100000e10,
2.100000e10,
1.000000e10,
1.400000e10,
1.000000e10,
7.000000e09,
5.400000e09,
5.400000e09,
],
[
0.000000e00,
0.000000e00,
1.400000e10,
1.100000e10,
9.200000e09,
7.600000e09,
6.000000e09,
5.300000e09,
5.100000e09,
5.800000e09,
2.900275e09,
3.850138e09,
4.800000e09,
4.700000e09,
4.400000e09,
4.300000e09,
4.800000e09,
],
[
0.000000e00,
0.000000e00,
8.000000e09,
6.600000e09,
6.300000e09,
5.800000e09,
5.000000e09,
4.400000e09,
4.600000e09,
4.800000e09,
4.700000e09,
4.600000e09,
4.800000e09,
5.200000e09,
5.500000e09,
5.600000e09,
6.000000e09,
],
[
0.000000e00,
0.000000e00,
7.400000e09,
7.400000e09,
7.100000e09,
6.300000e09,
6.000000e09,
6.000000e09,
6.100000e09,
6.300000e09,
6.100000e09,
6.100000e09,
6.100000e09,
6.200000e09,
6.600000e09,
7.000000e09,
7.200000e09,
],
[
0.000000e00,
0.000000e00,
7.000000e09,
6.800000e09,
6.600000e09,
6.700000e09,
6.700000e09,
7.000000e09,
7.100000e09,
7.100000e09,
7.100000e09,
7.000000e09,
6.900000e09,
6.800000e09,
6.800000e09,
6.800000e09,
6.700000e09,
],
[
0.000000e00,
0.000000e00,
5.100000e09,
6.000000e09,
5.800000e09,
6.000000e09,
6.200000e09,
6.500000e09,
6.600000e09,
6.500000e09,
6.200000e09,
6.100000e09,
5.900000e09,
5.700000e09,
5.500000e09,
5.300000e09,
4.900000e09,
],
[
0.000000e00,
1.300000e09,
3.200000e09,
3.700000e09,
4.000000e09,
4.300000e09,
4.500000e09,
4.700000e09,
4.800000e09,
4.800000e09,
4.700000e09,
4.700000e09,
4.400000e09,
4.100000e09,
3.700000e09,
3.400000e09,
3.000000e09,
],
[
0.000000e00,
3.000000e08,
1.000000e09,
1.300000e09,
1.500000e09,
1.700000e09,
2.000000e09,
2.100000e09,
2.300000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.000000e09,
1.800000e09,
1.600000e09,
1.300000e09,
1.000000e09,
],
[
0.000000e00,
1.000000e08,
2.500000e08,
2.900000e08,
3.200000e08,
3.900000e08,
4.700000e08,
5.400000e08,
5.900000e08,
6.000000e08,
5.800000e08,
5.500000e08,
5.100000e08,
4.600000e08,
4.000000e08,
3.300000e08,
2.800000e08,
],
],
[
[
5.300000e10,
3.500000e10,
3.800000e10,
4.100000e10,
4.200000e10,
4.200000e10,
4.000000e10,
3.900000e10,
3.800000e10,
3.800000e10,
3.800000e10,
3.500000e10,
3.800000e10,
3.700000e10,
3.500000e10,
3.300000e10,
3.200000e10,
],
[
4.500000e10,
4.800000e10,
5.200000e10,
5.600000e10,
5.800000e10,
5.700000e10,
5.500000e10,
5.300000e10,
5.100000e10,
5.100000e10,
5.200000e10,
5.200000e10,
5.200000e10,
5.000000e10,
4.800000e10,
4.600000e10,
4.400000e10,
],
[
6.600000e10,
7.100000e10,
7.600000e10,
8.100000e10,
8.300000e10,
8.200000e10,
7.900000e10,
7.500000e10,
7.400000e10,
7.300000e10,
7.500000e10,
7.600000e10,
7.600000e10,
7.400000e10,
7.100000e10,
6.700000e10,
6.500000e10,
],
[
1.000000e11,
1.170000e11,
1.200000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.200000e11,
1.100000e11,
1.100000e11,
1.000000e11,
],
[
1.500000e11,
1.600000e11,
1.800000e11,
1.000000e11,
2.000000e11,
2.100000e11,
2.200000e11,
2.200000e11,
2.200000e11,
2.100000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.500000e11,
1.600000e11,
1.500000e11,
1.500000e11,
],
[
2.200000e11,
2.300000e11,
2.500000e11,
2.800000e11,
3.000000e11,
3.200000e11,
3.300000e11,
3.300000e11,
3.300000e11,
3.200000e11,
3.100000e11,
2.900000e11,
2.700000e11,
2.500000e11,
2.300000e11,
2.200000e11,
2.100000e11,
],
[
2.900000e11,
3.000000e11,
3.300000e11,
3.600000e11,
3.800000e11,
4.000000e11,
4.100000e11,
4.000000e11,
4.000000e11,
3.900000e11,
3.900000e11,
3.800000e11,
3.500000e11,
3.300000e11,
3.100000e11,
2.900000e11,
2.900000e11,
],
[
3.000000e11,
3.100000e11,
3.300000e11,
3.500000e11,
3.700000e11,
3.700000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.300000e11,
3.200000e11,
3.100000e11,
3.100000e11,
],
[
1.800000e11,
1.800000e11,
2.900000e11,
2.000000e11,
2.100000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
1.900000e11,
1.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
2.000000e11,
],
[
4.300000e10,
4.500000e10,
4.700000e10,
4.900000e10,
4.900000e10,
4.800000e10,
4.600000e10,
4.300000e10,
4.000000e10,
4.100000e10,
4.500000e10,
4.700000e10,
4.900000e10,
5.100000e10,
5.300000e10,
5.400000e10,
5.500000e10,
],
[
2.500000e10,
2.800000e10,
3.100000e10,
3.300000e10,
3.200000e10,
2.900000e10,
2.500000e10,
2.200000e10,
2.100000e10,
2.000000e10,
2.200000e10,
2.300000e10,
2.100000e10,
1.900000e10,
1.500000e10,
1.500000e10,
1.200000e10,
],
[
1.400000e10,
1.600000e10,
1.400000e10,
1.200000e10,
9.400000e09,
7.300000e09,
6.200000e09,
5.500000e09,
5.500000e09,
5.600000e09,
5.700000e09,
6.000000e09,
6.000000e09,
6.100000e09,
5.700000e09,
6.300000e09,
7.500000e09,
],
[
7.800000e09,
8.100000e09,
7.200000e09,
6.700000e09,
6.300000e09,
5.800000e09,
5.200000e09,
4.400000e09,
4.300000e09,
4.600000e09,
4.700000e09,
4.600000e09,
4.750000e09,
4.900000e09,
5.400000e09,
5.600000e09,
6.100000e09,
],
[
6.000000e09,
7.300000e09,
7.300000e09,
7.300000e09,
6.900000e09,
6.100000e09,
5.700000e09,
5.600000e09,
5.800000e09,
5.800000e09,
5.700000e09,
5.600000e09,
5.500000e09,
5.600000e09,
5.700000e09,
6.000000e09,
6.500000e09,
],
[
0.000000e00,
7.000000e09,
6.800000e09,
6.700000e09,
6.600000e09,
6.500000e09,
6.500000e09,
6.700000e09,
6.700000e09,
6.600000e09,
6.700000e09,
6.500000e09,
6.300000e09,
6.350000e09,
6.400000e09,
6.200000e09,
6.100000e09,
],
[
0.000000e00,
5.800000e09,
5.800000e09,
5.800000e09,
5.800000e09,
6.000000e09,
6.100000e09,
6.000000e09,
6.400000e09,
6.200000e09,
6.100000e09,
6.000000e09,
5.800000e09,
5.650000e09,
5.500000e09,
5.400000e09,
5.300000e09,
],
[
1.000000e09,
3.100000e09,
3.700000e09,
4.000000e09,
4.300000e09,
4.500000e09,
4.500000e09,
4.600000e09,
4.700000e09,
4.800000e09,
4.700000e09,
4.700000e09,
4.900000e09,
4.500000e09,
4.100000e09,
3.800000e09,
3.100000e09,
],
[
2.500000e08,
1.000000e09,
1.400000e09,
1.600000e09,
1.800000e09,
2.000000e09,
2.100000e09,
2.200000e09,
2.300000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.100000e09,
3.000000e08,
1.600000e09,
1.300000e09,
9.000000e08,
],
[
1.100000e08,
2.900000e08,
3.700000e08,
3.700000e08,
4.000000e08,
4.700000e08,
5.400000e08,
6.000000e08,
6.400000e08,
6.300000e08,
6.000000e08,
5.400000e08,
5.000000e08,
1.600000e08,
3.600000e08,
2.700000e08,
1.800000e08,
],
],
[
[
3.200000e10,
3.400000e10,
3.800000e10,
4.200000e10,
4.400000e10,
4.400000e10,
4.300000e10,
4.200000e10,
4.200000e10,
4.200000e10,
4.300000e10,
4.400000e10,
4.400000e10,
4.300000e10,
4.100000e10,
4.000000e10,
3.900000e10,
],
[
4.400000e10,
4.800000e10,
5.300000e10,
5.700000e10,
6.000000e10,
6.000000e10,
5.800000e10,
5.700000e10,
5.600000e10,
5.700000e10,
5.800000e10,
5.900000e10,
6.000000e10,
5.900000e10,
5.700000e10,
5.400000e10,
5.300000e10,
],
[
6.600000e10,
7.100000e10,
7.800000e10,
8.400000e10,
8.700000e10,
8.700000e10,
8.400000e10,
8.200000e10,
8.000000e10,
8.100000e10,
8.300000e10,
8.500000e10,
8.600000e10,
8.500000e10,
8.200000e10,
7.900000e10,
7.700000e10,
],
[
1.000000e11,
1.100000e11,
1.200000e11,
1.300000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.200000e11,
1.200000e11,
],
[
1.500000e11,
1.600000e11,
1.800000e11,
2.000000e11,
2.100000e11,
2.200000e11,
2.300000e11,
2.400000e11,
2.000000e11,
2.400000e11,
2.300000e11,
2.200000e11,
2.100000e11,
2.000000e11,
1.900000e11,
1.800000e11,
1.800000e11,
],
[
2.100000e11,
2.300000e11,
2.500000e11,
2.800000e11,
3.100000e11,
3.300000e11,
3.500000e11,
3.600000e11,
3.600000e11,
3.600000e11,
3.500000e11,
3.300000e11,
3.100000e11,
2.900000e11,
2.700000e11,
2.600000e11,
2.500000e11,
],
[
2.900000e11,
3.000000e11,
3.300000e11,
3.600000e11,
4.000000e11,
4.200000e11,
4.300000e11,
4.300000e11,
4.300000e11,
4.300000e11,
4.300000e11,
4.200000e11,
4.000000e11,
3.700000e11,
3.500000e11,
3.400000e11,
3.300000e11,
],
[
3.200000e11,
3.300000e11,
3.400000e11,
3.700000e11,
3.900000e11,
4.000000e11,
4.000000e11,
3.900000e11,
3.900000e11,
3.900000e11,
3.900000e11,
3.900000e11,
3.800000e11,
3.600000e11,
3.500000e11,
3.300000e11,
3.200000e11,
],
[
2.200000e11,
2.200000e11,
2.200000e11,
2.200000e11,
2.300000e11,
2.200000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.100000e11,
2.000000e11,
2.000000e11,
1.000000e11,
1.500000e11,
],
[
6.300000e10,
6.200000e10,
6.000000e10,
5.900000e10,
5.700000e10,
5.300000e10,
4.900000e10,
4.500000e10,
4.300000e10,
4.600000e10,
4.900000e10,
4.900000e10,
4.900000e10,
4.800000e10,
4.700000e10,
4.500000e10,
4.300000e10,
],
[
1.500000e10,
2.400000e10,
2.800000e10,
2.900000e10,
2.800000e10,
2.500000e10,
2.300000e10,
2.400000e10,
2.500000e10,
2.300000e10,
2.100000e10,
2.400000e10,
2.800000e10,
2.800000e10,
2.600000e10,
2.300000e10,
1.600000e10,
],
[
9.100000e09,
8.900000e09,
9.000000e09,
8.800000e09,
7.600000e09,
6.600000e09,
6.200000e09,
6.500000e09,
6.600000e09,
6.000000e09,
6.300000e09,
6.600000e09,
7.300000e09,
9.700000e09,
1.400000e10,
1.600000e10,
1.500000e10,
],
[
8.000000e09,
7.700000e09,
7.700000e09,
7.100000e09,
6.350000e09,
5.600000e09,
5.300000e09,
4.700000e09,
4.400000e09,
4.400000e09,
4.700000e09,
4.800000e09,
5.100000e09,
5.200000e09,
7.000000e09,
1.000000e10,
1.100000e10,
],
[
7.000000e09,
6.900000e09,
7.000000e09,
6.800000e09,
6.400000e09,
6.000000e09,
5.700000e09,
5.300000e09,
5.100000e09,
5.100000e09,
6.100000e09,
6.000000e09,
6.200000e09,
6.300000e09,
5.700000e09,
7.300000e09,
0.000000e00,
],
[
6.800000e09,
6.700000e09,
6.600000e09,
6.600000e09,
6.500000e09,
6.500000e09,
6.400000e09,
6.500000e09,
6.500000e09,
6.400000e09,
6.500000e09,
6.300000e09,
6.100000e09,
6.000000e09,
1.000000e10,
6.700000e09,
0.000000e00,
],
[
5.100000e09,
5.600000e09,
6.200000e09,
5.800000e09,
5.900000e09,
6.000000e09,
6.100000e09,
6.300000e09,
6.500000e09,
6.400000e09,
6.800000e09,
6.000000e09,
5.900000e09,
6.000000e09,
6.000000e09,
6.000000e09,
0.000000e00,
],
[
2.700000e09,
3.600000e09,
4.000000e09,
4.300000e09,
4.500000e09,
4.600000e09,
4.600000e09,
4.700000e09,
4.700000e09,
4.700000e09,
4.600000e09,
4.600000e09,
4.600000e09,
4.500000e09,
4.400000e09,
3.900000e09,
3.100000e09,
],
[
9.900000e08,
1.400000e09,
1.600000e09,
1.900000e09,
2.100000e09,
2.200000e09,
2.200000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.200000e09,
2.100000e09,
2.000000e09,
1.800000e09,
1.500000e09,
1.100000e09,
6.500000e08,
],
[
3.100000e08,
4.100000e08,
4.300000e08,
4.600000e08,
5.100000e08,
5.500000e08,
6.000000e08,
6.300000e08,
6.400000e08,
6.200000e08,
5.700000e08,
5.100000e08,
4.500000e08,
3.800000e08,
3.000000e08,
2.000000e08,
1.000000e08,
],
],
[
[
2.900000e10,
3.100000e10,
3.500000e10,
3.900000e10,
4.100000e10,
4.200000e10,
4.200000e10,
4.100000e10,
4.200000e10,
4.300000e10,
4.500000e10,
4.600000e10,
4.700000e10,
4.600000e10,
4.600000e10,
4.500000e10,
4.400000e10,
],
[
4.100000e10,
4.400000e10,
4.900000e10,
5.400000e10,
5.700000e10,
5.800000e10,
5.700000e10,
5.700000e10,
5.700000e10,
5.800000e10,
6.100000e10,
6.300000e10,
6.400000e10,
6.400000e10,
6.200000e10,
6.100000e10,
6.000000e10,
],
[
6.100000e10,
6.700000e10,
7.400000e10,
8.100000e10,
8.500000e10,
8.600000e10,
8.400000e10,
8.200000e10,
8.200000e10,
8.400000e10,
8.700000e10,
9.000000e10,
9.200000e10,
9.200000e10,
9.000000e10,
8.800000e10,
8.600000e10,
],
[
9.900000e10,
1.100000e11,
1.200000e11,
1.300000e11,
1.400000e11,
1.600000e11,
1.400000e11,
1.300000e11,
1.300000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.300000e11,
],
[
1.400000e11,
1.500000e11,
1.700000e11,
1.900000e11,
2.100000e11,
2.200000e11,
2.300000e11,
2.400000e11,
2.500000e11,
2.500000e11,
2.400000e11,
2.400000e11,
2.300000e11,
2.200000e11,
2.100000e11,
2.000000e11,
2.000000e11,
],
[
1.900000e11,
2.100000e11,
2.400000e11,
2.700000e11,
3.000000e11,
3.300000e11,
3.500000e11,
3.600000e11,
3.700000e11,
3.700000e11,
3.700000e11,
3.500000e11,
3.400000e11,
3.200000e11,
3.000000e11,
2.900000e11,
2.800000e11,
],
[
2.700000e11,
2.900000e11,
3.100000e11,
3.500000e11,
4.000000e11,
4.300000e11,
4.500000e11,
4.500000e11,
4.500000e11,
4.500000e11,
4.600000e11,
4.500000e11,
4.300000e11,
4.100000e11,
3.900000e11,
3.700000e11,
3.500000e11,
],
[
3.200000e11,
3.200000e11,
3.400000e11,
3.700000e11,
4.000000e11,
4.200000e11,
4.200000e11,
4.200000e11,
4.100000e11,
4.100000e11,
4.200000e11,
4.100000e11,
4.000000e11,
3.800000e11,
3.600000e11,
3.000000e11,
3.200000e11,
],
[
2.500000e11,
2.300000e11,
2.300000e11,
2.400000e11,
2.400000e11,
2.400000e11,
2.300000e11,
2.200000e11,
2.300000e11,
2.300000e11,
2.300000e11,
2.200000e11,
2.200000e11,
2.100000e11,
1.900000e11,
1.800000e11,
1.700000e11,
],
[
8.400000e10,
7.700000e10,
7.100000e10,
6.800000e10,
6.400000e10,
5.800000e10,
5.300000e10,
4.800000e10,
4.600000e10,
5.000000e10,
5.300000e10,
5.200000e10,
5.000000e10,
4.600000e10,
4.200000e10,
3.800000e10,
3.600000e10,
],
[
6.100000e09,
7.400000e09,
9.800000e09,
1.200000e10,
1.600000e10,
1.800000e10,
2.000000e10,
2.400000e10,
2.800000e10,
2.900000e10,
2.900000e10,
2.900000e10,
3.200000e10,
3.300000e10,
2.900000e10,
0.000000e00,
0.000000e00,
],
[
5.700000e09,
5.500000e09,
5.900000e09,
5.500000e09,
5.900000e09,
5.800000e09,
6.000000e09,
6.400000e09,
6.600000e09,
6.500000e09,
6.800000e09,
8.300000e09,
9.900000e09,
1.400000e10,
1.600000e10,
0.000000e00,
0.000000e00,
],
[
6.700000e09,
6.500000e09,
6.400000e09,
6.200000e09,
5.800000e09,
5.400000e09,
4.700000e09,
4.700000e09,
4.400000e09,
4.300000e09,
4.500000e09,
5.300000e09,
6.100000e09,
7.700000e09,
1.100000e10,
0.000000e00,
0.000000e00,
],
[
6.700000e09,
6.600000e09,
6.500000e09,
6.300000e09,
6.000000e09,
5.800000e09,
5.600000e09,
5.400000e09,
5.200000e09,
5.300000e09,
5.900000e09,
6.600000e09,
6.900000e09,
7.400000e09,
8.100000e09,
0.000000e00,
0.000000e00,
],
[
6.600000e09,
6.600000e09,
6.600000e09,
7.200000e09,
6.600000e09,
6.700000e09,
6.600000e09,
6.900000e09,
6.800000e09,
6.700000e09,
6.300000e09,
6.500000e09,
6.200000e09,
6.300000e09,
7.300000e09,
0.000000e00,
0.000000e00,
],
[
5.000000e09,
5.300000e09,
5.700000e09,
5.900000e09,
6.200000e09,
6.200000e09,
6.200000e09,
6.400000e09,
6.600000e09,
6.500000e09,
5.900000e09,
6.100000e09,
6.100000e09,
6.200000e09,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
2.900000e09,
3.400000e09,
3.800000e09,
4.100000e09,
4.400000e09,
4.600000e09,
4.600000e09,
4.600000e09,
4.700000e09,
4.500000e09,
4.100000e09,
5.000000e08,
4.400000e09,
4.500000e09,
4.400000e09,
3.000000e09,
0.000000e00,
],
[
1.100000e09,
1.400000e09,
1.700000e09,
2.000000e09,
2.200000e09,
2.300000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.100000e09,
2.000000e09,
2.000000e09,
1.900000e09,
1.600000e09,
1.200000e09,
6.000000e08,
0.000000e00,
],
[
3.200000e08,
4.000000e08,
4.600000e08,
5.400000e08,
6.100000e08,
6.400000e08,
6.500000e08,
6.500000e08,
6.300000e08,
5.800000e08,
5.300000e08,
4.700000e08,
3.900000e08,
3.100000e08,
2.200000e08,
1.100000e08,
0.000000e00,
],
],
[
[
2.600000e10,
2.800000e10,
3.200000e10,
3.600000e10,
3.800000e10,
3.900000e10,
3.900000e10,
3.900000e10,
3.900000e10,
4.100000e10,
4.300000e10,
4.500000e10,
4.600000e10,
4.600000e10,
4.600000e10,
4.500000e10,
4.500000e10,
],
[
3.700000e10,
4.000000e10,
4.500000e10,
5.000000e10,
5.300000e10,
5.400000e10,
5.300000e10,
5.300000e10,
5.400000e10,
5.600000e10,
5.900000e10,
6.100000e10,
6.300000e10,
6.300000e10,
6.300000e10,
6.200000e10,
6.100000e10,
],
[
5.600000e10,
6.100000e10,
6.900000e10,
7.600000e10,
8.000000e10,
8.000000e10,
7.900000e10,
7.800000e10,
7.800000e10,
8.000000e10,
8.400000e10,
8.800000e10,
9.000000e10,
9.100000e10,
9.000000e10,
8.900000e10,
8.800000e10,
],
[
9.200000e10,
1.000000e11,
1.100000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.300000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
1.400000e11,
],
[
1.300000e11,
1.800000e11,
1.600000e11,
1.900000e11,
2.000000e11,
2.100000e11,
2.200000e11,
2.300000e11,
2.400000e11,
2.400000e11,
2.400000e11,
2.300000e11,
2.300000e11,
2.200000e11,
2.100000e11,
2.100000e11,
2.000000e11,
],
[
2.000000e11,
2.300000e11,
2.600000e11,
2.900000e11,
3.200000e11,
3.300000e11,
3.400000e11,
3.500000e11,
3.600000e11,
3.700000e11,
3.600000e11,
3.600000e11,
3.400000e11,
3.300000e11,
3.100000e11,
3.000000e11,
2.900000e11,
],
[
2.500000e11,
3.000000e11,
3.400000e11,
3.800000e11,
4.200000e11,
4.300000e11,
4.400000e11,
4.400000e11,
4.400000e11,
4.500000e11,
4.600000e11,
4.500000e11,
4.300000e11,
4.100000e11,
3.900000e11,
3.700000e11,
3.500000e11,
],
[
3.100000e11,
3.100000e11,
3.200000e11,
3.600000e11,
3.900000e11,
4.200000e11,
4.200000e11,
4.100000e11,
4.100000e11,
4.200000e11,
4.200000e11,
4.200000e11,
4.000000e11,
3.800000e11,
3.600000e11,
3.300000e11,
3.100000e11,
],
[
2.500000e11,
2.400000e11,
2.300000e11,
2.400000e11,
2.500000e11,
2.400000e11,
2.300000e11,
2.300000e11,
2.300000e11,
2.300000e11,
2.300000e11,
2.300000e11,
2.200000e11,
2.000000e11,
1.500000e11,
1.700000e11,
1.600000e11,
],
[
9.700000e10,
8.500000e10,
7.600000e10,
7.200000e10,
6.800000e10,
6.100000e10,
5.400000e10,
4.800000e10,
4.800000e10,
5.200000e10,
5.500000e10,
5.300000e10,
4.900000e10,
4.400000e10,
3.800000e10,
3.500000e10,
3.200000e10,
],
[
4.600000e09,
5.100000e09,
5.700000e09,
7.800000e09,
1.100000e10,
1.400000e10,
1.800000e10,
2.200000e10,
2.300000e10,
2.400000e10,
2.500000e10,
2.600000e10,
3.100000e10,
3.300000e10,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
5.000000e09,
4.600000e09,
4.800000e09,
4.600000e09,
4.600000e09,
4.600000e09,
5.300000e09,
5.700000e09,
5.800000e09,
5.900000e09,
6.800000e09,
9.000000e09,
1.200000e10,
1.600000e10,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
6.600000e09,
6.300000e09,
5.900000e09,
5.500000e09,
5.300000e09,
5.400000e09,
4.500000e09,
4.800000e09,
4.600000e09,
4.400000e09,
4.700000e09,
5.800000e09,
7.200000e09,
1.000000e10,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
7.100000e09,
7.100000e09,
6.800000e09,
6.500000e09,
6.100000e09,
6.000000e09,
6.100000e09,
6.100000e09,
5.900000e09,
5.800000e09,
5.800000e09,
2.000000e10,
8.100000e09,
8.600000e09,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
6.600000e09,
6.900000e09,
7.000000e09,
7.000000e09,
6.900000e09,
7.200000e09,
7.400000e09,
7.500000e09,
7.400000e09,
7.200000e09,
6.900000e09,
6.800000e09,
6.300000e09,
8.100000e09,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
4.800000e09,
5.200000e09,
5.600000e09,
6.000000e09,
6.200000e09,
6.400000e09,
6.500000e09,
6.800000e09,
6.900000e09,
6.700000e09,
6.200000e09,
5.800000e09,
5.700000e09,
6.600000e09,
0.000000e00,
0.000000e00,
0.000000e00,
],
[
2.900000e09,
3.300000e09,
3.700000e09,
4.000000e09,
4.300000e09,
4.500000e09,
4.700000e09,
4.800000e09,
4.800000e09,
4.600000e09,
4.300000e09,
4.100000e09,
3.800000e09,
3.400000e09,
3.200000e09,
0.000000e00,
0.000000e00,
],
[
1.200000e09,
1.400000e09,
1.700000e09,
1.900000e09,
2.100000e09,
2.200000e09,
2.300000e09,
2.300000e09,
2.200000e09,
2.100000e09,
1.900000e09,
1.700000e09,
1.500000e09,
1.100000e09,
7.100000e08,
0.000000e00,
0.000000e00,
],
[
3.300000e08,
4.000000e08,
4.800000e08,
5.500000e08,
6.100000e08,
6.400000e08,
6.400000e08,
6.300000e08,
6.000000e08,
5.500000e08,
5.000000e08,
4.200000e08,
3.200000e08,
2.200000e08,
1.300000e08,
0.000000e00,
0.000000e00,
],
],
],
dtype=float,
)
def _month_coordinate(time: Any) -> float:
"""Return a fractional month coordinate in the interval [1, 13)."""
if isinstance(time, np.datetime64):
day = time.astype("datetime64[D]").astype(object)
if isinstance(day, int):
msg = "np.datetime64 values must include a calendar date"
raise TypeError(msg)
time = day
if isinstance(time, _dt.datetime | _dt.date):
days_in_month = calendar.monthrange(time.year, time.month)[1]
return time.month + (time.day - 1) / days_in_month
value = float(time)
if 1.0 <= value < 13.0:
return value
msg = "time must be a date/datetime/np.datetime64 or a fractional month in [1, 13)"
raise ValueError(msg)
def _interpolate_month(table: np.ndarray, month_coordinate: float) -> np.ndarray:
lower = int(np.floor(month_coordinate))
fraction = month_coordinate - lower
lower_index = (lower - 1) % 12
upper_index = lower % 12
return (1.0 - fraction) * table[lower_index] + fraction * table[upper_index]
[docs]
def number_density(
altitudes_m: npt.ArrayLike,
time: _dt.datetime | _dt.date | np.datetime64 | float,
latitude_degrees: float,
lower_fill_value: float | None = 0.0,
upper_fill_value: float | None = None,
) -> np.ndarray:
"""Return atomic oxygen number density on an altitude grid.
Parameters
----------
altitudes_m : array-like
Altitudes in metres. Values must be within 40 km to 130 km.
time : datetime, date, np.datetime64, or float
Time used for interpolation. Numeric values are interpreted as fractional
month numbers, where 1 is January and 12 is December. Values between
month integers are linearly interpolated, with December wrapping to January.
latitude_degrees : float
Geographic latitude in degrees. Values must be within -80 to 80 degrees.
lower_fill_value : float or None, optional
Number density in m^-3 to use below 40 km. If None, the 40 km value is
used. The default is 0.
upper_fill_value : float or None, optional
Number density in m^-3 to use above 130 km. If None, the 130 km value is
used. The default is None.
Returns
-------
np.ndarray
Atomic oxygen number density in m^-3, with the same shape as ``altitudes_m``.
"""
altitude_array = np.asarray(altitudes_m, dtype=float)
latitude = float(latitude_degrees)
if latitude < LATITUDES_DEGREES.min() or latitude > LATITUDES_DEGREES.max():
msg = "latitude_degrees must be within the tabulated range -80 to 80 degrees"
raise ValueError(msg)
table_at_time = _interpolate_month(_TABLE_CM3, _month_coordinate(time))
profile_cm3 = np.empty(table_at_time.shape[0], dtype=float)
for altitude_index in range(table_at_time.shape[0]):
profile_cm3[altitude_index] = np.interp(
latitude, LATITUDES_DEGREES, table_at_time[altitude_index]
)
lower_fill_cm3 = (
profile_cm3[-1] if lower_fill_value is None else float(lower_fill_value) / 1.0e6
)
upper_fill_cm3 = (
profile_cm3[0] if upper_fill_value is None else float(upper_fill_value) / 1.0e6
)
density_cm3 = np.interp(
altitude_array,
ALTITUDES_M[::-1],
profile_cm3[::-1],
left=lower_fill_cm3,
right=upper_fill_cm3,
)
return density_cm3 * 1.0e6