# Iterations
max_step = 400

# Domain
amr.n_cell =  32 32 32
amr.max_grid_size = 16
amr.max_level = 0

# Geometry
geometry.dims = 3
geometry.prob_lo = -50.e-6 -50.e-6 -50.e-6
geometry.prob_hi =  50.e-6  50.e-6  50.e-6

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

# Order of particle shape factors
algo.particle_shape = 3

# Numerics
warpx.cfl = 0.25
warpx.do_divb_cleaning = 1
warpx.use_filter = 1

# External magnetic field
my_constants.qm = 1e-1
warpx.B_ext_grid_init_style = parse_B_ext_grid_function
warpx.Bx_external_grid_function(x,y,z) = qm * x / (x*x + y*y + z*z)
warpx.By_external_grid_function(x,y,z) = qm * y / (x*x + y*y + z*z)
warpx.Bz_external_grid_function(x,y,z) = qm * z / (x*x + y*y + z*z)

# Particle beam
particles.species_names = beam
beam.charge = -q_e
beam.mass = 1.e30
beam.injection_style = "gaussian_beam"
beam.x_rms = 2.e-6
beam.y_rms = 2.e-6
beam.z_rms = 2.e-6
beam.x_m = 0.
beam.y_m = 0.
beam.z_m = 0.e-6
beam.npart = 20000
beam.q_tot = -1.e-20
beam.momentum_distribution_type = "constant"
beam.ux = 0.0
beam.uy = 0.0
beam.uz = 0.5

# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 0, 398:400:1
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz divE divB G