In this directory there is simulation code for reproducing each panel of Figure 2.
Usage: ./runs.fig2[a-d].sh parameters.dat
Following this command, the code is compiled with gfortran and six simulation runs will start.
The simulation runs are for each of the six strain rates displayed in the panels of Figure 2.

The parameters in these files are for 'short' runs that should take only a few minutes to run.
For more accurate results, you can alter the 'parameters.dat' file to increase the number of atoms per streamline and the number of streamlines simulated. 
A longer simulation time can be requested by changing 'ntimecess' in the 'parameters.dat' file.

The runs.fig2[a-d].sh file overwrites some parameters in 'parameters.dat' to explore parameter space.
The runs.fig2[a-d].sh file also includes the facility for multiple repetitions of each run.
  
Results will be printed in subdirectories as follows: 

For Figure 2a:
sgr.short.fig2a.0.x.0.3.sr.[0.04-0.24].flowcess.a.0.0.homo.xdecayt.1.0D0 (or equivalently, any value of x decay time)

For Figure 2b:
sgr.short.fig2b.0.x.0.3.sr.[0.04-0.24].flowcess.a.8.0.xdecayt.1.0D0

For Figure 2c:
sgr.short.fig2c.0.x.0.3.sr.[0.04-0.24].flowcess.a.8.0.homo.xdecayt.100.0D0

For Figure 2d:
sgr.short.fig2d.0.x.0.3.sr.[0.04-0.24].flowcess.a.8.0.xdecayt.100.0D0


You can plot all the output files in this directory using:
gnuplot plotstress.gnu 
