kstp =0:߂ꂽ`dnt܂ŌvZp
     =1:Wl肵0.001
        _ŌvZI
kskp xzo(i_*_t.out)̃XLbv
     傫قǍõ͓R}
     rȂ
lp   z(PML)̑w
     قǍǖʔ˂͑傫Ȃ
clp  xyǖʂ̋E
     (1:,2:)
     =0:PML  =1:PBC
crn  N[w
    iN[lɑ΂鎞Ԕ䗦j
     ʏ0.99x,Uޗł͔U}
     ߏ(0.9x)ɐݒ
kfl  Wscntł̕\͈
     =0:PML܂܂Ȃ\  =1:PML܂ޕ\
kot  xAzAܗzȂǂ̏o͂ő
     5+kotŕ\
ity  xz̒`A=0:Poynting Vectoȓ
     , =1:dꋭx, =2:dꋭx, =3:
     ꋭx
---------------------------------------------
kpls ʌ̔U
     = 0:CW
         <twcos, =twpeak, >twň
         ʕ͂ɍœK
     = 1:Gauss pulse
         =twpeak, tw/21/elS
     = 2:sin**2 pulse
         =twpeak, tw/21/2lS
     = 3:sin**3 pulse
         =0.5*tw,1.5*twŁ}peak,
         epeaktw/21/8lS
     = 4:Gauss-env.sinϒ
         =twpeak, tw/21/elS
          Ȃg͂ɓK
     = 5:Gauss-env.cosϒ
         =twpeak, tw/21/elS
