ncy  伝搬の反復数
    =1:往路
    =2:往路＋復路
    >2:往路＋復路+２回目往路・・
wb   吸収境界幅(μm)。=0 は吸収壁なし。
     大きいほど理想的な無反射吸収になるが、
     計算量増大。
kfl  =0:吸収壁を含まない表示
     =1:吸収壁を表示
kot  強度、吸収、屈折率分布などの出力を最大
     5+kot桁で表現
ms   グラフ分布出力の圧縮指数
     x軸,y軸上で描画ポイント数を1/2^msに圧縮
ity  往復計算の場合の光強度分布の定義。
    =0:非干渉光,=1:電磁場強度,=2:電場強度,
    =3:磁場強度
--------------------------------------------
wdx  ｘ方向解析幅(μm)。wdx=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
wdx  ｙ方向解析幅(μm)。wdy=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
dxy  x,y方向グリッド間隔の目安。グリッド数を
     優先しているため若干のずれがある。実際
     の間隔はx,y方向ともこれに近い値で最適化
     され、wsb.outのdx/dyに表示される。実際
     値が波長の1/10以下となることが好ましい。
dz   ｚ方向グリッド間隔(μm)
     tk(層厚)/dzが整数にならない層ではtk/dzを
     切り上げた整数でtkを等分割する間隔になる
---------------------------------------------
Lam  波長(μm)
th   入射光線の方位角(deg)、ｚ軸となす角
fi   入射光線の方位偏角(deg)
     xy面内でx軸となす角
kps  光源の偏光条件
     =0:ランダム偏光
     =1:P偏光
     =2:S偏光
     =3:45度偏光
     =4:135度偏光
--------------------------------------------
wx0  光源のｘ方向広がり幅(μm)
wy0  光源のｙ方向広がり幅(μm)
xrm  =1:一様強度
     =0:cos強度分布,強度半値全幅=wx0/2
     =0〜1:Gaussian分布のx方向リム強度比       
yrm  =1:一様強度
     =0:cos強度分布,強度半値全幅=wy0/2
     =0〜1:Gaussian分布のy方向リム強度比       
sx0  光源中心のｘ座標(μm)
sy0  光源中心のｙ座標(μm)
dfc  光源のﾃﾞｨﾌｫｰｶｽ収差(μm)、
     dfc>0で下側に上記定義の焦点がずれる
kap  光源の開口形状
     =0:方形, =1:楕円形
--------------------------------------------
stx  全体構造ｘ方向中心位置(μm)
     光源位置には不適用
sty  全体構造ｙ方向中心位置(μm)
     光源位置には不適用
csx  グラフ断面のｘ方向位置(μm)
csy  グラフ断面のｙ方向位置(μm)

---------km指定欄(光学材料指定)-------------
     先頭4桁は通し行番号,200行まで入力可能
Name 材料名（8桁内）
     SiO2,Ag,Al,Au,Be,Cr,Cu,Ni,Pd,Pt,Ti,Wは
     内部データが存在。それら以外は外部デー
     タとしてnk.datのファイルに波長、屈折率、
     消衰係数の列を記載することで屈折率、消
     衰係数が自動的に内挿される。nk.datにデ
     ータが存在しない場合にはan以下の定義値
     を優先。nk.datは各自作成の上、wsb.exeと
     同じフォルダに格納のこと。
ko   wsb1.outへの検出光量出力有無の指定
     =0:結果出力なし
     =1:結果出力あり
an   屈折率
ab   アッベ数、=0の場合,分散なし(屈折率固定)
ak   消衰係数

---------kr指定欄（領域制限形状）------------
    先頭4桁は通し行番号,1000行まで入力可能
kd  構造形状の入力方式
    =0:内部定義
    =1,-1:外部構造データによる定義（sub.dat）
       wx,wy以外が適用,外部構造データは400種
       まで（各種1000ラインまで）入力可能
kt  構造形状の選択
    kd=1の場合、kt=sub.datにおけるパターンNo.
    kd=-1の場合、kd=1の反転パターン
    kd=0の場合（-ktはktの反転パターン)
      kt=0 領域制限なし
        =1 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の長方形で領域制限
        =2 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の楕円内で領域制限
           xpは楕円指数、-2.0<xp<-1.0で星形,
           -1.0で菱形,0.0で円(楕円),>0.0で方
           形
        =3 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の６角形(上下頂角)
           で領域制限
        =4 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の６角形(左右頂角)
           で領域制限
        =5 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の菱形で領域制限
        =6 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角３角形(斜辺
           第１象限)で領域制限
        =7 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角３角形(斜辺
           第２象限)で領域制限
        =8 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角３角形(斜辺
           第３象限)で領域制限
        =9 ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角３角形(斜辺
           第４象限)で領域制限
