hm   Harmonic Number Ratio. The larger the 
     absolute value, the more accurate, 
     but the greater the computational 
     load. Normally set to 4.0 to 7.0.
     >0 : Harmonics number is an odd value,
          rounding up |hm|*(wdx+2*wb)/Lam 
      or |hm|*(wdy+2*wb)/Lam.
     <0: Harmonics number is an even value
         (i.e., cited odd value -1).
trc  Truncation factor (>=0) that works only 
     for 3-D problems.
     If =0, it is treated as no truncation.
     The smaller the coefficient, the more 
     accurate, but the greater the 
     computational load. Normally set to 1.0 
     to 3.0.
wb   Absorbing boundary width (um).
     =0 is for no absorbing boundaries.
     The larger the value, the more non-
     reflective the boundaries are, but
     the amount of calculation increases.
kfl  =0: Drawing without absorbing 
         boundaries.
     =1: Drawing with absorbing boundaries.
kot  Distributions such as intensity, 
     absorption, and refractive index are 
     output in a maximum of (kot+5) digits.
ity  Definition of intensity distribution.
    =0 : magnitude of Poynting Vector,
    =1 : electric & magnetic fielf intensity, 
    =2 : electric fielf intensity,
    =3 : magnetic fielf intensity.
--------------------------------------------
wdx  Analysis width in the x direction (um).
     wdx=0 becomes a 2-dimensional problem.
     The center of the width is the 
     positional basis for the light source 
     and structures. The larger the width, 
     the larger the number of harmonics 
     required.
wdx  Analysis width in the y direction (um).
     wdy=0 becomes a 2-dimensional problem.
     The center of the width is the 
     positional basis for the light source 
     and structures. The larger the width, 
     the larger the number of harmonics 
     required.
dxy  Grid interval in x/y direction (um). 
     The actual interval is optimized to be 
     close to that and displayed in wsr.out.
     The value should be less than 1/10 of
     the wavelength.
dz   Grid interval in-z direction (um).
---------------------------------------------
Lam  Wavelength (um).
th   Azimuth angle of incident light with 
     z-axis (deg).
fi   Argument angle of incident light in 
     xy-plane (deg).
gm   Polarization direction of light 
     source (deg). Angle that the electric 
     field vector makes with the x-axis
--------------------------------------------
alx  Light source spread in x-direction (um).
      =1.0 : uniform distribution.
      =0.0 : Minimum x-width.
      =0-1 : full width half maximum 
             in x direction = wdx*alx
aly  Light source spread in y-direction (um).
      =1.0 : uniform distribution.
      =0.0 : Minimum y-width.
      =0-1 : full width half maximum 
             in y direction = wdx*alx.
sx0  Shift length of light source center in 
     x-direction (um).
sy0  Shift length of light source center in 
     y-direction (um).
--------------------------------------------
stx  Shift length of overall structure 
     center in x-direction (um).
     Not applicable for light source 
     position.
sty  Shift length of overall structure 
     center in y-direction (um).
     Not applicable for light source 
     position.
csx  Cross sectional position of graphics
     in x-direction (um).
csy  Cross sectional position of graphics
     in y-direction (um).

--------- km designation field -------------
           (for optical materials)
     The first 4 digits are serial line 
     numbers, up to 200 lines can be input.
Name Material name (within 8 digits)
     SiO2,Ag,Al,Au,Be,Cr,Cu,Ni,Pd,Pt,Ti,W 
     have internal data. For others, by 
     entering the wavelength, refractive 
     index, and extinction coefficient in 
     the file of nk.dat as external data, 
     the refractive index and extinction 
     coefficient are automatically 
     interpolated. If no data exists in 
     nk.dat, the values defined by the 
     right-side parameters of 'an' are given 
     priority. The nk.dat should be created 
     by each user and stored in the same 
     folder as wsr.exe.
ko   Whether to output detected light amount 
     to wsb1.out or not.
     =0 : not output, =1 : output.
an   Refractive index.
ab   Abbe number, if =0, no dispersion 
     (fixed to refractive index).
ak   Extinction coefficient.

--------- kr designation field -------------
         (for restricting shapes)
     The first 4 digits are serial line 
     numbers, up to 1000 lines can be input.
kd   How to input shape data of structures.
     =0: by internal definition.
     =1,-1: by external data using sub.dat.
         Applied to all except for wx and 
         wy, sub.dat can be input up to 400 
         types (up to 1000 lines for each 
         type).
kt   Selection of shape type.
     When kd=1, kt=Pattern No. in sub.dat.
     kd=-1 shows an inverted shape for kd=1.
     When kd=0,-kt shows an inverted shape
          for kt.
      kt=0 No area restriction.
        =1 Restricted by rectangular areas 
           of width wx*wy centered on a 
           square grid position of period 
           px*py.
        =2 Restricted by elliptic shape of 
           width wx*wy and elliptic index xp
           centered on a square grid 
           position of period px*py,
           where xp=-2.0 to -1.0 for star, 
                   =-1.0 for diamond, 
                   = 0.0 for ellipse,
                   > 0.0 for square.
        =3 Restricted by hexagons shape (top
           /bottom vertex angles) of width 
           wx*wy centered on a square grid 
           position of period px*py.
        =4 Restricted by hexagon shape (left
           /right vertex angles) of width 
           wx*wy centered on a square grid 
           position of period px*py.
        =5 Restricted by diamond shape of 
           width wx*wy centered on a square 
           grid position of period px*py.
        =6 Restricted by a right-angled 
           triangular shape (diagonal fist 
           quadrant) of width wx*wy centered 
           on a square grid position of 
           period px*py.
        =7 Restricted by a right-angled 
           triangular shape (diagonal second 
           quadrant) of width wx*wy centered 
           on a square grid position of 
           period px*py.
        =8 Restricted by a right-angled 
           triangular shape (diagonal third 
           quadrant) of width wx*wy centered 
           on a square grid position of 
           period px*py.
        =9 Restricted by a right-angled 
           triangular shape (diagonal fourth 
           quadrant) of width wx*wy centered 
           on a square grid position of 
           period px*py.
