#################################
####### GENERAL PARAMETERS ######
#################################
max_step = 1
amr.n_cell = 64 32 32
amr.max_grid_size = 16
amr.blocking_factor = 16
amr.max_level = 0
geometry.dims = 3
geometry.prob_lo = -20.e-6   -20.e-6   -20.e-6    # physical domain
geometry.prob_hi =  20.e-6    20.e-6    20.e-6

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

#################################
############ NUMERICS ###########
#################################
warpx.serialize_initial_conditions = 1
warpx.verbose = 1
warpx.cfl = 1.0

# Order of particle shape factors
algo.particle_shape = 1

#################################
############ PLASMA #############
#################################
particles.species_names = electrons

electrons.species_type = electron
electrons.injection_style = "NUniformPerCell"
electrons.num_particles_per_cell_each_dim = 1 1 2
#electrons.profile = "predefined"
#electrons.predefined_profile_name = "parabolic_channel"
#electrons.predefined_profile_params = 0.0 0.02 0.297 0.003 40.e-6 1.7e23
#electrons.density = predefined #1.e25  # number of electrons per m^3
electrons.profile = "fromfile"
electrons.read_density_from_path = "diags/asymm/openpmd_000000.h5"
electrons.density_name = "rho"
electrons.momentum_distribution_type = "gaussian"
electrons.ux_th  = 0.01 # uth the std of the (unitless) momentum
electrons.uy_th  = 0.01 # uth the std of the (unitless) momentum
electrons.uz_th  = 0.01 # uth the std of the (unitless) momentum

# Diagnostics
diagnostics.diags_names = asym #chk
asym.intervals = 1
asym.diag_type = Full
asym.electrons.variables = ux uy uz w
#asym.fields_to_plot = rho
asym.fields_to_plot = rho
asym.format = openpmd
asym.openpmd_backend = h5
#asym.openpmd_encoding = g


#chk.intervals = 6
#chk.diag_type = Full
#chk.format = checkpoint