DOI: 10.5281/zenodo.10702853
title: BAS-PRO Model Solution for "Modeling Field Line Curvature Scattering Loss of 1 to 10MeV Protons during Geomagnetic Storms"

description: The file ModelingFLC_BASPRO_solution.zip is a BAS-PRO model solution output archived as a zip file. After extracting the zip file, the solution will be spread across multiple plaintext files. The solution is a grid of proton phase space density multiplied by proton rest mass cubed, f, with units km-6 s3. f is specified in terms of the first, second and third adiabatic invariants μ, K and L as well as time.

The solution files can be loaded using the BAS-PRO plotting library, available at https://github.com/atmosalex/BAS-PRO_plotting. A copy of the BAS-PRO plotting library has also been bundled with this dataset (BAS-PRO_plotting-main.zip) to prevent potential compatibility issues arising from future updates to the online repository. It is recommend to follow the steps in the "Getting started" section of the plotting library README.md file, as this will result in plots of the solution, and will also convert the plaintext solution files into a single file in binary .cdf format which allows for faster loading.

The plaintext solution included in this dataset is made up of two sets of files which correspond to different grid resolutions:



Files ending in 'dyn.txt' are 'dynamic output' files, containing the sampled time evolution of f throughout the simulation period. The dynamic output grid is lower resolution than the original BAS-PRO simulation grid in order to save disk space. These files are useful for producing plots.

Files not ending in 'dyn.txt' are 'simulation grid' files, containing f at the final simulation epoch only, at the original simulation grid resolution. These files are useful for loading into BAS-PRO as an initial condition, or for plotting the final epoch at higher resolution.


The coordinate range of the 'simulation grid' ('dynamic output') files is as follows:



log10(μ/ (1MeV/G)) ranges from: 0.029384425 to 4.2519649 (0.17108176 to 4.1952860)

K ranges from: 0 to 5.729029 (0 to 5.729029) in units G0.5 RE

L ranges from: 1.13 to 4.0 (1.13 to 4.0)

time ranges from 1388534400 to 1517443200, given in terms of seconds passed since January 1, 1970 UTC, and this time range is from January 1, 2014 to February 1, 2018.


The following table gives a description of each file included:




axis_mu.txt
first dimension axis: a list of log10(μ/ (1MeV/G)) for the μ of each simulation grid point


axis_K.txt
second dimension axis: a list of K for each simulation grid point, with units G0.5 RE


axis_L.txt
third dimension axis: a list of L at each simulation grid point


axis_t.txt
time axis: a list of each simulation epoch, showing the history of timestepping


map_iK-aeq.txt
a 2D grid of equatorial pitch angle (degrees) corresponding to each L (rows) and K (columns) listed in the corresponding simulation axis files. A fill value of -1 is used to signify coordinates outside the trapping region.


iK-0001_2D_en.txt
a 2D grid of energy, with units of megaelectron volt, at each μ (rows) and L (columns) coordinate defined in the simulation axis files, at the K corresponding to the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K.txt file. 


iK-0001_2D_f.txt
a 2D grid of f, with units km-6 s3, at each μ (rows) and L (columns) coordinate defined in the simulation axis files, at the K corresponding to the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K.txt file.


iK-0001_axis_aeq.txt
a list of equatorial pitch angle (degrees) at each L in the axis_L.txt file, at the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K.txt file. A fill value of -1 is used to signify coordinates outside the trapping region.


...
...


axis_mu_dyn.txt
first dimension axis: a list of log10(μ/ (1MeV/G)) for the μ of each grid point in the dynamic output of the model


axis_K_dyn.txt
second dimension axis: a list of K for each grid point in the dynamic output of the model, with units G0.5 RE


axis_L_dyn.txt
third dimension axis: a list of L at each grid point in the dynamic output of the model


axis_t_dyn.txt
time axis: a list of each dynamic output epoch


map_iK-aeq_dyn.txt
a 2D grid of equatorial pitch angle (degrees) corresponding to each L (rows) and K (columns) listed in the corresponding axis files for the dynamic output. A fill value of -1 is used to signify coordinates outside the trapping region.


iK-0001_2D_en_dyn.txt
a 2D grid of energy, with units of megaelectron volt, at each μ (rows) and L (columns) coordinate defined in the dynamic output axis files, at the K corresponding to the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K_dyn.txt file. 


iK-0001_2D_f_dyn.txt
a 2D grid of f, with units km-6 s3, at each μ (rows) and L (columns) coordinate defined in the dynamic output axis files, at the K corresponding to the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K_dyn.txt file. The 2D grid is output for every timestep and appended to the file, so subsequent 2D grids correspond to subsequent timesteps at the same K.


iK-0001_axis_aeq_dyn.txt
a list of equatorial pitch angle (degrees) at each L in the axis_L_dyn.txt file, at the K index listed in the file name. For example, iK-0001... means the first K on the 3D model grid, corresponding to the first value of K listed in the axis_K_dyn.txt file. A fill value of -1 is used to signify coordinates outside the trapping region.


...
...


progress.txt
a file used by the BAS-PRO model to continue from partially complete simulations. It contains three values (one per line): epoch of the simulation start time; total simulation time elapsed (seconds); and a mode select value (1 for dynamic, 0 for steady state)


resume.config
a backup of the original configuration options used to execute the BAS-PRO simulation, used only by the model