ncy  Number of propagation iterations.
   =1:forward.
   =2:forward + backward.
   >2:forward + backward + 2nd forward..
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.
ms   Drawing point number in x/y directions 
     is compressed to 1/2^ms.
ity  Definition of intensity distribution
     for a round trip calculation.
    =0 : for incoherent light
         The interference between forward and 
         backward light is neglected and 
         overlapped light is treated as 
         intensity sum of them.
    =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.
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.
dxy  Rough grid interval in x/y-direction 
     (um). The actual interval is optimized 
     to be close to that and displayed in 
     wsb.out. The actual value should be 
     less than 1/10 of the wavelength.
dz   Grid interval in z-direction (um).
     For layers where the layer thickness tk 
     divided by dz is not an integer, the 
     grid interval becomes tk divided by an 
     integer obtained by rounding up tk/dz.
---------------------------------------------
Lam  Wavelength (um).
th   Azimuth angle of incident light with 
     z-axis (deg).
fi   Argument angle of incident light in 
     xy-plane (deg).
kps  Source polarization.
     = 0 : random polarization.
     = 1 : p-polarization.
     = 2 : s-polarization.
     = 3 : 45-degree polarization.
     = 4 : 135-degree polarization.
--------------------------------------------
wx0  Light source spread in x-direction (um).
wy0  Light source spread in y-direction (um).
xrm  =1 : uniform intensity in x-direction.
     =0 : cos-intensity distribution,
          full width half maximum = wx0/2
     =0～1 : rim intensity ratio of Gaussian 
          distribution in x-direction.          
yrm  =1 : uniform intensity in y-direction.
     =0 : cos-intensity distribution,
          full width half maximum = wy0/2.
     =0～1 : rim intensity ratio of Gaussian 
          distribution in y-direction.      
sx0  Shift length of light source center in 
     x-direction (um).
sy0  Shift length of light source center in 
     y-direction (um).
dfc  Defocus aberration of light source (um).
     If dfc>0, the focal point shifts to the 
     lower side.
kap  Apperture shape of light source.
     =0 : rectangular, =1 : elliptic shape.
--------------------------------------------
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 wsb.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   Not operated (operated in wsf and wsr).
km   Construction material number referred in 
     km designation field.
     km=0 means vacuum (n=1.0).
kp  =0:Automatic setting by definition of dz. 
    >0:The layer thickness tk is divided by
     kp. In case of homogeneous 
     configuration, setting of kp=1 enables 
     skip calculation.
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.