tw   s[NʒuApX(m,`ŊZ)
kdip ʌ̕˕
     =0:(EHU)
     =1:o(EU)
     =2:o(HU)
kdr  ʌ̍L
     =0:z(xy)
     =1:x(yz)
     =2:y(xz)
dnt  `d(^um), Xebv~ԍ
nd1  CWUikpls=0j̏ꍇA̓`ɔ
     [ߖTŐUWl[
     񐔁B
nd2  PulseUikpls<>0j̏ꍇA̓`
     [ɑQ߂Wl̉l
     (10^nd2)B
--------------------------------------------
ksct UE̗L(kdr=0̏ꍇɋ@\)
     = 0:Ȃ
     = 1:
lx   xUË̕(w)
ly   yUË̕(w)
lz   zUË̕(w)
--------------------------------------------
(CWU(kpls=0)̏ꍇɋ@\)
kff  360xE̗L
     =0:͂Ȃ   =1:͂B
nff  ͂鉓E̕ʊpi])B
     ʂ360far.outɏo
thf  Eʊp̏lixj
fif  EΊpixj

(PulseU(kpls0)̏ꍇɋ@\)
krm  gXyNg̗͂L
     =0:͂Ȃ   =1:͂B
knm  gXyNg𑜐
     UEʂƎwޗ̂Uʂł̎Ԏ
     Ǝg͂flow_t.out ,flow_f.out
     o
rm1  ͔gl(m)
rm2  ͔gŏIl(m)

--------------------------------------------
wdx  ͕(m)B wx=02ɁB
     ̒Sƍ\̈ʒůɂȂB
     ObhԐnx=int(wx/dx)
wdy  ͕(m)Bwy=02ɁB
     ̒Sƍ\̈ʒůɂȂB
     ObhԐny=int(wy/dy)
dx   ObhԊu(m)
dy   ObhԊu(m)
dz   ObhԊu(m)
     tk(w)/dzɂȂȂwłtk/dz
     ؂グtk𓙕ԊuɂȂ
---------------------------------------------
Lam  ̔g(m)
th   ˕ʔg̕ʊp(deg)AƂȂp
fi   ˕ʔg̕Ίp(deg)
     kdr=0 : xyʓxƂȂp
     kdr=1 : yzʓyƂȂp
     kdr=2 : xzʓxƂȂp
gm   ˕ʔg̕ΌŁA
     dExNgEƂȂp(deg)
     kdr=1: xyʏEx̂Ȃp
     kdr=2: yzʏEŷȂp
     kdr=3: xzʏEx̂Ȃp
--------------------------------------------
wx0  ̂g(ʂ)C=0Ő
     kdr=2̏ꍇ͂g
wy0  ̂g(ʂ), =0Ő
     kdr=2̏ꍇ͂g
xrm  =1:lx
     =0:cosxz,xlS=wx0/2
     =0`1:Gaussianzxx       
yrm  =1:lx
     =0:cosxz,xlS=wy0/2
     =0`1:Gaussianzyx
sx0  ŜW(m)
sy0  ŜW(m)
kpx  ÛLȐ̎ށB
     kpx=0:O[h`A
     kpx=1:P[h`A
     kpx=2:Q[h`B
kpy  ÛLȐ̎ށB
     kpy=0:O[h`A
     kpy=1:P[h`A
     kpy=2:Q[h`B
--------------------------------------------
stx  S̍\Sʒu(m)
     ʒuɂ͕sKp
sty  S̍\Sʒu(m)
     ʒuɂ͕sKp
csx  Otfʂ̂ʒu(m)
csy  Otfʂ̂ʒu(m)

---------kmw藓(wޗw)-------------
     擪4͒ʂsԍ,200s܂œ͉\
Name ޗ(8)B
     SiO2,Ag,Al,Au,Be,Cr,Cu,Ni,Pd,Pt,Ti,W
     f[^݂w萔QƂ邪
     FDTDZŔ񕪎UA"|"ꍇ
     gUޗƂČvZBȊO
     Of[^Ƃnk.dat̃t@Cɔg
     AܗAW̗Lڂ邱
     ŋܗAWIɓ}B
     nk.datɃf[^݂Ȃꍇɂan
     ̒`lDBnk.dat͊e쐬̏A
     wsf.exeƓtH_Ɋi[̂ƁB
ko   wsf1.outւ̌oʏo͗L̎w
     =0:ʏo͂Ȃ   =1: ʏo͂
an   ܗ
ab   AbxA=0̏ꍇ,UȂ(ܗŒ)
ak   W

---------kdw藓ï搧`j------------
    擪4͒ʂsԍ,1000s܂œ͉\
kd  \`͕̓
    =0:`
    =1,-1:O\f[^ɂ`isub.datj
       wx,wyȊOKp,O\f[^400
       ܂Łie1000C܂Łj͉\
kt  \`̑I
    kd=1̏ꍇAkt=sub.datɂp^[No.
    kd=-1̏ꍇAkd=1̔]p^[
    kd=0̏ꍇi-ktkt̔]p^[)
      kt=0 ̈搧Ȃ
        =1 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒`ŗ̈搧
        =2 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̑ȉ~ŗ̈搧
           xp͑ȉ~wA-2.0<xp<-1.0Ő`,
           -1.0ŕH`,0.0ŉ~(ȉ~),>0.0ŕ
           `
        =3 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̂Up`(㉺p)
           ŗ̈搧
        =4 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̂Up`(Ep)
           ŗ̈搧
        =5 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̕H`ŗ̈搧
        =6 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒pRp`(Ε
           Pی)ŗ̈搧
        =7 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒pRp`(Ε
           Qی)ŗ̈搧
        =8 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒pRp`(Ε
           Rی)ŗ̈搧
        =9 sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒pRp`(Ε
           Sی)ŗ̈搧

ps   ̈撆S̐}`̉]p(deg)
px   \̂sb`(m)
     =0̎͌Ǘp^[
py   \̂sb`(m)
     =0̎͌Ǘp^[
wx   \̂̕(m) 
wy   \̂̕(m) 
sx   \̂WVtg(m)
sy   \̂WVtg(m)
xp   kt=2ł͑ȉ~w
     kt=10`17ł͊iqduty

---------kfw藓iOi\j------------
    擪4͒ʂsԍ,9999s܂œ͉\
km  \̍\ޗwikmw藓̍ޗNo.Q)
    km=0͐^(n=1.0)
kr  ̈搧wikrw藓̒ʂNo.Q)
    kr=0͎QƂ
kd   \`͕̓
   =0:`
   =1,-1:O\f[^ɂ`isub.datj
      wx,wyȊOKp,O\f[^400
     ie1000Cj܂œ͉\
kt  \`̑I
    kd=1,-1̏ꍇAkt=sub.datɂp^[No.
    kd=-1̏ꍇAkd=1̔]p^[
    kd=0̏ꍇi-ktkt̔]p^[)
      kt=0:̈Ȃ
        =1:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒`
        =2:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̑ȉ~Axp͑ȉ~
           wA-2.0<xp<-1.0Ő`,-1.0ŕH
           `,0.0ŉ~(ȉ~),>0.0ŕ`
        =3:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̂Up`(㉺p)
        =4:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̂Up`(Ep)
        =5:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̕H`
        =6:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒p3p`(Εӑ
           1ی)
        =7:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒p3p`(Εӑ
           2ی)
        =8:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒p3p`(Εӑ
           3ی)
        =9:sb`(px*py)̕`iqʒu𒆐S
           Ƃ镝(wx*wy)̒p3p`(Εӑ
           4ی)
        =10:sb`(px*py)̕`iq̊e`
           Ƀsb`(wx),px(wy),dutyxp,
           N_xq̒iq
        =11:sb`(px*py)̕`iq̊e`
           Ƀsb`(wx*wy),dutyxp,N_xq
           ̓Sȉ~iq
        =12:sb`(px*py)̕`iq̊e`
           Ƀsb`(wx*wy),dutyxp,N_xq
           ̓S12p`
        =13:1215x]}`
        =14:sb`(px*py)̕`iq̊e`
           Ƀsb`(wx*wy),dutyxp,N_xq
           ̓S18p`
        =15:1410x]}`
        =16:sb`(px*py)̕`iq̊e`
           Ƀsb`(wx*wy),dutyxp,N_xq
           ̓S6p`
        =17:1630x]}`
ps   ̈撆S̐}`̉]p(deg)
px   \̂sb`(m)
     =0̎͌Ǘp^[
py   \̂sb`(m)
     =0̎͌Ǘp^[
wx   \̂̕(m) 
wy   \̂̕(m) 
sx   \̂WVtg(m)
sy   \̂WVtg(m)
xp   kt=2ł͑ȉ~w
     kt=10`17ł͊iqduty
xq   kt=10`17Ŋiq̋N_

---------wiw-----------------------------
     10000C܂œ͉\AcC
     Ƃ܂œǂݍ܂Bŏw͍̏
     wAŉw͍̉ŉwƓw萔
     ̍\ŏ㉺ʂ̋E˂ȂB
kl   =1Fʒuw,w̒Ԑʒu
     S0̏ꍇcos(th)>0͑1w̏[A
     cos(th)<0͍ŏIw̉[ʒu
     =2Fw[xyfʋxzi_xy_t.out
     ɏo(ʒufʂ̋xzɉ)
km   w̍\ޗwikmw藓̍ޗNo.Q)
     km=0͐^(n=1.0)
kp   @\(wsbŎgp) 
tk   w̌(m)
kf   =0:QƂȂ
     >0:KfkfԖڂ̍sQƁBwɂ̍s
        łꂽ\㏑B4J\
        1s100܂Ŏw\
