#
#

<node_configuration>

sizey = 16# number of cpus in the y direction

<time>

last = 2000# last timestep

c= .45# velocity of light in comp. units (this defines the timestep)

timespan= 86400# time, in seconds, available to run the problem

<grid>

mx0 = 800 # number of actual grid points in the x direction
my0 = 256# number of actual grid points in the y direction
mz0 = 256 # ... (ignored for 2D simulations)

<algorithm>

conserv = 1 # charge-conservative current deposition -- the only available option
highorder= 0 # 0 -- 2nd order FDTD field integrateion; 1 -- 4th order;
# don't use 1 for non-relativistic flows

Corr= 1.025 # correction for the speed of light

ntimes= 4# number of passes for smoothing filter (current)
cleanfld= 0# number of passes for smoothing filter (fields). don't use.
cleanint= 10# interval for field cleaning; don't use.

cooling= 0# cool particles? ; not implemented
acool= 10.# cooling parameter for particles
splitparts = 0# split particles to improve statistics?

<restart>

irestart= 0# 1 to restart the simulation from saved restart/*.d files.
intrestart= 1000# how often to save restart files. They overwrite previous *d files.
laprestart= 0# if different from 0, restart from a named restart file, saved at timestep laprestart
namedrestartint = 1000000 # interval for saving named restart files, like restart.lap01234.d

<output>

interval = 20# plot interval
torqint= 2000000# interval for outputs at different resolution (currently broken)
pltstart= 0# starting iteration for first plot

istep= 1# downsampling factor for grid output
istep1= 4# downsampling factor for grid output every torqint steps
stride= 10# particle stride for particle output

writetestpart= 0# write test particles?
selectprt= 0# re-trace the same selected particles?

<boundaries>

periodicx= 1# periodic boundaries in the x direction? Choose 0 for radiative boundaries.
periodicy= 1# periodic boundaries in the y direction?
periodicz= 1# periodic boundaries in the z direction?

<domain>

enlarge= 0# if 1, enlarge box in the x direction if injector is close to right wall?
movwin = 0# if 1, use moving window
shiftinterval = 20# how often to apply moving window (in steps)
shiftstart = 1000# at what step to start shifting moving window
movwingam = 5. # gamma factor of moving window. If > 10000, it moves at c.
# if < 1, it is interpreted as v/c.
<fields>

btheta= 85# bfield angle bphi=0 -> bz, bph=90 in x-y plane, bth=0-> parallel
bphi= 90#

<particles>

sigma= 0.# magnetization number (omega_c/omega_p)^2, including gamma0
maxptl0 = 6e7# max number of particles in the simulation
ppc0 = 10# number of particles per cell

delgam = 1.e-4# delta gamma (temperature control)
me= 1.# electron mass
mi= 1.# ion mass (actually mass to charge ratio)

gamma0= 100.# flow drift gamma. If < 1, interpreted as v/c.

c_omp= 8# electron skin depth in cells

<problem>
caseinit = 1 #can be used to select subcases of the problem. Not used.

#density_profile=0. #x,y,z dependent function distributing initial particle weight
temperature_ratio = 1 # T_e/T_i

external_fields = 1 # if nonzero, add external nonevolving fields to mover
sigma_ext = 0.1 # strength of external magnetization,(omega_c/omega_p)^2,
# including gamma0


user_part_bcs=1 # call particle_bc_user routine from user file, specify particle bcs like walls

wall = 0 # use reflecting wall? Position is set in user file.
wallgam = 0. # gamma of the moving reflecting wall. If < 1, read as v/c. Set to 0 or 1 to not move the wall.

mx0, my0, and mz0

MAXPTL0