# stop_time = 5e-11
max_step = 1000
amr.n_cell = 128 128 128
amr.check_input = 0
warpx.numprocs = 2 1 2
amr.max_level = 0
geometry.dims = 3

my_constants.kp_inv = 10.e-6
my_constants.kp = 1. / kp_inv
my_constants.wp = clight * kp
my_constants.ne = wp^2 * m_e * epsilon0 / q_e^2

collisions.collision_names = collision1
collision1.type = pairwisecoulomb
collision1.species = plasma plasma
collision1.CoulombLog = 20

#my_constants.ne = 1e22
#my_constants.kp_inv = 5.314093261582035e-05
#my_constants.R = 0.3*kp_inv
#my_constants.Z = 1.41*kp_inv
#my_constants.zinj = 100.e-6
#my_constants.Rw = 0.1*kp_inv
#my_constants.Zw = 0.2*kp_inv
#my_constants.Zz = -5*kp_inv
#my_constants.nd = 10.*ne
#my_constants.nw = 100.*ne

geometry.prob_lo     = -100.e-6   -100.e-6   -75.e-6
geometry.prob_hi     =  100.e-6    100.e-6    20.e-6

boundary.field_lo = periodic periodic pml
boundary.field_hi = periodic periodic pml
boundary.particle_lo = periodic periodic Absorbing
boundary.particle_hi = periodic periodic Absorbing
# boundary.particle_lo = Absorbing Absorbing Absorbing
# boundary.particle_hi = Absorbing Absorbing Absorbing

#################################
############ NUMERICS ###########
#################################
algo.maxwell_solver = ckc
algo.particle_shape = 3
warpx.use_filter = 1

warpx.cfl = .999
warpx.do_moving_window = 1
warpx.moving_window_dir = z
warpx.moving_window_v = 1. # in units of the speed of light

particles.species_names = plasma beam

beam.charge = -q_e
beam.mass = m_e
beam.injection_style = "NUniformPerCell"
beam.zmin = -10.e-6
beam.zmax =  10.e-6
beam.xmin = -30.e-6
beam.xmax =  30.e-6
beam.ymin = -30.e-6
beam.ymax =  30.e-6
beam.profile = parse_density_function
beam.density_function(x,y,z) = "if((x**2+y**2)<30.e-6,ne,0.)"
# beam.profile = constant
# beam.density = ne
beam.num_particles_per_cell_each_dim = 2 2 1
beam.momentum_distribution_type = "constant" # "gaussian"
beam.uz = 1000000000.

plasma.charge = -q_e
plasma.mass = m_e
plasma.injection_style = "NUniformPerCell"
plasma.zmin = 19.e-6
plasma.zmax = 100.
plasma.xmin = -7.*kp_inv
plasma.xmax =  7.*kp_inv
plasma.ymin = -7.*kp_inv
plasma.ymax =  7.*kp_inv
plasma.profile = constant
plasma.density = ne
plasma.num_particles_per_cell_each_dim = 4 4 1
plasma.momentum_distribution_type = "constant"
plasma.do_continuous_injection = 1

###########################
####### DIAGNOSTICS #######
###########################
diagnostics.diags_names = xz_slice beams

xz_slice.intervals = 100
xz_slice.format = openpmd
xz_slice.diag_type = Full
# data type to save openpmd files as (obviously h5 is the only good option): 
xz_slice.openpmd_backend = bp
xz_slice.fields_to_plot = Ex Ez By
xz_slice.write_species = 1
xz_slice.species = beam
xz_slice.diag_lo = -1. 0. -1.
xz_slice.diag_hi =  1. 0.  1.

beams.intervals = 100
beams.format = openpmd
beams.diag_type = Full
# data type to save openpmd files as (obviously h5 is the only good option): 
beams.openpmd_backend = bp
beams.species = beam
beams.coarsening_ratio = 8 8 8
beams.fields_to_plot = part_per_cell