##############################################
# Scriptfile: silicon_AR.lsf
#
# Copyright 2008, Lumerical Solutions, Inc.
##############################################

# plot transmission and reflection
f = getdata("T","f");
T = -transmission("T");
R = transmission("R");
plot(c/f*1e9,R,T,1-(R+T),"wavelength (nm)","Normalized power");
legend("Total reflection","Transmitted to substrate","Absorption in moth eyes");

# image |E|^2 at the silicon surface at several wavelengths
if (true) {#false not to plot
E2 = getelectric("T");
x = getdata("T","x");
y = getdata("T","y");
image(x*1e9,y*1e9,pinch(E2,4,1),"x (nm)","y (nm)",
        "|E|^2 at " + num2str(c/f(1)*1e9) + " nm");
p = round(length(f)/2);
image(x*1e9,y*1e9,pinch(E2,4,p),"x (nm)","y (nm)",
        "|E|^2 at " + num2str(c/f(p)*1e9) + " nm");
p = length(f);
image(x*1e9,y*1e9,pinch(E2,4,p),"x (nm)","y (nm)",
        "|E|^2 at " + num2str(c/f(p)*1e9) + " nm");
}

# compare with no structures
if (true) {
  load("silicon_motheye_noAR.fsp");
  abs_noAR=1-transmission("R");
load("silicon_motheye_AR.fsp");
abs_AR=1-transmission("R");
plot(c/f*1e9,abs_AR,abs_noAR,"wavelength (nm)","Total absorbed power","Absorption");
legend("Moth eye structure","Flat structure");
}