/* Fundamental dimensions */
set step=um/20;
set Tpml = 8*step
set xsize=5.1*um
set ysize=xsize
set zsize=6.2*um;
set dipoleheight=0.2*um;
set dipolethickness=0.30*um;
set phi=0; /* phase of excitation */
set ky0=step/lambda;
set kx0=0; /* direction of plane wave */
set loadlength=2*step;
set loadThickness=dipolethickness
set loadheight=dipoleheight
set Zvertex=zsize-Tpml-loadthickness/2-1.5*um /* back end of load 1.5 um from PML */
set DipoleLengthUM = 1.40 ; /* optimal values from dipole2f.par */
set ZL = 80 /* characteristic impedance of load */
set ZRakeUM=.8
set sigmaL =loadlength/(ZL *loadThickness*loadheight)
set epsppL=sigmaL/(eps0*omega)
set n_load=sqrt(1+epsppL*epsppL)*cos(0.5*atan(epsppL))
set k_load=sqrt(1+epsppL*epsppL)*sin(0.5*atan(epsppL))
END
MATERIAL
DEFINE matname= Air type=dielectric n=1.0003 k=0
DEFINE matname= AirPML type=PML n=1.0003 k=0
DEFINE matname= PCL type=black
DEFINE matname= Gold type=metal n=Au_N(lambda) k=Au_K(lambda)
53