This repository contains the data necessary to regenerate the plots in 
the paper preprint "A higher-order finite-element implementation of the exact
Landau Fokker-Planck collision operator for charged particle collisions
in a low density plasma" by M. R. Hardman, M. Abazorius, J. Omotani, 
M. Barnes, S. L. Newton, J. W. S. Cook, P. E. Farrell, and F. I. Parra.

To generate the pdfs containing the thematically organised plots in 

 /plots/

run 

 $ python plot_error_data.py 
 $ python plot_integration_error_data.py
 $ python plot_many_collisions.py

These plotting scripts use processed data from the original runs.
Requirements for a python environment can be found in the following file.
  
  requirements.txt

To extract the plots used in the paper, run

  $ ./pdfsplit.sh

Simulation and error data is stored within

 /data/ 

with the structure of the simulation data as follows.
For runs using the collision operator and numerical-error-correcting terms

 /data/runs/fokker-planck-relaxation-flux-bc-only*long/fokker-planck-relaxation-flux-bc-only*long.*

For runs using the collision operator without numerical error correction

 /data/runs/fokker-planck-relaxation-flux-bc-only-no-conserve*long/fokker-planck-relaxation-flux-bc-only-no-conserve*long.*

The *.toml files in /data/runs/ are the input files for the simulations.
The source code used for the simulations is 
available and documented at 

 https://github.com/mabarnes/moment_kinetics

with the commit 8bff6a5658e8e3012fb8c0bef81c01c02760fb22 being used for the simulations here.
