This is a guide to using CHEM1D to perform the ozone simulations presented in Siskind et al. paper for
a hypothetical supernova-induced soft Xray input

For compiling, it seems to require ifortran; I tried gfortran and it crashed

ifort -c -80 *.f
ifort *.o

then mv a.out chem1d_out_66
to run:
./chem1d_run_flare

It will ask for the name of your initial file ("restart", *.res). There are three files provided
nox_base.res, nox_x13.res and nox_x27.res

Then, the code will create a big unformatted OUT file with lots of extraneous reaction rates etc. 
You can use read_output.pro to read it, but easier is that
here this has been condensed into python readable IDL save files
A problem is that it doesn't
calculate pressure- you need the input temperature.
Can see fort.15, timegcm_t_m_p.txt for that
Anyway, fig11_plot.pro reads the 3 save files and shows how Figure 11 was made