# scriptfile: grating_silicon_nitride.lsf
# 
# description: This script performs a transmission
# 	and reflection analysis of a sub-wavelength
#	grating made in silicon nitride on epoxy,
#	simulated with grating_silicon_nitride.fsp and
#	described in Cunnignham et al,. "A plastic 
#	colorimetric resonant optical biosensor for 
#	multiparallel detection of label-free biochemical 
#	interactions", Sensors and Actuators, B, 85, 
#	219-226 (2002). 
#
# copyright 2006, Lumerical Solutions, Inc.

# calculate the transmission and reflection
# note the negative sign when calculating the
# reflection. This is because we want a positive
# reflection when power propagates in the 
# negative y direction for this problem.
T = transmission("transmission");
R = -transmission("reflection");
f = getdata("transmission","f");

# plot the results vs wavelength in microns
plot(c/f*1e6,R,T,R+T,"wavelength (microns)","Normalized power");
legend("R","T","R+T");

# plot only the reflection
plot(c/f*1e6,R,"wavelength (microns)","Reflection");