ps   領域中心回りの図形の回転角(deg)
px   構造のｘ方向ピッチ(μm)
     =0の時は孤立パターン
py   構造のｙ方向ピッチ(μm)
     =0の時は孤立パターン
wx   構造のｘ方向の幅(μm) 
wy   構造のｙ方向の幅(μm) 
sx   構造のｘ座標シフト(μm)
sy   構造のｙ座標シフト(μm)
xp   kt=2では楕円指数
     kt=10〜17では格子のduty比

---------kf指定欄（前景構造）------------
    先頭4桁は通し行番号,9999行まで入力可能
km  構造の構成材料指定（km指定欄の材料No.参照)
    km=0は真空(n=1.0)
kr  領域制限指定（kr指定欄の通しNo.参照)
    kr=0は参照せず
kd   構造形状の入力方式
   =0:内部定義
   =1,-1:外部構造データによる定義（sub.dat）
      wx,wy以外が適用,外部構造データは400種
     （各種1000ライン）まで入力可能
kt  構造形状の選択
    kd=1,-1の場合、kt=sub.datにおけるパターンNo.
    kd=-1の場合、kd=1の反転パターン
    kd=0の場合（-ktはktの反転パターン)
      kt=0:領域なし
        =1:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の長方形
        =2:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の楕円内、xpは楕円
           指数、-2.0<xp<-1.0で星形,-1.0で菱
           形,0.0で円(楕円),>0.0で方形
        =3:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の６角形(上下頂角)
        =4:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の６角形(左右頂角)
        =5:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の菱形
        =6:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角3角形(斜辺第
           1象限)
        =7:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角3角形(斜辺第
           2象限)
        =8:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角3角形(斜辺第
           3象限)
        =9:ピッチ(px*py)の方形格子位置を中心
           とする幅(wx*wy)の直角3角形(斜辺第
           4象限)
        =10:ピッチ(px*py)の方形格子の各方形
           内にピッチ(wx),角度(wy),duty比xp,
           起点xqの直線格子
        =11:ピッチ(px*py)の方形格子の各方形
           内にピッチ(wx*wy),duty比xp,起点xq
           の同心楕円格子
        =12:ピッチ(px*py)の方形格子の各方形
           内にピッチ(wx*wy),duty比xp,起点xq
           の同心12角形
        =13:12の15度回転図形
        =14:ピッチ(px*py)の方形格子の各方形
           内にピッチ(wx*wy),duty比xp,起点xq
           の同心18角形
        =15:14の10度回転図形
        =16:ピッチ(px*py)の方形格子の各方形
           内にピッチ(wx*wy),duty比xp,起点xq
           の同心6角形
        =17:16の30度回転図形
ps   領域中心回りの図形の回転角(deg)
px   構造のｘ方向ピッチ(μm)
     =0の時は孤立パターン
py   構造のｙ方向ピッチ(μm)
     =0の時は孤立パターン
wx   構造のｘ方向の幅(μm) 
wy   構造のｙ方向の幅(μm) 
sx   構造のｘ座標シフト(μm)
sy   構造のｙ座標シフト(μm)
xp   kt=2では楕円指数
     kt=10〜17では格子のduty比
xq   kt=10〜17で格子の起点

---------背景層-----------------------------
     10000ラインまで入力可能、cラインが入る
     ところまで読み込まれる。最上層の上は最
     上層、最下層の下は最下層と同じ光学定数
     の構成で上下面の境界反射ない。
kl   未機能(wsf,wsrで使用)
km   層の構成材料指定（km指定欄の材料No.参照)
     km=0は真空(n=1.0)
kp   =0:dzの定義による自動設定
     >0:tkの厚さをkp分割で刻む、均一媒質の場
        合はkp=1の設定で飛ばし計算できる
tk   層の厚さ(μm)
kf   =0:参照なし
     >0:Kf欄でkf番目の行を参照。層上にこの行
        でされた構造を上書き。4カラム表現で
        1行100個まで指定可能
