fini
/cle

/fil,hf_prdc

/vie,1,2,1,1
/vup,1,z
/pnu,mat,1
/num,1
 
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  


C**************************************************************
C*** PARAMETERS (DIMENSIONS IN METERS)
C**************************************************************
frqncy=10e9		! FREQUENCY (Hz)
phi=90			! INCOMING PLANEWAVE (ANGLE FROM X TO Y AXIS)
thta=45			! INCOMING PLANEWAVE (ANGLE FROM Z AXIS)

efx=1			! INCOMING PLANEWAVE POLARITY (V/m)
efy=0
efz=0

lambda=3e8/frqncy	! FREE SPACE WAVELENGTH (m)

l=lambda/2		! DIPOLE LENGTH
w=l/5			! DIPOLE WIDTH
prd=1.25*l			! SIZE OF PERIODIC CELL

t_sbstrt=w		! SUBSTRATE THICKNESS

n_cellsx=20		! # OF PERIODIC CELLS IN X DIRECTION
n_cellsy=20		! # OF PERIODIC CELLS IN Y DIRECTION

h_eqv_src=prd		! HEIGHT OF EQUIVALENT SOURCE ABOVE SUBSTRATE
h_port=1.5*prd		! HEIGHT OF PORT ABOVE SUBSTRATE
h_dmn=2.0*prd		! HEIGHT OF PML ABOVE SUBSTRATE
d_PML=lambda/4		! PML DEPTH

esz_xy=lambda/20	! ELEMENT SIZE IN XY PLANE
dvz01=3			! # OF DIVISIONS THRU SUBSTRATE
dvz02=15		! # OF DIVISIONS TO EQV SRC
dvz03=5			! # OF DIVISIONS TO PORT
dvz04=5			! # OF DIVISIONS TO PML
dvz05=4			! # OF DIVISIONS THROUGH PML


C**************************************************************
C*** MATERIAL PROPERTIES
C**************************************************************
!per_2r=7.5		! SUBSTRATE RELATIVE PERMITTIVITY (REAL PART)
!per_2i=0		! SUBSTRATE RELATIVE PERMITTIVITY (IMAG PART)

per_2=7.5		! SUBSTRATE RELATIVE PERMITTIVITY
mur_2=1			! SUBSTRATE RELATIVE PERMEABILITY
rsv_2=4.5e-2		! SUBSTRATE ELECTRICAL RESITIVITY (ohm-m)

!lsst_2=per_2i/per_2r	! SUBSTRATE LOSS TANGENT


C**************************************************************
C*** CREATE PLANAR GEOMETRY AND MESH
C**************************************************************
/prep7
rect,,l/2,,w/2
rect,,prd/2,,prd/2
aovl,all

et,200,200,7
esiz,esz_xy
ames,all


C**************************************************************
C*** DEFINE ATTRIBUTES
C**************************************************************
et,1,120	! AIR ATTRIBUTES
mp,perx,1,1
mp,murx,1,1

et,2,120	! SUBSTRATE ATTRIBUTES
mp,perx,2,per_2
mp,murx,2,mur_2
mp,rsvx,2,rsv_2

!mp,lsst,2,lsst_2

et,11,120,,,,1	! PML ATTRIBUTES
mp,perx,11,1
mp,murx,11,1


C**************************************************************
C*** EXTRUSIONS
C**************************************************************
exto,esiz,dvz01
asel,s,loc,z
mat,2 $real,2 $type,2
vext,all,,,0,0,-t_sbstrt

exto,esiz,dvz02
asel,s,loc,z
mat,1 $real,1 $type,1
vext,all,,,0,0,h_eqv_src

exto,esiz,dvz03
asel,s,loc,z,h_eqv_src
mat,1 $real,1 $type,1
vext,all,,,0,0,h_port-h_eqv_src

exto,esiz,dvz04
asel,s,loc,z,h_port
mat,1 $real,1 $type,1
vext,all,,,0,0,h_dmn-h_port

exto,esiz,dvz05
asel,s,loc,z,h_dmn
mat,11 $real,11 $type,11
vext,all,,,0,0,d_PML


C**************************************************************
C*** REFLECTION
C**************************************************************
alls
vsym,x,all
vsym,y,all
numm,node,1e-3*w,1e-3*w
numm,kp,1e-3*w,1e-3*w
acle,all
etde,200


C**************************************************************
C*** DEFINE COMPONENTS FOR BC APPLICATION
C**************************************************************
asel,s,loc,z,0
asel,r,loc,x,-l/2,l/2
asel,r,loc,y,-w/2,w/2
cm,dipole_a,area	! USER TO SET AX=0 (PEC) ON DIPOLE

asel,s,loc,z,-t_sbstrt
cm,grnd_a,area		! USER TO SET AX=0 (PEC) ON GROUND PLANE

asel,s,loc,z,h_dmn+d_PML
cm,PML_a,area		! USER TO SET AX=0 (PEC) ON TOP OF PML

vsel,s,loc,z,0,h_eqv_src
alls,belo,volu
nsel,r,loc,z,h_eqv_src
esln,r
cm,eqvsrc_n,node
cm,eqvsrc_e,elem
cmgr,eqv_src,eqvsrc_n,eqvsrc_e	! USER TO SET MXWF FLAGS ON THESE NODES

asel,s,loc,z,h_port
cm,port_a,area		! USER TO DEFINE PORT ON THESE AREAS

asel,none
cmse,all
alls,belo,area
cmse,a,eqv_src
cmplot

/eof			! USER TO DO ALL THAT FOLLOWS



C**************************************************************
C*** BOUNDARY CONDITIONS: PEC
C**************************************************************
alls
cmse,all
da,PML_a,ax
da,grnd_a,ax
da,dipole_a,ax	! COMMENT OUT THIS LINE TO SUPPRESS DIPOLE PEC


C**************************************************************
C*** BOUNDARY CONDITIONS: EQUIV SOURCE SURFACE
C**************************************************************
alls
cmse,s,eqv_src
sf,all,mxwf


C**************************************************************
C*** BOUNDARY CONDITIONS: PERIODICITY
C**************************************************************
csys

nsel,s,loc,x,-prd/2
nsel,a,loc,x,prd/2
cpcyc,ax,,,prd,0,0,1

nsel,s,loc,y,-prd/2
nsel,a,loc,y,prd/2
cpcyc,ax,,,0,prd,0,1


C**************************************************************
C*** EXCITATION: PLANE WAVE PORT
C**************************************************************
alls
cmse,all
bf,port_a,port,1
hfport,1,plan,,,soft,1,0,0,phi,thta


fini


C**************************************************************
C*** SOLVE
C**************************************************************
/solu
eqslv,sparse
antype,harmic
harfrq,frqncy
eplo
save
solve
fini


C**************************************************************
C*** POSTPROCESS
C**************************************************************
/post1
set

plve,ef,,,,vect,,on

hfarray,n_cellsx,n_cellsy,prd,prd   
plhffar,field,ef,phi,phi,0,-90,0,90,0.1,0,x 

fssparm,1,1
hfpower,1,1

/eof