hm   ハーモニックス数比, 絶対値が大きいほど
     正確だが計算負荷大。通常4.0〜7.0に設定。
     >0:ハーモニックス数は奇数値になる。
     ハーモニック数は|hm|*(wdx+2*wb)/Lam又は
     |hm|*(wdy+2*wb)/Lamの切り上げ値。
     <0:ハーモニックス数は偶数値(上記奇数値-1)
     になる。
trc  切り捨て係数>=0, 3次元問題でのみ機能する。
     0では切り捨てなしとして処理。
     大きいほど正確だが計算負荷大。
     通常1.0〜3.0に設定。
wb   吸収境界幅(μm)。 =0は吸収壁なし。
     大きいほど理想的な無反射吸収になるが、
     計算量増大。
kfl  =0:吸収壁を含まない表示
     =1:吸収壁を表示
kot  強度、吸収、屈折率分布などの出力を最大
     5+kot桁で表現
ity  光強度分布の定義、=0:Poynting Vectorの大
     きさ, =1:電磁場強度, =2:電場強度, =3:磁
     場強度
--------------------------------------------
wdx  ｘ方向解析幅(μm)。wdx=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
     大きいほど大きなハーモニックス数が必要。
wdx  ｙ方向解析幅(μm)。wdy=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
     大きいほど大きなハーモニックス数が必要。
dxy  グリッド間隔の目安。実際の間隔はｘ、ｙ
     方向ともこれに近い値で最適化され、
     wsr.outのdx/dyに表示される。実際値が波
     長の1/10以下となることが好ましい。
dz   ｚ方向グリッド間隔(μm)
---------------------------------------------
Lam  波長(μm)
th   入射光線の方位角(deg)、ｚ軸となす角
fi   入射光線の方位偏角(deg)
     xy面内でx軸となす角
gm   光源の偏光方向(deg)
     電界ベクトルがｘ軸となす角
--------------------------------------------
alx  光源のｘ方向広がり
      =1.0：一様分布
      =0.0：x幅極小
      =0〜1：x方向の強度半値全幅=wdx*alx
aly  光源のy方向広がり
      =1.0：一様分布
      =0.0：y幅極小
      =0〜1：y方向の強度半値全幅=wdy*alx
sx0  光源中心のｘ座標(μm)
sy0  光源中心のｙ座標(μm)
--------------------------------------------
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は各自作成の上、wsr.exeと
     同じフォルダに格納のこと。
ko   wsr1.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で格子の起点

---------背景層-----------------------------
     先頭4桁は通し行番号,10000ラインまで入力
     可能、cラインが入るところまで読み込まれ
     る。最上層の上は最上層、最下層の下は最
     下層と同じ光学定数の構成で上下面の境界
     反射ない。
kl   =1：光源位置の層を指定
     全て0の場合はcos(th)>0は第1層の上端、
     cos(th)<0は最終層の下端が光源位置
km   層の構成材料指定（km指定欄の材料No.参照)
     km=0は真空(n=1.0)
kp   未機能(wsbで機能)
tk   層の厚さ(μm)
kf   =0:参照なし
     >0:Kf欄でkf番目の行を参照。層上にこの行
        でされた構造を上書き。4カラム表現で
        1行100個まで指定可能
