Datasets For "Estimating Maximum Extent of Auroral Equatorward Boundary using Historical and Simulated Surface Magnetic Field Data", Blake et al. (2020), JGR

Datasets and sample Python codes for the 2020 paper "Estimating Maximum Extent of Auroral Equatorward Boundary using Historical and Simulated Surface Magnetic Field Data", by Blake et al., submitted to the Journal of Gephysical Research, Space Physics. 

Up-to-date Python codes can be found at https://github.com/TerminusEst/Auroral_Boundary_Geomag

The complete SWMF simulation folders (including parameter and log files etc.) can be requested from NASA's Community Coordinated Modeling Center.

#########

Data/ contains the following:

Data/HIST_DATA.txt contains the minimum Dst values and calculated maximum extents of the auroral equatorward boundaries for 25 years of INTERMAGNET data (1991-2016). The fourth column is the standard deviation of the calculated auroral boundary in degrees. 

Data/Boundary_Fits.csv contains the calculated minimum Dst values, and calculated auroral boundaries using Method 1 and Method 2 (see main paper's ttext), for each of the 15 SWMF simulations. Also included are the uncertainties for each calculation.

Data/SWMF_outputs/ contains 15 .txt files, each of which correspond to an SWMF simulation of the same name given in Table 1 in the main text. These data are for the magnetic longitude, magnetic latitude and maximum calculated EH (V/km) for each simulation.

#########

Codes/ contains two python scripts, and some sample data. These scripts correspond to Section 2 in the main text:

1) Boundary_Calc.py calculates the extent of the auroral boundary using magnetic latitudes and maximum calculated EH values from multiple INTERMAGNET sites.  

2) Efield_Calc.py calculates the E-field for a single INTERMAGNET site using the Quebec 1-D resistivity model.

A more detailed description of these codes can be found here: https://github.com/TerminusEst/Auroral_Boundary_Geomag