BAS-PRO Model Solution for "Modeling Field Line Curvature Scattering Loss of 1 to 10MeV Protons during Geomagnetic Storms"

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文件ModelingFLC_BASPRO_solution.zip为BAS-PRO模型求解结果的归档压缩包。解压该压缩包后，求解结果将分散存储于多个纯文本文件中。该求解结果为质子相空间密度与质子静质量三次方的乘积构成的网格数据，记为$f$，单位为km⁻⁶·s³，$f$由第一、第二、第三绝热不变量μ、K、L以及时间共同定义。 可通过BAS-PRO绘图库加载该求解结果文件，该库的开源地址为https://github.com/atmosalex/BAS-PRO_plotting。本数据集同时附带了该绘图库的离线副本（BAS-PRO_plotting-main.zip），以避免因线上仓库后续更新引发的兼容性问题。建议按照绘图库README.md文件的"Getting started"（快速入门）章节操作，既可直接生成求解结果的可视化图像，还可将纯文本格式的求解结果文件转换为单一二进制.cdf格式文件，从而实现更快的数据加载速度。 本数据集包含的纯文本求解结果分为两组，分别对应两种不同的网格分辨率： 以"dyn.txt"结尾的文件为「动态输出文件」，存储了整个模拟周期内$f$的采样时间演化数据。为节省磁盘空间，动态输出的网格分辨率低于原始BAS-PRO模拟网格的分辨率，此类文件可用于生成可视化图像。 不以"dyn.txt"结尾的文件为「模拟网格文件」，仅存储最终模拟时刻的$f$数据，且采用原始模拟网格的分辨率。此类文件可用于作为初始条件导入BAS-PRO模型，或以更高分辨率绘制最终时刻的可视化结果。 「模拟网格文件」（「动态输出文件」）的坐标范围如下： - log₁₀(μ/(1MeV/G)) 的取值范围为：0.029384425 至 4.2519649（动态输出对应范围为0.17108176 至 4.1952860） - K的取值范围为：0 至 5.729029（动态输出与模拟网格范围一致），单位为G^0.5·R_E（R_E为地球半径） - L的取值范围为：1.13 至 4.0（动态输出与模拟网格范围一致） - 时间戳范围为1388534400至1517443200（以协调世界时1970年1月1日起的秒数计），对应实际时间范围为2014年1月1日至2018年2月1日。 下表对本数据集包含的各文件进行说明： 1. axis_mu.txt：第一维度轴文件，存储模拟网格各节点对应的log₁₀(μ/(1MeV/G))数值列表 2. axis_K.txt：第二维度轴文件，存储模拟网格各节点对应的K值列表，单位为G^0.5·R_E 3. axis_L.txt：第三维度轴文件，存储模拟网格各节点对应的L值列表 4. axis_t.txt：时间轴文件，存储各模拟时刻的数值列表，用于体现模拟的时间步进历程 5. map_iK-aeq.txt：二维网格文件，存储对应模拟轴文件中各L（行）与K（列）对应的赤道投掷角（单位：度）。若坐标处于捕获区域外，则以填充值-1标记 6. iK-0001_2D_en.txt：二维网格文件，存储模拟轴文件中定义的各μ（行）与L（列）坐标对应的能量数值，单位为兆电子伏特，其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K.txt文件中列出的第一个K值相对应 7. iK-0001_2D_f.txt：二维网格文件，存储模拟轴文件中定义的各μ（行）与L（列）坐标对应的f数值，单位为km⁻⁶·s³，其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K.txt文件中列出的第一个K值相对应 8. iK-0001_axis_aeq.txt：数值列表文件，存储axis_L.txt文件中各L对应的赤道投掷角（单位：度），其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K.txt文件中列出的第一个K值相对应。若坐标处于捕获区域外，则以填充值-1标记 9. （后续iK-xxxx系列文件命名规则与上述一致，此处省略） 针对动态输出组的文件说明如下： 1. axis_mu_dyn.txt：第一维度轴文件，存储模型动态输出中各网格节点对应的log₁₀(μ/(1MeV/G))数值列表 2. axis_K_dyn.txt：第二维度轴文件，存储模型动态输出中各网格节点对应的K值列表，单位为G^0.5·R_E 3. axis_L_dyn.txt：第三维度轴文件，存储模型动态输出中各网格节点对应的L值列表 4. axis_t_dyn.txt：时间轴文件，存储各动态输出时刻的数值列表 5. map_iK-aeq_dyn.txt：二维网格文件，存储动态输出对应轴文件中各L（行）与K（列）对应的赤道投掷角（单位：度）。若坐标处于捕获区域外，则以填充值-1标记 6. iK-0001_2D_en_dyn.txt：二维网格文件，存储动态输出轴文件中定义的各μ（行）与L（列）坐标对应的能量数值，单位为兆电子伏特，其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K_dyn.txt文件中列出的第一个K值相对应 7. iK-0001_2D_f_dyn.txt：二维网格文件，存储动态输出轴文件中定义的各μ（行）与L（列）坐标对应的f数值，单位为km⁻⁶·s³，其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K_dyn.txt文件中列出的第一个K值相对应。该二维网格会随每一步时间步进输出并追加至文件中，因此后续的二维网格依次对应同一K参数下的不同模拟时刻 8. iK-0001_axis_aeq_dyn.txt：数值列表文件，存储axis_L_dyn.txt文件中各L对应的赤道投掷角（单位：度），其K参数对应文件名中指定的K索引。例如，iK-0001...代表三维模型网格中的第一个K值，与axis_K_dyn.txt文件中列出的第一个K值相对应。若坐标处于捕获区域外，则以填充值-1标记 9. （后续动态输出组的iK-xxxx系列文件命名规则与上述一致，此处省略） 10. progress.txt：BAS-PRO模型用于从部分完成的模拟任务中恢复运行的配置文件。文件按行存储三个数值：模拟开始时刻的时间戳、已完成的模拟总时长（单位：秒），以及模式选择值（1代表动态模式，0代表稳态模式） 11. resume.config：本次BAS-PRO模拟所使用的原始配置选项的备份文件，仅用于模型运行