# Maximum number of time steps
max_step = 2

# number of grid points
amr.n_cell = 16 16 16

# The lo and hi ends of grids are multipliers of blocking factor
amr.blocking_factor = 8

# Maximum allowable size of each subdomain in the problem domain;
#    this is used to decompose the domain for parallel calculations.
amr.max_grid_size = 16

# Maximum level in hierarchy (for now must be 0, i.e., one level in total)
amr.max_level = 0

# Geometry
geometry.dims = 3
geometry.prob_lo     = -4 -4 -4
geometry.prob_hi     =  4  4  4

# Deactivate Maxwell solver
algo.maxwell_solver = none
warpx.const_dt = 7e-9

# Boundary condition
boundary.field_lo = periodic periodic periodic
boundary.field_hi = periodic periodic periodic

# particles
particles.species_names = electron proton
algo.particle_shape = 3

electron.charge = -q_e
electron.mass = m_e
electron.injection_style = NFluxPerCell
electron.num_particles_per_cell = 100
electron.surface_flux_pos = -4.
electron.flux_normal_axis = y
electron.flux_direction = +1
electron.flux_profile = parse_flux_function
electron.flux_function(x,y,z,t) = "1."
electron.momentum_distribution_type = gaussianflux
electron.ux_th = 0.1
electron.uy_th = 0.1
electron.uy_m = 0.07
electron.uz_th = 0.1

proton.charge = +q_e
proton.mass = m_p
proton.injection_style = NFluxPerCell
proton.num_particles_per_cell = 100
proton.surface_flux_pos = 4.
proton.flux_normal_axis = x
proton.flux_direction = -1
proton.flux_profile = constant
proton.flux = 1.
proton.momentum_distribution_type = gaussianflux
proton.ux_th = 0.1
proton.ux_m = 0.05
proton.uy_th = 0.1
proton.uz_th = 0.1

# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 1000
diag1.diag_type = Full
diag1.fields_to_plot = none