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