BAS-PRO solutions J81, SA19 and S16

Overview

This set of directories contains three solutions from the BAS-PRO proton radiation belt model. These solutions are saved as .cdf files, making available unidirectional, differential proton flux j as a function of time, first invariant mu, second invariant K and third invariant L or phi. The value of the distribution function f is also available (see below format information).

The time period modelled is from 1st March 2014, to 1st February 2018.

The three solutions correspond to different radial diffusion coefficients ('DLL') like so:

    - solution 20220203_220342 corresponds to the DLL from Jentsch, 1981 (Equation 18, Selesnick et al., 2007)
    - solution 20220203_220420 corresponds to the post-1st Jan 2015 DLL from Equation 5, Selesnick & Albert, 2019 (SA19)
    - solution 20220203_220436 corresponds to the DLL from Equation 12, Selesnick et al., 2016 (S16)

There are two .cdf files for each solution: one file is a high resolution static solution, corresponding to 1st February 2018, which is the final output of each model run; the second (larger) file is a lower resolution dynamic (time-dependent) solution, corresponding to the ~four year period above. Both solutions were generated from a high resolution model run, but the dynamic solution was output at a lower resolution to save disk space.

These solutions can be considered updates to the work shown in Lozinski et al., 2021 (https://doi.org/10.1029/2021JA029777)

A set of Python scripts to create plots of the data are available at https://github.com/AlexisNaN/BAS-PRO_plotting. Example plots of each dynamic solution are included in the attached data, produced using this set of Python scripts.

Format of the dynamic solution cdf files (i.e. 20220203_220420_solution_dyn.cdf):

 axis_mu
  shape: (482,)
  units: log10(mu/ (1 MeV/G) )

 axis_K
  shape: (35,)
  units: G^0.5 RE

 axis_L
  shape: (34,)

 axis_t
  shape: (206,)
  units: seconds since Jan 01 1970 (UTC)

 axis_t_date
  shape: (206,)
  units: datetime CDF object

 map_KL-aeq
  shape: (34, 35)
  units: degrees

 axis_phi
  shape: (34,)
  units: T m2

 f
  shape: (206, 482, 35, 34)
  units: km-6 s3
  desc: distribution function: relativistic phase space density multiplied by proton rest mass cubed

 energy
  shape: (1, 482, 35, 34)
  units: MeV
  desc: energy at each data coordinate, assumed constant in time (ignoring secular variation)

 j
  shape: (206, 482, 35, 34)
  units: cm-2 s-1 str-1 MeV-1
  desc: unidirectional differential proton flux

Format of the static solution cdf files (i.e. 20220203_220420_solution.cdf):

 axis_mu
  shape: (700,)
  units: log10(mu/ (1 MeV/G) )

 axis_K
  shape: (35,)
  units: G^0.5 RE

 axis_L
  shape: (171,)

 axis_t
  shape: (1,)
  units: seconds since Jan 01 1970 (UTC)

 axis_t_date
  shape: (1,)
  units: datetime CDF object

 map_KL-aeq
  shape: (171, 35)
  units: degrees

 axis_phi
  shape: (171,)
  units: T m2

 f
  shape: (1, 700, 35, 171)
  units: km-6 s3
  desc: distribution function: relativistic phase space density multiplied by proton rest mass cubed

 energy
  shape: (1, 700, 35, 171)
  units: MeV
  desc: energy at each data coordinate, assumed constant in time (ignoring secular variation)

 j
  shape: (1, 700, 35, 171)
  units: cm-2 s-1 str-1 MeV-1
  desc: unidirectional differential proton flux