########################################################################
#This is an example of the input file for F2P calculations. 
#A 90 degree waveguide bend in a square lattice of dielectric 
#rods in air.
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#------------------------------------------------Prepared by Dr. Min Qiu
########################################################################

########################################################################
#Calculate the TM or TE modes??
#              ==0: TE modes
#              ==1: TM modes
  1
########################################################################

########################################################################
#Define the materials used in the calculations.
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#The number of materials 
   2
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#The dielectric constant, and conductivity of each material
#Index---dielectric constant---conductivity
   1           1.0000            0.00E+0  #air
   2          11.5600            0.00E+0  #Silicon (n=3.4)
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#End of material definition.
########################################################################

########################################################################
#define the size of the computational domain
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#--- lattice constant (For normalized purpose)
          1.0    
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#The size of the computational domain.
#---Sx-------Sy--
    60.0    60.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Background materials index 
     1                                
########################################################################

########################################################################
#Absorbing boundary condition (ABC)
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#ABC type: =0 All surrounded by PML; =1 Only the Y direction surrounded 
#by PML, while periodic BC in the X direction.
#---ABC type----------------------------------------------------------
      0              
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#---number of PML layers in X----in Y---------The exponent index for PML
             101                   101                 2.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
########################################################################

########################################################################
#Define parameters for FDTD
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#----Dx---------Dy------
    0.05       0.05
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Dt coefficient. Dt=Dt_coe*1/sqrt(1/(Dx*Dx)+1/(Dy*Dy))/vc
#Should be less than 1.0. 
#For example, 0.95 or 0.866 are all good choices.
#---Dt_Coe---
    0.95     
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Number of total time steps 
    10000
########################################################################

########################################################################
#Output the z-field component in some time steps. 
#For snapshot purpose.
#If interval is less than 1, no output. 
# The current z-field distribution is always stored in "Fieldz_Bin.dat".
#----Start Time steps---------End Time steps------Interval--------------
          20                     10000              -20
########################################################################

########################################################################
#define the inclusions, could be waveguides, cylinders, ...
#-----------------------------------------------------------------------
#Inclusion type:
#  == 1: ellipse
#  == 2: quadrilateral
#  == 3: triangular
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#The total number of inclusions
    5
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Inclusion No. 1, square lattice of dielectric rods
#Inclusion Type ----Material index----
    1                    2  
#-----Rx---------Ry---------Xc---------Yc--------Theta--
     0.18      0.18        0.50       0.50       0.0   
#-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat)
      60       1.00      60      1.00     0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Inclusion No. 2, circles, waveguide down
#Inclusion Type ----Material index----
    1                    1  
#-----Rx---------Ry---------Xc---------Yc--------Theta--
     0.50      0.50        0.50       10.5       0.0   
#-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat)
      50       1.00      1      1.00     0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Inclusion No. 3, circles, waveguide after bend
#Inclusion Type ----Material index----
    1                    1  
#-----Rx---------Ry---------Xc---------Yc--------Theta--
     0.50      0.50        49.50      10.5       0.0   
#-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat)
      1       1.00       50      1.00     0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Inclusion No. 4, circles, corner 1
#Inclusion Type ----Material index----
    1                    2  
#-----Rx---------Ry---------Xc---------Yc--------Theta--
     0.18      0.18        49.50      10.5       0.0   
#-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat)
      1       1.00        1      1.00     0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Inclusion No. 5, circles, corner 2
#Inclusion Type ----Material index----
    1                    1  
#-----Rx---------Ry---------Xc---------Yc--------Theta--
     0.50      0.50        48.50      11.5       0.0   
#-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat)
      1       1.00       1      1.00     0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Examples of Triangular 
#Inclusion Type ----Material index----
#    3                    1  
#---X1------Y1------X2------Y2------X3-------Y3---
#   10.5    0.0    24.5    3.964   16.0    9.964   
#---Mx-----Lx-----My-----Ly----Alpha---(For Repeat)
#    1     0.00     1    0.00   0.0
########################################################################

########################################################################
#Source part
#-----------------------------------------------------------------------
#  For source type: ==1 point source
#                   ==2 line source, direction, from P1 to P2, in P3
#                   ==3 waveguide mode source
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#----Number of Source------------
     1
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Point Source Examples
#---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)-
#     1                    0.305                        0.007
#----XS-----YS-------- (Point source position)
#    4.0    6.06218
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Gaussian Line Source sample
#For line source, direction == +1(X); -1(-X); +2(+Y); -2(-Y);  
#---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)-
#     2                    0.25                       0.1
#----direction---XS-----YS-----Beam Width---Material Type----Theta------- 
#     1          9.00   6.06218     2.0            2           0.0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Waveguide Source sample
#For Waveguide source, direction == +1(X); -1(-X); +2(+Y); -2(-Y);  
#---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)-
#     3                   0.295                          0.007
#----direction---XS-----YS-----Waveguide Width---Mode Number-----------
#     1          1.00  6.06218   2.60            1
#---Material Up (Left)----Material Core----Material Down (Right)------- 
#     1                        2                 1 
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#A Point Source
#---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)-
     1                    0.38                        -0.10
#----XS-----YS-------- (Point source position)
    6.0    10.5
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
########################################################################

########################################################################
#detectors definition
#-----------------------------------------------------------------------
#  For detector type: ==1 point detector
#                     ==2 line detector
#                            direction == +1(X); -1(-X); +2(+Y); -2(-Y);
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#FT frequency range and points in between
#-W_min---------W_max-------NF---
   0.10        0.50       201
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#----Number of Detector------------
     0
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Line Detector Example
#---Detector Type-----Starting Time Steps-----Ending Time Steps--------
#     2                      0                    30000
#----direction-----Xc-------Yc---------Length------- 
#     1            2.0     6.06218       4.00
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
#Point Detector Example
#---Detector Type-----Starting Time Steps-----Ending Time Steps--------
#     1                      0                    500000
#-----Xc--------Yc-------
#   11.00     6.06218     
########################################################################