/
input
112 lines (73 loc) · 3.82 KB
/
input
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#
#
# Tristan-mp input file
#
#
<node_configuration>
sizey = 8 # number of cpus in the y direction
<time>
last = 100000 # 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 = 1000 # number of actual grid points in the x direction
my0 = 256 # number of actual grid points in the y direction
mz0 = 1 # ... (ignored for 2D simulations)
<algorithm>
conserv = 1 # charge-conservative current deposition -- the only available option
highorder = 1 # 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 = 10000 # 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 = 500 # plot interval
torqint = 2000000 # interval for outputs at different resolution (currently broken)
pltstart = 0 # starting iteration for first plot
istep = 2 # downsampling factor for grid output
istep1 = 2 # downsampling factor for grid output every torqint steps
stride = 100 # particle stride for particle output
writetestpart = 0 # write test particles?
selectprt = 0 # re-trace the same selected particles?
<boundaries>
periodicx = 0 # periodic boundaries in the x direction?
periodicy = 1 # periodic boundaries in the y direction?
periodicz = 1 # periodic boundaries in the z direction?
<domain>
enlarge = 1 # 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.1 # magnetization number (omega_c/omega_p)^2, including gamma0
maxptl0 = 3e7 # max number of particles in the simulation
ppc0 = 4 # 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 = 15. # flow drift gamma. If < 1, interpreted as v/c.
c_omp = 10 # 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 = 0 # if nonzero, add external nonevolving fields to mover
sigma_ext = 0.008 # 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 = 1 # 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.