Data and plotting code for "Energetic electron scattering by kinetic Alfvén waves at strong magnetic field gradients of dipolarization front"

This publication includes the dataset and plotting code for reproducing results and figures of the manuscript "Energetic electron scattering by kinetic Alfvén waves at strong magnetic field gradients of dipolarization front", which was submitted to Physics of Plasma for potential publication. In this paper, we use test particle simulations to study energetic electron pitch angle and momentum scattering driven by kinetic Alfvén waves, at the magnetic field gradients associated with the dipolarization front, current sheet and plasma sheet fields in the magnetotail. 

The simulation data contains time history of 1000 test particle trajectories, momentum, energy, and KAW wave phases recorded every 0.001 second. Each simulation for different electron energies and pitch angles has a size of over 1 GB, so here we have only included one example of test particle simulation run associated with Figure 3. Only the final results of pitch angle diffusion coefficients are provided associated with Figure 4. A full set of time history results of all test particle simulation runs for Figure 4 will be available upon reasonable request ([e-mail]). 

Detailed descriptions:

 PoP_figure_1.m --- Matlab plotting code for reproducing magnetic field model and geometry for Figure 1.

 PoP_figure_2.m --- Matlab plotting code for reproducing resonant energy as a function of wave normal angle, electron pitch angle, and field model parameters in Figure 2.

 PoP_figure_3.m --- Reading test particle simulation data and plotting code for reproducing Figure 3. Please ensure proper directory setting to read test particle data files. 

 PoP_figure_4_5.m --- Reading saved bounce-averaged diffusion coefficient data and plotting Figure 4 and Figure 5.

 Daa_Energy_alpha_ad3xp15_e1.csv --- Electron diffusion coefficients driven by KAW scattering as a function of energy and pitch angle obtained from test particle simulations. KAW wave amplitude is set to 1 mV/m. Electrons are initiated from a position of x=150 km within the dipolarization front.

 Tau_Energy_alpha_xp15.csv --- Electron bounce periods as a function of energy and pitch angle based on numerical test particle simulations without KAW input.

 pKAW_1000p_dify_E100keV_PA5_Daa_E4_pop.zip --- Data files containing 1000 test particle trajectories, momentum and adiabatic invariant. This one single run is for electrons with an energy of 100 keV and a pitch angle of 5 degree, under influence of KAW wave electric field of 4 mV/m. This dataset is used to plot Figure 3. 

 pKAW_1000p_dify_init_Ppara_mue_E100_PA4_Daa_E4_pop.dat --- initial parallel momentum and first adiabatic invariant of 1000 test particles. Read this file for plotting Figure 3. 

 Bxyz.m --- calculate model magnetic field for given positions in x, y, and z

 FAC.m --- convert local x, y, z coordinates to field-aligned coordinates

 gradB.m --- calculate the gradient of magnetic field for given positions in x, y, and z

 E_to_pnorm_hsr.m --- calculate normalized momentum based on relativistic electron energy

 mysubplot.m --- plotting function to replace subplot routine of Matlab