Suprathermal Proton Spectra at Interplanetary Shocks in 3D Hybrid Simulations

This repository contains the source code and input files required to reproduce the simulation runs used in Young, et al. (2020) "Suprathermal Proton Spectra at Interplanetary Shocks in 3-D Hybrid Simulations", published in the Astrophysical Journal.

The source code is in src.tar. The remaining files are input files. Each simulation run requires six input files: dat.rol, dat.spe, dat.cod, dat.list, dat.dia, and dat.bdist. The name of each input file encodes the shock-normal angle in degrees as "tbn_##" and the upstream flow speed as "v1_#".

Instructions for reproducing the simulations:

1. Modify the makefile for your system.
2. Compile the source code with 

   make ych.out

    

3. Place each set of six input files with a common "tbn_##-v1_#" prefix into a separate directory.
4. Remove the "tbn_##-v1_#" from each input file name.
5. Run the executable in each directory. The exact procedure will depend on your batch system. The simulation runs used in Young et al. (2020) ran on 128 processors, on a CRAY XE6m-200.