blob: e5fe21c73cc7051f8a9244de83d95f18386ca752 (
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
|
stop_time = 1.5e-13
amr.n_cell = 32 64
amr.max_grid_size = 256
amr.blocking_factor = 16
amr.max_level = 0
geometry.dims = 2
geometry.prob_lo = -50.e-6 -10.e-6
geometry.prob_hi = 50.e-6 0.e-6
# Boundary condition
boundary.field_lo = periodic pml
boundary.field_hi = periodic pml
warpx.cfl = .999
warpx.do_moving_window = 1
warpx.moving_window_dir = z
warpx.moving_window_v = 1.0 # in units of the speed of light
warpx.serialize_initial_conditions = 1
warpx.use_filter = 0
# Order of particle shape factors
algo.particle_shape = 1
particles.species_names = beam
particles.rigid_injected_species = beam
beam.charge = -q_e
beam.mass = m_e
beam.injection_style = "gaussian_beam"
beam.x_rms = 1.e-6
beam.y_rms = 1.e-6
beam.z_rms = .5e-6
beam.x_m = 0.
beam.y_m = 0.
beam.z_m = -5.e-6
beam.npart = 2000
beam.q_tot = -1.e-20
beam.momentum_distribution_type = "gaussian"
beam.ux_m = 0.0
beam.uy_m = 0.0
beam.uz_m = 1000.
beam.ux_th = 100.
beam.uy_th = 100.
beam.uz_th = 0.
beam.zinject_plane = 20.e-6
beam.rigid_advance = true
beam.addRealAttributes = orig_z
beam.addIntegerAttributes = center
beam.attribute.orig_z(x,y,z,ux,uy,uz,t) = z
beam.attribute.center(x,y,z,ux,uy,uz,t) = 1*(sqrt(x*x + y*y)<5.e-7)
# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 10000
diag1.diag_type = Full
|
