EUVpy

A Python wrapper for empirical solar extreme ultraviolet (EUV) models (NEUVAC, EUVAC, HEUVAC, HFG) that generates, compares, and exports EUV spectra for thermosphere-ionosphere modeling.

EUVpy is a Python package that brings the most widely used empirical models of solar extreme ultraviolet (EUV) irradiance together behind a single, consistent interface, primarily so that they can drive thermosphere-ionosphere models. It provides NEUVAC (the Nonlinear Extreme Ultraviolet Irradiance Model for Aeronomic Calculations), EUVAC, HEUVAC, and HFG (also referred to as SOLOMON). Each model is parameterized by the 10.7 cm solar radio flux (F10.7) and quantities derived from it - the 81-day average centred on the current day and a 54-day backward-looking average - so a full EUV spectrum can be produced from a readily available ground-based index alone. Spectra can be generated in several binning schemes, including the 59-bin scheme used by the Global Ionosphere-Thermosphere Model (GITM) and the Aether model, the 37 bins of EUVAC and of Torr, and the 22 overlapping standard bands of Solomon and Qian (2005). Alongside running the models, EUVpy retrieves and cleans the data they need and are judged against. It can download F10.7 from NASA OMNIWeb and flare-corrected, Sun-Earth-distance-adjusted radio flux from Collecte Localisation Satellites, read Penticton F10.7 files, and gap-fill, filter and average those series. It can download and read solar spectral irradiance from FISM2 (through LISIRD, in both native and standard-band form), TIMED/SEE Level 3, SDO/EVE, and NRLSSI2, rebin any of those spectra into the supported band schemes, and convert between spectral irradiance and photon flux. It computes comparison statistics between models and observations, including root-mean-square error, correlation, mean absolute percentage error and percent deviation, and can generate NEUVAC irradiance ensembles from a residual correlation matrix so that quantified uncertainties can be passed to downstream models. Results can be written directly as a GITM EUV input file or as an Aether coefficient table. Two assumptions are worth knowing before use. NEUVAC is the model developed with this package and is the only one for which running ensembles is recommended; the EUVAC, HEUVAC and HFG implementations are included primarily so that NEUVAC can be compared against them. HEUVAC is a Fortran code that must be compiled once with gfortran before it can be called, and the package has so far been tested only on Ubuntu Linux. Users should cite the original authors of whichever model or dataset they use, as set out in the repository's README.