fig1b.csv contains the data shown in Fig. 1 (b). The columns correspond 
to the detuning, and the linewidths of P_ya, P_ys, P_xa, and P_sx, in order.

fig3a_r.csv and fig3a_i.csv form a 2D array of the real and imaginary parts
of the incident field shown in Fig. 3 (a), respectively. fig3b_r.csv and 
fig3b_i.csv form a 2D array of the real and imaginary parts of the 
scatterd field shown in Fig. 3 (b), respectively. The total field, shown in
Fig. 3 (c) is the sum of the fields in (a) and (b).

fig4a.csv and fig4c.csv have columns corresponding to the time, P_in, P_Delta, 
and P_0 shown in Fig. 4 (a) and (c), respectively. fig4b.csv and fig4d.csv have
columns corresponding to the time, E_0, E_x, and E_y shown in Fig. 4 (b) and (d).
respectively. 

fig5a.csv has columns corresponding to the time, P_f, and P_b shown in Fig. 5 (a).
fig5b.csv has columns corresponding to the time, E_f, and E_b shown in Fig. 5 (b).

fig6a.csv has columns corresponding to the time, L_xa, and L_xs shown in
Fig. 6 (a). fig6b.csv shows the data markers plotted in Fig. 6 (b), with columns
corresponding to the detuning, max(L_xs), and Omega_osc, respectively.

fig7a.csv has columns corresponding to time t, incident power P_in, transmitted
power P_0, and reflected power P_r.

fig7b.csv has columns corresponding to time t, energy of the quadrupole mode E_q,
and energy of the dipole mode E_d.