run;

f=getdata("Monitor1","f");
x=getdata("Monitor1","x");
y=getdata("Monitor1","y");
simf=f(1);

E2=getelectric("Monitor1");

# plot cross section of nearfield
image(x*1e9,y*1e9,E2,"x (nm)","y (nm)","Nearfield intensity (logscale)","logplot");

# plot y-slice of near field in middle of photoresist 
yslice=46;
?y(yslice)*1e9;
title="Near field intensity @"+num2str(c/simf*1e9)+" nm";
plot(x*1e9,pinch(E2,2,yslice),"x (nm)",title);