blob: f053d717860fdb33ae357717a63b8ef03742ad5a (
plain) (
blame)
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
|
amr.n_cell = 64 64 128
amr.max_grid_size = 128
amr.max_level = 0
max_step = 60
geometry.dims = 3
geometry.prob_lo = -125.e-6 -125.e-6 -149.e-6
geometry.prob_hi = 125.e-6 125.e-6 1.e-6
warpx.cfl = 0.99
# Domain Boundary condition
boundary.field_lo = none none none
boundary.field_hi = none none none
# Use parser to build EB
# Note that for amrex EB implicit function, >0 is covered, =0 is boundary and <0 is regular.
# This sets the region from -12.5e-6 < x < 12.5e-6, -12.5e-6 < y < 12.5e-6, -8.65e-6 < z < -6.15e-6 to be "covered"
warpx.eb_implicit_function = "-max(max(max(x-12.5e-6,-12.5e-6-x), max(y-12.5e-6,-12.5e-6-y)), max(z-(-6.15e-5),-8.65e-5-z))"
diagnostics.diags_names = diag1
diag1.intervals = 20
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz
# Order of particle shape factors
algo.particle_shape = 1
#################################
############ PLASMA #############
#################################
particles.species_names = electrons
electrons.species_type = electron
electrons.injection_style = "NUniformPerCell"
electrons.num_particles_per_cell_each_dim = 1 1 1
electrons.xmin = -1.e-5
electrons.xmax = 1.e-5
electrons.ymin = -1.e-5
electrons.ymax = 1.e-5
electrons.zmin = -149.e-6
electrons.zmax = -129.e-6
electrons.profile = constant
electrons.density = 1.e23 # number of electrons per m^3
electrons.momentum_distribution_type = "constant"
electrons.uz = 2000. # uth the std of the (unitless) momentum
electrons.save_particles_at_xhi = 1
electrons.save_particles_at_eb = 1
|
