mname = "signal::reflection";
# plot far field without bump
switchtolayout;
setnamed("DVD","addbump",0);
run;
E_near = getresult(mname,"E");
E2_near = E_near.E2;
x = E_near.x;
y = E_near.y;
image(x,y,E2_near,"","","flat surface");
setplot("x min",-2e-6); setplot("x max",2e-6); setplot("y min",-2e-6); setplot("y max",2e-6);
E2 = farfield3d(mname,1,200,200,1,1,1,1);
ux = farfieldux(mname,1,200,200,1);
uy = farfielduy(mname,1,200,200,1);
image(ux,uy,E2,"","","flat surface","polar");

# plot far field with bump
switchtolayout;
setnamed("DVD","addbump",1);
setnamed("DVD","bumpwidth",200e-9);
setnamed("DVD","bumplength",800e-9);
setnamed("DVD","bump_x",0);
setnamed("DVD","bump_y",0);
run;
E_near = getresult(mname,"E");
E2_near = E_near.E2;
x = E_near.x;
y = E_near.y;
image(x,y,E2_near,"","","with bump");
setplot("x min",-2e-6); setplot("x max",2e-6); setplot("y min",-2e-6); setplot("y max",2e-6);
E2 = farfield3d(mname,1,200,200,1,1,1,1);
ux = farfieldux(mname,1,200,200,1);
uy = farfielduy(mname,1,200,200,1);
image(ux,uy,E2,"","","with bump","polar");