PyIRI

Python implementation of the International Reference Ionosphere (IRI) that evaluates the electron density and associated ionospheric parameters on the entire given global grid and for the entire day simultaneously.

Installation

PyIRI may be installed from PyPi, which will handle all dependencies.

pip install PyIRI

You may also clone and install PyIRI from the GitHub repository.

git clone https://github.com/victoriyaforsythe/PyIRI.git cd /path/to/PyIRI python -m build . pip install .

For more details about the dependencies, see the [examples] and Jupyter Notebooks [tutorials].

Example

PyIRI can calculate monthly mean ionospheric parameters for the user provided grid. The estimation of the parameters occurs simultaneously at all grid points and for all desired diurnal time frames,

``` import numpy as np import PyIRI import PyIRI.edp_update as ml import PyIRI.plotting as plot

year = 2020 month = 4 ccirorursi = 0

Create any horizontal grid (regular or irregular, global or regional).

The grid arrays (alon and alat) should be flattened to be 1-D arrays.

This is an example of a regular global grid:

dlon = 5 dlat = 5 alon2d, alat2d = np.mgrid[-180:180 + dlon:dlon, -90:90 + dlat:dlat] alon = np.reshape(alon2d, alon2d.size) alat = np.reshape(alat2d, alat2d.size)

Create any temporal array expressed in decimal hours (regular or irregular).

For this example we use regularly spaced time array:

hrres = 1 ahr = np.arange(0, 24, hrres)

Find ionospheric parameters for F2, F1, E, and Es regions by

calling IRImonthlymean_par function:

f2, f1, epeak, espeak, sun, mag = ml.IRImonthlymeanpar(year, month, ahr, alon, alat, PyIRI.coeffdir, ccirorursi)

```

PyIRI can calculate daily ionospheric parameters for the user provided grid. The estimation of the parameters occurs simultaneously at all grid points and for all desired diurnal time frames. The parameters are obtained by linearly interpolating between min and max levels of solar activity, and by interpolation between median monthly values to the day of interest. The user should provide the F10.7 value for the day of interest.

f107 = 100 f2, f1, e_peak, es_peak, sun, mag, edp = ml.IRI_density_1day(year, month, day, ahr, alon, alat, aalt, f107, PyIRI.coeff_dir, ccir_or_ursi)

PyIRI does not calculate the Total Electron Content (TEC) automatically, as the altitude array is provided by the user. If the user supplies an array with low resolution, the resulting TEC values may be inaccurate. The edptovtec function can be used to calculate the vTEC from the edp output.

```

Find TEC in shape [Ntime, Ngrid]

TEC = main.edptovtec(edp, aalt, minalt=0.0, maxalt=202000.0)

```

In case you are interested in a single location (as opposed to the grid) PyIRI can evaluate parameters at this location if it is passed as a 1-element NumPy array.

``` alon = np.array([10.]) alat = np.array([20.])

f2, f1, epeak, espeak, sun, mag, edp = ml.IRIdensity1day(year, month, day, ahr, alon, alat, aalt, f107, PyIRI.coeffdir, cciror_ursi) ```