kstp =0:決められた伝搬距離dntまで計算を継続
     =1:安定係数値が安定して0.001を下回った
        時点で計算を終了
kskp 光強度分布出力(i_*_t.out)のスキップ数
     大きいほど高速だが出力動画のコマ送りは
     荒くなる
lp   吸収壁(PML)の層数
     小さいほど高速だが壁面反射は大きくなる
clp  xy両壁面の境界条件
     (1桁目:ｘ方向,2桁目:ｙ方向)
     =0:PML  =1:PBC
crn  クーラン指数
    （クーラン条件値に対する時間比率）
     通常は0.99程度,分散材料では発散を抑える
     ため小さめ(0.9程度)に設定
kfl  Wscntでの表示範囲
     =0:PMLを含まない表示  =1:PMLを含む表示
kot  強度、吸収、屈折率分布などの出力を最大
     5+kot桁で表現
ity  光強度分布の定義、=0:Poynting Vectorの大
     きさ, =1:電磁場強度, =2:電場強度, =3:磁
     場強度
---------------------------------------------
kpls 面光源の発振条件
     = 0:CW光源
         <twでcos増加, =twでpeak, >twで一定
         光量分析に最適
     = 1:Gauss pulse光源
         =twでpeak, tw/2の1/e値全幅
     = 2:sin**2 pulse光源
         =twでpeak, tw/2の1/2値全幅
     = 3:sin**3 pulse光源
         =0.5*tw,1.5*twで±peak,
         各peakでtw/2の1/√8値全幅
     = 4:Gauss-env.sin変調光源
         =twでpeak, tw/2の1/e値全幅
          直流がなく周波数特性分析に適切
     = 5:Gauss-env.cos変調光源
         =twでpeak, tw/2の1/e値全幅
tw   ピーク位置、パルス幅(μm,伝搬距離で換算)
kdip 面光源の放射方向
     =0:一方向放射(EH励振)
     =1:双方向放射(E励振)
     =2:双方向放射(H励振)
kdr  面光源の広がり方位
     =0:z軸直交面(xy面)
     =1:x軸直交面(yz面)
     =2:y軸直交面(xz面)
dnt  伝播長(真空中um), ステップ数×時間刻み
nd1  CW発振（kpls=0）の場合、光の伝搬に伴い
     ゼロ近傍で振動する安定係数値がゼロを跨
     ぐ回数。
nd2  Pulse発振（kpls<>0）の場合、光の伝搬に
     伴いゼロに漸近する安定係数値の下限値
     (10^nd2)。
--------------------------------------------
ksct 散乱界の有無(kdr=0の場合に機能)
     = 0:なし
     = 1:あり
lx   x軸方向散乱界領域の幅(層数)
ly   y軸方向散乱界領域の幅(層数)
lz   z軸方向散乱界領域の幅(層数)
--------------------------------------------
(CW発振(kpls=0)の場合に機能)
kff  360度遠方界の有無
     =0:分析なし   =1:分析あり。
nff  分析する遠方界の方位角数（一回転分割数)。
     結果は360far.outに出力
thf  遠方界方位角の初期値（度）
fif  遠方界偏角（度）

(Pulse発振(kpls≠0)の場合に機能)
krm  周波数スペクトル分析の有無
     =0:分析なし   =1:分析あり。
knm  周波数スペクトル分割解像数
     ６境界面と指定材料の６側面での時間軸分析
     と周波数軸分析をflow_t.out ,flow_f.outに
     出力
rm1  分析波長初期値(μm)
rm2  分析波長最終値(μm)

--------------------------------------------
wdx  ｘ方向解析幅(μm)。 wx=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
     グリッド間数はnx=int(wx/dx)
wdy  ｙ方向解析幅(μm)。wy=0で2次元問題に。
     幅の中心が光源と構造の位置の基準になる。
     グリッド間数はny=int(wy/dy)
dx   ｘ方向グリッド間隔(μm)
dy   ｙ方向グリッド間隔(μm)
dz   ｚ方向グリッド間隔(μm)
     tk(層厚)/dzが整数にならない層ではtk/dzを
     切り上げた整数でtkを等分割する間隔になる
---------------------------------------------
Lam  光源の波長(μm)
th   入射平面波の方位角(deg)、ｚ軸となす角
fi   入射平面波の偏角(deg)
     kdr=0 : xy面内でx軸となす角
     kdr=1 : yz面内でy軸となす角
     kdr=2 : xz面内でx軸となす角
gm   入射平面波の偏光方向で、
     電界ベクトルEが軸となす角(deg)
     kdr=1: xy面上でEとx軸のなす角
     kdr=2: yz面上でEとy軸のなす角
     kdr=3: xz面上でEとx軸のなす角
--------------------------------------------
wx0  光源のｘ方向拡がり(μｍ)，=0で線光源
     kdr=2の場合はｚ方向拡がり
wy0  光源のｙ方向拡がり(μｍ), =0で線光源
     kdr=2の場合はｚ方向拡がり
xrm  =1:一様強度
     =0:cos強度分布,強度半値全幅=wx0/2
     =0〜1:Gaussian分布のx方向リム強度比       
yrm  =1:一様強度
     =0:cos強度分布,強度半値全幅=wy0/2
     =0〜1:Gaussian分布のy方向リム強度比
sx0  光源中心のｘ座標(μm)
sy0  光源中心のｙ座標(μm)
kpx  光源振幅のｘ方向広がり曲線の種類。
     kpx=0:０次モード形状、
     kpx=1:１次モード形状、
     kpx=2:２次モード形状。
kpy  光源振幅のｙ方向広がり曲線の種類。
     kpy=0:０次モード形状、
     kpy=1:１次モード形状、
     kpy=2:２次モード形状。
--------------------------------------------
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は
     内部データが存在し光学定数を参照するが
     FDTD演算で非分散処理、"−"がつく場合は
     周波数分散材料として計算。それら以外は
     外部データとしてnk.datのファイル内に波
     長、屈折率、消衰係数の列を記載すること
     で屈折率、消衰係数が自動的に内挿される。
     nk.datにデータが存在しない場合にはan以
     下の定義値を優先。nk.datは各自作成の上、
     wsf.exeと同じフォルダに格納のこと。
ko   wsf1.outへの検出光量出力有無の指定
     =0:結果出力なし   =1: 結果出力あり
an   屈折率
ab   アッベ数、=0の場合,分散なし(屈折率固定)
ak   消衰係数

---------kd指定欄（領域制限形状）------------
    先頭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   =1：光源位置を指定,層の中間線が光源位置
     全て0の場合はcos(th)>0は第1層の上端、
     cos(th)<0は最終層の下端が光源位置
     =2：層上端のxy断面強度分布をi_xy_t.out
     に出力(光源位置断面の強度分布に加えて)
km   層の構成材料指定（km指定欄の材料No.参照)
     km=0は真空(n=1.0)
kp   未機能(wsbで使用) 
tk   層の厚さ(μm)
kf   =0:参照なし
     >0:Kf欄でkf番目の行を参照。層上にこの行
        でされた構造を上書き。4カラム表現で
        1行100個まで指定可能
