In this directory there is simulation code for reproducing Figure 3.
Usage: ./runs.sh parameters.dat
Following this command, the code is compiled with gfortran and ten simulation runs will start.
The simulation runs are ten repetitions of the same system and results should be averaged.

The runs.sh file overwrites some parameters in 'parameters.dat' to explore parameter space.
To reproduce Figure 1b, multiple values of a and multiple values of the x delay time should be explored.
  
Results will be printed in subdirectories as follows: 

sgr.[1-10].x.0.3.sr.0.16.flowcess.a.[0.0-32.0].xdecayt.[1.0D0-5000.0D0] 

