readme

一维光子晶体的能带结构模拟

 Here is the input data file format
 Build V1.0.1
 Date: 6/26/02
 By: Shangping Guo
 [e-mail]
 1-757-4231507

 PML2DTM v1.01
 PML2DTE v1.01
 are compiled and run on Solaris, copy and run.
 tar -xvf *.tar to extract all files in tar
 Note: 
 1. follow format strictly, else the data will not be read
 correctly and program will not run correctly
 2. Verification is written in echo2dtd.prb
 3. Lower/upper case matters
 4. no blank lines between data
 5. no comments allowed between data

 !**sampletstep, integer, the time interval for writing into GIF file
 !**negative means no writing**

 -1

 !**nb,nu,R0, 
 !**nb: integer, number of PML layers**
 !**nu: double, the order of PML function
 !**R0: double, the normal reflection coefficient, affect a lot

 10 2.0 1.0E-10

 !**nx,dx,ny,dy***
 !**nx,ny: integer, number of cells along x and y direction
 !**dx,dy: double, size of cell along x and y direction

 101 0.01 101 0.01

 !**dt,nt**
 !**dt: double, the time step, should follow stability check
 !**nt: integer, number of time steps to be calculated

 1.6667E-11 400

 !**x1,y1,x2,y2,srcshape,wave,amp,freq,t0,tw,srctype,tsrcstop**
 !**to describe the source
 !**x1,y1,x2,y2: integer, the 2 points of a rectangle.
 !**x1,y1: the coordinate of the lower-left point
 !**x2,y2: the coordinate of the upper-right point
 !**srcshape: how the source is distributed in the rectangle,
 !**          UN: uniform
 !**wave: only support these 3
	sine: a sine wave
	gauss_lp: gauss pulse
	gauss_bp: gauss pulse modulated with a sine wave
!**amp: double, amplitude
!**freq: double, frequency
!**t0: time delay, any for sine,should put 0
!**tw: pulse width, any for sine,should put 0
!**srctype: only have 2 value:
	hard: hard source
	soft: treat as a soft transparent source
!**tsrcstop: double, when to stop the source, (in unit of tw)
!** 		0: will not stop during simulation

 gauss_bp        1000.0 477460000.0 1.3333E-9 5.0E-10 soft           

 !**numofmats: integer, number of materials**

 3

 !**mater,epsr,mur**
 !**mater: material number, 100 is for source, 200 is for PEC
 !**epsr: double, epsr
 !**mur: double, mur

 1 1.0 1.0
 100 1.0 1.0
 200 1.0 1.0

 !**numofregions: number of regions**
 1
 !**mater,shape,x1,y1,x2,y1,patchtype,radius**
 !**mater: material number, refer above
 !**shape: the shape of this region
	   R: rectangle
	   D: disk, including circle and ellipse
 !**x1,y1,x2,y2, define the rectangle region
 !**patchtype: the sub-type in the rectangle region
	   NA: rectangle
	   LL: low-left triangle
	   UL: upper-left triangle
	   LR: lower right triangle
	   UR: upper-right triangle
!**radius: only used for disk, I am now considering to remove it
 1   R 1 1 101 101  NA 0.0E+0
 !**number of source rect
 1
 !**srcpatch, x1,y1,x2,y2,patchtype**
 R  51 51 51 51 UN
 !**PCsupport: only 2 valuse support: PC:YES and PC:NO**
 PC:NO   
 !**if photonic crystal is present, ie. PC:YES, then the following data
 !**is needed for input
 !**1. numofatoms: integer
 !**2. mater, llx,lly,urx,ury,type,x_ext,y_ext
 !**mater: material number of the atom
 !**llx,lly,urx,ury: defines the rectangle region for the first atom 
 !**			on the lower-left 
 !**type: the shape of the atom,"R" and "D" is now supported
 !**x_ext,y_ext: how many same atoms are going to place in u1 and 
 !** u2 direction (U1 and U2 are the unit axis of the photonic crystal)
 !**numofplotfiles**
 11
 !**plottype,field,x1,y1,x2,y2,tstep**
 !**plottype: T: time series
	      R: field in a rectangle
	      M: material distribution in a region
	      P: power flux
 !**fieldtype: Ex, Ey, Ez, Hx, Hy,Hz
 !**x1,y1,x2,y2,define the region to output
 !**tstep: the frequency to record
 T Ez 1 51 1 1 1
 T Ez 11 51 1 1 1
 T Ez 21 51 1 1 1
 T Ez 31 51 1 1 1
 T Ez 41 51 1 1 1
 T Ez 51 51 1 1 1
 T Ez 61 51 1 1 1
 T Ez 71 51 1 1 1
 T Ez 81 51 1 1 1
 T Ez 91 51 1 1 1
 T Ez 101 51 1 1 1