fig2.csv contains the data shown in Fig .2, namely d/\lambda, the linewidth of the out-of-plane uniform mode \upsilon_P, the linewidth of the in-plane uniform mode \upsilon_I, and the linewidth of the out-of-plane antiferromagnetic mode \upsilon_{AF}, all in units of \gamma.

fig3a.csv, fig3c.csv, and fig3e.csv each contain each contain the time in units of 1/gamma, then the occupation of the mode plotted in each subfigure, then the occupation of all other modes.
fig3b.csv, fig3d.csv, and fig3f.csv each contain a column giving the linewidth of each mode, and a column giving the occupation of each mode.

For figure four, each of the images has a separate file, e.g. fig4a.csv etc., which gives the data shown as a 31x31 array, corresponding to position in the atomic array. 

fig5.csv gives the total atom number N, followed by the time, in units of 1/gamma, taken to 30%$, 50%, 70%, and 90% respectively.

fig6a.csv gives the time in units of 1/\gamma, followed by the occupation of the out-of-plane mode L_P, the occupation of the in-plane mode L_I, and the total remaingin polulation L'.
fig6b.csv, fig6c.csv, and fig6e.csv each contain a 194x194 2D array which gives the log of the Intensity for a regularly spaced square grid where ky/kx and kz/kx range between -0.3 and 0.3.
fig6d.csv gives the total atom number N followed the the angle theta_max within wihc 99% of the integrated light intensity is emittedd.

fig7a.csv gives the time in units of 1/\gamm, and the remaining probability for no disorder, disorder l=0.02d, l=0.05d, and no disorder with in-plane excitation in that order.
fig7b.csv gives the central value of each histogram bin, the total density of occurences of modes within each bin, and the density of modes which which account for 30% of the initial occupation, built up from 5000 realizations.