Figs 2(a) and 2(b) from the manuscript are reproduced here. The following parameter values had to be changed. The original values are given in brackets. 
alpha = 0.55 (0.47)
S(0) = 78841214.13 (81800000)
tau_1 = 1.47e-9 (1.41e-9)