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.