# This is a 2D electrostatic simulation containing a circle embedded boundary at
# the center of the grid domain, pre-seeded with equal density of argon ions
# and electrons.

max_step = 11
warpx.const_dt = 3.99e-13
warpx.do_electrostatic = labframe
warpx.self_fields_required_precision = 1e-06
warpx.eb_implicit_function = -((x-0.00005)**2+(z-0.00005)**2-1e-05**2)
warpx.eb_potential(x,y,z,t) = -10
warpx.self_fields_absolute_tolerance = 0.02

algo.load_balance_intervals = 5
algo.load_balance_costs_update = timers
algo.load_balance_efficiency_ratio_threshold = 1.001
algo.load_balance_with_sfc = 0
algo.load_balance_knapsack_factor = 2

amr.n_cell = 128 128
amr.max_grid_size = 16
amr.max_level = 0
geometry.dims = 2
geometry.prob_lo = 0.0 0.0
geometry.prob_hi = 0.0001 0.0001
boundary.field_lo = pec periodic
boundary.field_hi = pec periodic
boundary.particle_lo = absorbing periodic
boundary.particle_hi = absorbing periodic
algo.particle_shape = 1

particles.species_names = electrons ar_ions
electrons.species_type = electron
electrons.injection_style = nuniformpercell
electrons.initialize_self_fields = 1
electrons.num_particles_per_cell_each_dim = 2 2
electrons.density = 10.2e16
electrons.profile = constant
electrons.momentum_distribution_type = maxwell_boltzmann
electrons.theta = (kb*30000/(m_e*clight^2))
electrons.save_particles_at_xlo = 1
electrons.save_particles_at_xhi = 1
electrons.save_particles_at_eb = 1

ar_ions.species_type = argon
ar_ions.charge = q_e
ar_ions.injection_style = nuniformpercell
ar_ions.initialize_self_fields = 1
ar_ions.num_particles_per_cell_each_dim = 2 2
ar_ions.profile = constant
ar_ions.density = 10.2e16
ar_ions.momentum_distribution_type = gaussian
ar_ions.ux_m = 0.0
ar_ions.uy_m = 0.0
ar_ions.uz_m = 0.0
ar_ions.ux_th = 2.6285641070031447e-06
ar_ions.uy_th = 2.6285641070031447e-06
ar_ions.uz_th = 2.6285641070031447e-06
ar_ions.save_particles_at_xlo = 1
ar_ions.save_particles_at_xhi = 1
ar_ions.save_particles_at_eb = 1

collisions.collision_names = coll_electrons coll_ar_ions
coll_electrons.type = background_mcc
coll_electrons.species = electrons
coll_electrons.background_density = 2.15e+23
coll_electrons.background_temperature = 450
coll_electrons.scattering_processes = excitation1 elastic ionization
coll_electrons.excitation1_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/excitation_1.dat
coll_electrons.excitation1_energy = 11.5
coll_electrons.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/electron_scattering.dat
coll_electrons.ionization_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ionization.dat
coll_electrons.ionization_energy = 15.7596112
coll_electrons.ionization_species = ar_ions
coll_ar_ions.type = background_mcc
coll_ar_ions.species = ar_ions
coll_ar_ions.background_density = 2.15e+23
coll_ar_ions.background_temperature = 450
coll_ar_ions.scattering_processes = back elastic charge_exchange
coll_ar_ions.back_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ion_back_scatter.dat
coll_ar_ions.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ion_scattering.dat
coll_ar_ions.charge_exchange_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/charge_exchange.dat

diagnostics.diags_names = diag1 diag3
diag1.diag_type = Full
diag1.format = plotfile
diag1.intervals = 1
diag1.fields_to_plot = phi rho_electrons rho_ar_ions

diag3.diag_type = BoundaryScraping
diag3.format = openpmd
diag3.intervals = 5