######### COMPILE-TIME PARAMETERS ########## ##### ARCHITECTURE arch= # bg for BlueGene, none for Linux ##### OPTIMIZATIONS make dbug version # debug, production: enable optimizations opts= # lap for lapack optimization in 3d maxwell bloch ##### FDTD COMPILE DATA N=2 # dimensionality of the problem (2,3) problem_type=fdtd+periodic_x#+dispersive#+maxbloch # fdtd (fdtd run), dispersive (add dispersion), all (allows 2D TE+TM runs at the same time, e.g., for casimir), maxbloch (enable four level maxwell bloch system, it requires thermal noise to be enabled), periodic_x (periodic along x, y and z also supported) ##### I/O output=raw # (raw) vtk to be implemented ##### SOURCES spatial=tfsf+Ey+single # + # point + field (Ex, Ey and Ez) # tfsf + field (Ey or Ex) + opts [single (single plane at z0)]; the tfsf assumes a nondispersive background # field (Ex or Ey) alone selects just the polarization # thermal + field (Ex or Ey in 2D, in 3D no need as all of them are enabled) ##### MEASUREMENT measurements=energy+probes+Ey#+farfield+fscatt #+poynting # probes (Ei, Hi, Ni) + farfield (+ bscatt to omit the plane at z1, + fscatt to omit z0) + energy + poynting ########## RUN-TIME PARAMETERS ########## ##### MPI OPTIONS procs=1 9 1 ##### FDTD GEOMETRY number_of_cells=200 200 150 # x y z number of unit cells box_size=1 1 4e-6 # physical size of simulation box along x y z media_filename=geometry.in ##### TIME & I/O base_dir=./res # directory where to save results CFL=2. # Courant factor, Dt=Dx/C0/CFL, with CFL>sqrt(number of dimensions)/n_min;, for air, CFL=2 is ok in 2-3 dimensions total_steps=20000 # steps used for performing final averages (total steps = equilibration + measure) io_screen=10000 # I/O on screen chkpoint_every=1000000 # checkpoint step chk_load=0 # 1 true, 0 false chk_numb=0 # chkpoint number: 0 or 1 ##### SOURCES waveform=spl+0.5+0.#cw+0.43+3+1e6#spl+0.2+0.#fcb+0.4+1+10+3+100#cw+2.54e-07+3+100000 # cw + wavelength + m-duration + n-duration # gaussian + wavelength + waist + peak time position # gaussian2 + wavelength + waist in wl # spl + wavelength + starting time # cw + wavelength + m-duration + n-duration # fcb + wl0 + wl1 + N + m-duration + n-duration (fill wl0, wl1 with N mnm cw) # TFSF SOURCE SPECIFIERS tfsf_box=0. 0.999 0.25 0.75 1e-6 3e-6 # xsta xend ysta yend zsta zend separated by just one space tfsf_power=1e-3 # average power of the source [W] at central frequency # POINT SOURCE SPECIFIERS Ex_origin=2e-6 0.7e-6 2e-6 # Ex source points lists x y z separated by a single space Ey_origin=0.5 3.5 2e-6 # Ey source points lists x y z separated by a single space Ez_origin=0.7 0.2 1 # Ez source points lists x y z separated by a single space # THERMAL SOURCE SPECIFIERS thermal_temperature=3000 # environment temperature in Kelvin thermal_wavelength_min=100e-9 # minimum wavelength generated by the source thermal_period=5000 # period of the source emitted by each wall; after 2*this value the source emitted by the wall repeats itself; ##### MEASUREMENTS # FAR FIELD !!!! IMPORTANT, IF ENABLED ALWAYS USE A TOTALSTEPS MULTIPLE OF CHKPOINTSTEP farfield_box=0.1 1.9 0.5 1.5 0.9 1.1 # far field box xsta xend ysta yend zsta zend separated by just one space n_angles=1 # number of angles (theta,phi) to process ff_file=angles.in # binary filename containing angles [theta,phi, theta,phi, ...] # PROBES probe_list=0.5.+.5 # probe point list (x y z) separated by +; probe are saved as p1_.bin,p2_.bin,... # probes owned by one processor are grouped in p_<0>_<1>_...__.bin and saved in [0-t0,1-t0,...,N-t0; # 0-t1,1-t1,...,N-t1, # ..................] # ENERGY & POYNTING erg_io=50000 # frequency of io save ##### UPML rmax=-16 # reflection error, designated as R(0) in Equation 7.62 orderbc=4 # order of the polynomial grading, designated as m in Equation 7.60a,b upml=0 20 20 # number of UPML layer along x y and z