# RZ setup
warpx.n_rz_azimuthal_modes = 2
geometry.coord_sys = 1
geometry.is_periodic = 0 0
warpx.do_pml = 0

# Mesh refinement
amr.max_level = 0

# Box config
max_step = 60 # 70500 to get to 1 GeV
warpx.cfl = 0.99
amr.n_cell = 416 256
amr.max_grid_size = 64
amr.blocking_factor = 32
geometry.prob_lo =  0.0 -3.11e-5
geometry.prob_hi =  1.55e-5  0.0

# Diagnostics
diagnostics.diags_names = diag
diag.diag_type = Full
diag.period = 20
diag.species = e-_primary e- e+ #electrons ions
diag.e-_primary.variables = w ux uy uz
diag.e-.variables = w ux uy uz
diag.e+.variables = w ux uy uz
diag.electrons.variables = w
diag.ions.variables = w
warpx.verbose = 1

# Moving window
warpx.do_moving_window=1
warpx.moving_window_dir = z
warpx.moving_window_v = 1.0

# Algorithms and numerics
particles.use_fdtd_nci_corr = 1
algo.current_deposition = esirkepov
algo.charge_deposition = standard
algo.field_gathering = energy-conserving
algo.particle_pusher = vay
algo.maxwell_solver = yee
warpx.use_filter = 1
warpx.filter_npass_each_dir = 1 1 1
interpolation.nox = 3
interpolation.noy = 3
interpolation.noz = 3
warpx.verbose = 1
warpx.do_dive_cleaning = 0

# Species
particles.nspecies = 5
particles.species_names = e-_primary e- e+ #electrons ions

# openPMD beams
particles.rigid_injected_species = e-_primary e- e+
e-_primary.charge = -q_e
e-_primary.mass=m_e
e-_primary.z_shift=-2.002e-03
e-_primary.injection_style = external_file
e-_primary.q_tot = 9.9438e-11
e-_primary.injection_file = e-_primary_opmd.h5
e-_primary.zinject_plane=-2.e-06
#e-.charge = -q_e
#e-.mass=m_e
#e-.z_shift=-2.002e-03
#e-.injection_style = external_file
#e-.q_tot = 1.0581e-11
#e-.injection_file = e-_opmd.h5
#e-.zinject_plane=-2.e-06
#e+.charge = q_e
#e+.mass=m_e
#e+.z_shift=-2.002e-03
#e+.injection_style = external_file
#e+.q_tot = 1.0426e-11
#e+.injection_file = e+_opmd.h5
#e+.zinject_plane=-2.e-06

# Plasma species
#electrons.charge=-q_e
#electrons.mass=m_e
#electrons.injection_style=NUniformPerCell
#electrons.num_particles_per_cell_each_dim=2 4 2
#electrons.profile=parse_density_function
#electrons.density_function(x,y,z)="((x*x+y*y)< 2.299e-10 )*((z< 2.428e-07 )*0.0 +(z> 2.428e-07 )*(z< 2.502e-04 )*(z* 1.847e+28 + -4.485e+21 ) +(z> 2.502e-04 )*(z< 4.590e-03 )* 4.618e+24 ) +(z> 4.590e-03 )*(z< 4.840e-03 )*( -1.847e+28 *(z- 4.590e-03 )+ 4.618e+24 ) +(z< 4.840e-03 )*0.0)"
#electrons.momentum_distribution_type = constant
#electrons.ux=0.0
#electrons.uy=0.0
#electrons.uz=0.0
#electrons.do_continuous_injection=1
#ions.charge=q_e
#ions.mass=m_p
#ions.injection_style=NUniformPerCell
#ions.num_particles_per_cell_each_dim=2 4 2
#ions.profile=parse_density_function
#ions.density_function(x,y,z)="((x*x+y*y)< 2.299e-10 )*((z< 2.428e-07 )*0.0 +(z> 2.428e-07 )*(z< 2.502e-04 )*(z* 1.847e+28 + -4.485e+21 ) +(z> 2.502e-04 )*(z< 4.590e-03 )* 4.618e+24 ) +(z> 4.590e-03 )*(z< 4.840e-03 )*( -1.847e+28 *(z- 4.590e-03 )+ 4.618e+24 ) +(z< 4.840e-03 )*0.0)"
#ions.momentum_distribution_type = constant
#ions.ux=0.0
#ions.uy=0.0
#ions.uz=0.0
#ions.do_continuous_injection=1

# Mirrors
warpx.num_mirrors=1
warpx.mirror_z=-4.370e-06
warpx.mirror_z_width=4.855e-06
warpx.mirror_z_npoints=40