inputs_mr_2d

max_step = 300
amr.n_cell = 128 128

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

# Geometry
geometry.dims = 2
geometry.prob_lo     = -32.e-6    -32.e-6      # physical domain
geometry.prob_hi     =  32.e-6     32.e-6
warpx.ref_patch_function(x,y,z) = " ((x > -16.e-6)*(x<-10.e-6) + (x > 9.e-6)*(x<16.e-6))*(z>-8.e-6)*(z<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 = -12.e-6 0. 0.
electron.single_particle_u = 2. 0. 0.
electron.single_particle_weight = 1.

proton.charge = q_e
proton.mass = m_p
proton.injection_style = "singleparticle"
proton.single_particle_pos = -12.e-6 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 = 300
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz jx jy jz