In this directory there is simulation code for reproducing Figure 1a.
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 1a, multiple values of a and multiple strain rates should be explored, as well as homogeneous and heterogeneous systems.
In addition to a=8 -  a=0,1,2,4 and 16.
In addition to strain rate=0.16  - strain rate=0.02-0.36
  
Results will be printed in subdirectories as follows: 

sgr.[1-10].x.0.3.sr.[0.04-0.32].steadystate.a.[0.0-32.0].[homo.]xdecayt.0.0D0 (or equivalently, any value of x decay time)

