fig2.csv: Contains the data plotted in Fig. 2. THe first column gives L, the occupation of the magnetic mode, the second column gives L', the occupation of all other modes, the third column gives the complex phase of r and the fourth column gives the absolute vvalue of r. Data is for a regular grid of detuning from -7.4 to -6.2 in steps of 0.012

fig3a.csv: Gives the 2D array of complex total field values with y varying from -70/(2*pi) to 70/(2*pi) in steps of 0.6/(2*pi) and x varying from -25/(2*pi) to 15/(2*pi) in steps of 0.0 6/(2*pi) for the magnetic mirror shown in Fig. 3(a)

fig3b.csv: Gives the 2D array of complex total field values with y varying from -70/(2*pi) to 70/(2*pi) in steps of 0.6/(2*pi) and x varying from -25/(2*pi) to 15/(2*pi) in steps of 0.0 6/(2*pi) for the electric mirror shown in Fig. 3(b)

fig4.csv: Column one is the laser detuning, and column two gives the complex transmitted field shown in Fig. 4

fig5a.csv: Gives the 2D array of complex total field values with y varying from -70/(2*pi) to 70/(2*pi) in steps of 0.3/(2*pi) and x varying from -25/(2*pi) to 125/(2*pi) in steps of 0.6/(2*pi)
fig5b.csv: Gives the complex total field shown in Fig. 5b for y varying from -40/(2*pi) to 40/(2*pi) in steps of 0.05/(2*pi)
fig5c.csv: Gives the complex total field shown in Fig. 5c for x varying from 0.1/(2*pi) to 120/(2*pi) in steps of 0.05/(2*pi)

fig6.csv: Gives the 2D array of complex transmitted field values with x varying from -40/(2*pi) to 100/(2*pi) in steps of 0.4/(2*pi) and z varying from -70/(2*pi) to 70/(2*pi) in steps of 0.4/(2*pi)

fig7ab.csv: Gives the 2D array of complex transmitted field values with  y varying from -50/(2*pi) to 500/(2*pi) in steps of 0.6/(2*pi) and z varying from -50/(2*pi) to 50/(2*pi) in steps of 0.6/(2*pi) for the case shown in Fig. 7(a) and (b)
fig7cd.csv: Gives the 2D array of complex transmitted field values with  y varying from -50/(2*pi) to 500/(2*pi) in steps of 0.6/(2*pi) and z varying from -50/(2*pi) to 50/(2*pi) in steps of 0.6/(2*pi) for the case shown in Fig. 7(c) and (d)
fig7ef.csv: Gives the 2D array of complex transmitted field values with  y varying from -50/(2*pi) to 500/(2*pi) in steps of 0.6/(2*pi) and z varying from -50/(2*pi) to 50/(2*pi) in steps of 0.6/(2*pi) for the case shown in Fig. 7(e) and (f)