ps   Rotation angle of all structures around 
     the region center (deg).
px   Structure period in x-direction (um).
     When =0, it is an isolated pattern.
py   Structure period in y-direction (um).
     When =0, it is an isolated pattern.
wx   Structure width in y-direction (um). 
wy   Structure width in y-direction (um). 
sx   Shift length of the structure center 
     in x-direction (um).
sy   Shift length of the structure center 
     in y-direction (um).
xp   Elliptic exponent index for kt=2.
     Lattice duty ratio for kt=10 to 17. 

--------- kf designation field -------------
         (for foreground structures) 
     The first 4 digits are serial line 
     numbers, up to 9999 lines can be input.
km   Construction material number referred in 
     km designation field.
     km=0 means vacuum (n=1.0).
kr   Restriction shape number referred in 
     kr designation field.
     kr=0 means restriction free.
kd   How to input shape data of structures.
     =0: by internal definition.
     =1,-1: by external data using sub.dat.
         Applied to all except for wx and 
         wy, sub.dat can be input up to 400 
         types (up to 1000 lines for each 
         type).
kt   Selection of shape type.
     When kd=1, kt=Pattern No. in sub.dat.
     kd=-1 shows an inverted shape for kd=1.
     When kd=0,-kt shows an inverted shape
          for kt.
      kt=0 No area definition.
        =1 Rectangular areas of width wx*wy 
           centered on a square grid position
           of period px*py.
        =2 Elliptic shape of width wx*wy and
           elliptic index xp centered on a 
           square grid position of period 
           px*py, where xp=-2.0 to -1.0 for
           star, 
                   =-1.0 for diamond, 
                   = 0.0 for ellipse,
                   > 0.0 for square.
        =3 Hexagons shape (top/bottom vertex
           angles) of width wx*wy centered 
           on a square grid position of 
           period px*py.
        =4 Hexagon shape (left/right vertex 
           angles) of width wx*wy centered 
           on a square grid position of 
           period px*py.
        =5 Diamond shape of width wx*wy 
           centered on a square grid position
           of period px*py.
        =6 Right-angled triangular shape 
           (diagonal 1st quadrant) of 
           width wx*wy centered on a square 
           grid position of period px*py.
        =7 Right-angled triangular shape 
           (diagonal 2nd quadrant) of 
           width wx*wy centered on a square 
           grid position of period px*py.
        =8 Right-angled triangular shape 
           (diagonal 3rd quadrant) of width
           wx*wy centered on a square grid 
           position of period px*py.
        =9 Right-angled triangular shape 
           (diagonal 4th quadrant) of 
           width wx*wy centered on a square 
           grid position of period px*py.
       =10 Linear lattice of period wx, angle 
           wy, duty ratio xp, starting point 
           xq included in each square grid 
           of period px*py.
       =11 Concentrically elliptic lattice 
           of period wx, angle wy, duty 
           ratio xp, starting point xq 
           included in each square px*py of 
           period px*py.
       =12 Concentrically dodecagonal 
           lattice of period wx, angle wy, 
           duty ratio xp, starting point xq 
           included in each square grid of 
           period px*py.
       =13 15-degrees-rotated lattice for 
           kt=12.
       =14 Concentrically 18-corner 
           polygonal lattice of period wx, 
           angle wy, duty ratio xp, starting 
           point xq included in each square 
           grid of period px*py.
       =15 10-degrees-rotated lattice for 
           kt=14.
       =16 Concentrically hexagonal lattice 
           of period wx, angle wy, duty 
           ratio xp, starting point xq 
           included in each square grid of 
           period px*py.
       =17 30-degrees-rotated lattice for 
           kt=16.
ps   Rotation angle of all structures around 
     the region center (deg).
px   Structure period in x-direction (um).
     When =0, it is an isolated pattern.
py   Structure period in y-direction (um).
     When =0, it is an isolated pattern.
wx   Structure width in y-direction (um). 
wy   Structure width in y-direction (um). 
sx   Shift length of the structure center 
     in x-direction (um).
sy   Shift length of the structure center 
     in y-direction (um).
xp   Elliptic exponent index for kt=2.
     Lattice duty ratio for kt=10 to 17. 
xq   Starting point of lattice for kt=10 
     to 17.

--------- Background layers -----------
     Up to 10000 lines can be input as far 
     as the last line or the line starting 
     from "c" appears. Optical constants 
     above the top layer or below the bottom 
     layer is the same ones as the top or the 
     bottom layer, respectively, and then no 
     boundary reflections from there.
kl   =1: light source position.
         If all of kl are 0, the top surface 
         of the first layer for cos(th)>0, or 
         the lower surface of the last layer 
         for cos(th)<0 is the light source 
         position.
km   Construction material number referred in 
     km designation field.
     km=0 means vacuum (n=1.0).
kp   Not operated (operated in wsb).
tk   Layer thickness (um)
kf   =0: No reference
     >0: Structure shape number referred in 
     kf designation field. The referred shape 
     structures are overwritten on the layer. 
     This numbers are represented by four 
     digits, up to 100 set per line.
