max_step = 200
amr.n_cell = 128 128 128

amr.blocking_factor = 16
amr.max_grid_size = 1024
amr.max_level = 1

# Geometry
geometry.dims = 3
geometry.prob_lo     = -32.e-6  -32.e-6  -32.e-6      # physical domain
geometry.prob_hi     =  32.e-6   32.e-6   32.e-6
warpx.fine_tag_lo     = -8.e-6  -8.e-6  -8.e-6      # physical domain
warpx.fine_tag_hi     =  8.e-6   8.e-6   8.e-6

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

# PML
warpx.pml_ncell = 12
warpx.pml_delta = 6
warpx.pml_has_particles = 1
warpx.do_pml_in_domain = 1
warpx.do_pml_j_damping = 1


# Algorithms
algo.current_deposition = esirkepov
algo.charge_deposition = standard
algo.particle_pusher = vay
algo.maxwell_solver = ckc
warpx.cfl = 0.9
warpx.use_filter = 1

# Particle species
particles.species_names = electron proton

electron.charge = -q_e
electron.mass = m_e
electron.injection_style = "singleparticle"
electron.single_particle_pos = 0. 0. 0.
electron.single_particle_u = 2. 0. 0.
electron.single_particle_weight = 1.

proton.charge = q_e
proton.mass = m_p    # Very heavy ; should not move
proton.injection_style = "singleparticle"
proton.single_particle_pos = 0. 0. 0.
proton.single_particle_u = -2. 0. 0.
proton.single_particle_weight = 1.

# Order of particle shape factors
algo.particle_shape = 1

# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 200
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz jx jy jz