####################################################
# File: Thermal_emission.lsf
#
# Description:  This file will calculate thermal emission 
# from periodic photonic crystal and plot the results
# Copyright 2013 Lumerical Solutions
####################################################

switchtolayout;
setnamed("::model","a",2e-6);
save("thermal_emission2");
run;
te2=getresult("thermal_emission","te");

switchtolayout;
setnamed("::model","a",3e-6);
save("thermal_emission3");
run;
te3=getresult("thermal_emission","te");

plot(c/te2.f,te2.emissivity,te3.emissivity,"Wavelength","Emissivity");
legend("a=2 micron","a=3 micron");

plot(c/te2.f,te2.thermal_emission*1e-10,te3.thermal_emission*1e-10,te2.blackbody_spectrum*1e-10,"Wavelength","Power [W/m^2/micron]","T=1200K");
legend("a=2 micron","a=3 micron","Blackbody");