blob: fe6d307d618831f187cfb26af2cdb7f074969567 (
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
49
50
51
52
53
54
55
|
# Iterations
max_step = 400
# Domain
amr.n_cell = 32 32 32
amr.max_grid_size = 16
amr.max_level = 0
# Geometry
geometry.dims = 3
geometry.prob_lo = -50.e-6 -50.e-6 -50.e-6
geometry.prob_hi = 50.e-6 50.e-6 50.e-6
# Boundary condition
boundary.field_lo = pml pml pml
boundary.field_hi = pml pml pml
# Order of particle shape factors
algo.particle_shape = 3
# Numerics
warpx.cfl = 0.25
warpx.do_divb_cleaning = 1
warpx.use_filter = 1
# External magnetic field
my_constants.qm = 1e-1
warpx.B_ext_grid_init_style = parse_B_ext_grid_function
warpx.Bx_external_grid_function(x,y,z) = qm * x / (x*x + y*y + z*z)
warpx.By_external_grid_function(x,y,z) = qm * y / (x*x + y*y + z*z)
warpx.Bz_external_grid_function(x,y,z) = qm * z / (x*x + y*y + z*z)
# Particle beam
particles.species_names = beam
beam.charge = -q_e
beam.mass = 1.e30
beam.injection_style = "gaussian_beam"
beam.x_rms = 2.e-6
beam.y_rms = 2.e-6
beam.z_rms = 2.e-6
beam.x_m = 0.
beam.y_m = 0.
beam.z_m = 0.e-6
beam.npart = 20000
beam.q_tot = -1.e-20
beam.momentum_distribution_type = "constant"
beam.ux = 0.0
beam.uy = 0.0
beam.uz = 0.5
# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 0, 398:400:1
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz divE divB G
|
