## Script to calculate incoherent image field from two simulations
## with illumination sources having orthogonal polarization 
##
## Copyright Lumerical Solutions Inc.
##
####

sn="litho_cross";
#sn="litho_4sq";

setnamed("image_field","M",1);
setnamed("image_field","plot_fields",1);

runsweep(sn); #run the two orthogonal simulations


TransE2_inc=getsweepdata(sn,"TransE2_image"); #get the summed incoherent data
x_image=getsweepdata(sn,"x_image");
y_image=getsweepdata(sn,"y_image");
E2norm_inc = TransE2_inc/max(TransE2_inc); #normalize the data

image(x_image,y_image,E2norm_inc,"x (microns)","y (microns)","Image Intensity");
#thresh=0.4;
#image(x_image,y_image,E2norm_inc>thresh,"x (microns)","y (microns)","Image Intensity");