amr.n_cell = 128 128 128

hipace.normalized_units = 1

my_constants.kp_inv = 1. # 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

amr.blocking_factor = 2
amr.max_level = 0

max_step = 0
hipace.output_period = 1
# hipace.do_tiling = 0

hipace.numprocs_x = 1
hipace.numprocs_y = 1

hipace.depos_order_xy = 2

geometry.coord_sys   = 0                  # 0: Cartesian
geometry.is_periodic = true  true  false  # Is periodic?
geometry.prob_lo     = -10*kp_inv   -10*kp_inv   -7.5*kp_inv    # physical domain
geometry.prob_hi     =  10*kp_inv    10*kp_inv    2*kp_inv

beams.names = beam
beam.injection_type = fixed_ppc
beam.profile = flattop
beam.zmin = -kp_inv
beam.zmax = kp_inv
beam.radius = 3*kp_inv
beam.density = 0.01 * ne
beam.u_mean = 0. 0. 2000
beam.u_std = 0. 0. 0.
beam.ppc = 1 1 1

plasmas.names = plasma # plasma2
plasma.density(x,y,z) = ne
plasma.ppc = 8 8
plasma.u_mean = 0.0 0.0 0.
plasma.element = electron
plasma2.density(x,y,z) = ne
plasma2.ppc = 8 8
plasma2.u_mean = 0.0 0.0 0.
plasma2.element = proton

# plasmas.collisions = collision1
# collision1.species = plasma plasma

diagnostic.diag_type = xz