# get monitor data
x=getdata("field","x");
y=getdata("field","y");
f=getdata("field","f");
Ez=getdata("field","Ez");

# plot actual source angle vs wavelength
theta=getsourceangle("source1",f);
plot(f/1e12,theta,"frequency (THz)","actual source angle","Nominal source angle is 30 deg");

#plot field profile at min and max wavelength
fi=1;           # first frequency
Ez_plot=pinch(Ez,4,fi);  
Ez_plot=real(Ez_plot);
image(x*1e6,y*1e6,Ez_plot,"x (um)","y (um)","Ez at "+num2str(f(fi)/1e12)+"THz");

fi=length(f);   # last frequency
Ez_plot=pinch(Ez,4,fi);  
Ez_plot=real(Ez_plot);
image(x*1e6,y*1e6,Ez_plot,"x (um)","y (um)","Ez at "+num2str(f(fi)/1e12)+"THz");