# algo
algo.load_balance_intervals = 1:1:1,100
# amr.restart = diags/checkpoint00005462
warpx.write_diagnostics_on_restart
algo.maxwell_solver = Yee
algo.particle_shape = 3

# amr
amr.blocking_factor = 16
amr.max_grid_size = 4000
amr.max_grid_size_x = 1000
amr.max_grid_size_y = 4000
amr.max_level = 0
amr.n_cell = 16000 20000

# boundary
boundary.field_hi = pml pml
boundary.field_lo = pml pml
boundary.particle_hi = absorbing absorbing
boundary.particle_lo = absorbing absorbing

# c_HH
c_HH.species = hydrogen hydrogen
c_HH.type = pairwisecoulomb

# c_eH
c_eH.species = electrons hydrogen
c_eH.type = pairwisecoulomb

# c_ee
c_ee.species = electrons electrons
c_ee.type = pairwisecoulomb

# checkpoint
checkpoint.diag_type = Full
checkpoint.file_min_digits = 8
checkpoint.format = checkpoint
checkpoint.intervals = ::53263,8053:8053:,8053:8053:

# collisions
collisions.collision_names = c_ee c_eH c_HH

# diagnostics
diagnostics.diags_names = particles fields checkpoint

# electrons
electrons.charge = -q_e
electrons.density_function(x,y,z) = "n0*((x*x+z*z)<=(r0*r0))"
electrons.initialize_self_fields = 0
electrons.injection_style = nrandompercell
electrons.mass = m_e
electrons.momentum_distribution_type = constant
electrons.num_particles_per_cell = 128
electrons.profile = parse_density_function
electrons.ux = 0.0
electrons.uy = 0.0
electrons.uz = 0.0

# fields
fields.diag_type = Full
fields.fields_to_plot = By Ex Ez
fields.file_min_digits = 8
fields.format = openpmd
fields.intervals = ::26631,4324:20303:533,:4324:5326,9650:9650:,20303::5326,5656:17640:533,9650:9650:,17640:36282:5326
fields.openpmd_backend = bp
fields.particle_fields.average_gamma(x,y,z,ux,uy,uz) = sqrt(1 + ux**2+uy**2+uz**2)
fields.particle_fields.average_gamma.do_average = 1
fields.particle_fields.average_momentum_x(x,y,z,ux,uy,uz) = ux
fields.particle_fields.average_momentum_x.do_average = 1
fields.particle_fields.average_momentum_z(x,y,z,ux,uy,uz) = uz
fields.particle_fields.average_momentum_z.do_average = 1
fields.particle_fields.particle_number(x,y,z,ux,uy,uz) = 1
fields.particle_fields.particle_number.do_average = 0
fields.particle_fields_to_plot = average_gamma average_momentum_x average_momentum_z particle_number
fields.write_species = 0

# geometry
geometry.dims = 2
geometry.prob_hi = 6.4e-05 0.000128
geometry.prob_lo = -6.4e-05 -3.2e-05

# hydrogen
hydrogen.charge = q_e
hydrogen.density_function(x,y,z) = "n0*((x*x+z*z)<=(r0*r0))"
hydrogen.initialize_self_fields = 0
hydrogen.injection_style = nrandompercell
hydrogen.mass = m_p
hydrogen.momentum_distribution_type = constant
hydrogen.num_particles_per_cell = 128
hydrogen.profile = parse_density_function
hydrogen.ux = 0.0
hydrogen.uy = 0.0
hydrogen.uz = 0.0

# lasers
lasers.names = prepulse main_pulse

# main_pulse
main_pulse.direction = 0 0 1
main_pulse.do_continuous_injection = 0
main_pulse.e_max = 200668818412505.94
main_pulse.polarization = 1 0 0
main_pulse.position = 0 0 -3.1834e-05
main_pulse.profile = Gaussian
main_pulse.profile_duration = 1.2739827004487725e-14
main_pulse.profile_focal_distance = 3.1834e-05
main_pulse.profile_t_peak = 7.499999999999998e-14
main_pulse.profile_waist = 2.208236722976703e-06
main_pulse.wavelength = 8e-07

# my_constants
my_constants.n0 = 5.226e+28
my_constants.r0 = 2.2e-06

# particles
particles.diag_type = Full
particles.fields_to_plot = none
particles.file_min_digits = 8
particles.format = openpmd
particles.intervals = ::106525,8053:11248:533,9650:9650:,8053:11248:533,9650:9650:
particles.openpmd_backend = bp
particles.species = electrons hydrogen
particles.species_names = electrons hydrogen
particles.write_species = 1

# prepulse
prepulse.direction = 0 0 1
prepulse.do_continuous_injection = 0
prepulse.e_max = 3525349801870.9043
prepulse.polarization = 1 0 0
prepulse.position = 0 0 -3.1834e-05
prepulse.profile = Gaussian
prepulse.profile_duration = 1.2739827004487725e-14
prepulse.profile_focal_distance = 3.1834e-05
prepulse.profile_t_peak = 7.499999999999998e-14
prepulse.profile_waist = 2.123304541323753e-05
prepulse.wavelength = 8e-07

# stop_time
stop_time = 1.1811867940653797e-12

# warpx
warpx.break_signals = HUP
warpx.checkpoint_signals = USR1
warpx.const_dt = 1.8774888520527483e-17
warpx.sort_bin_size = 1 1 1
warpx.sort_intervals = 1
warpx.sort_particles_for_deposition = 0
warpx.verbose = 1

# Non PICMI options

# Reduced diagniostics:
warpx.reduced_diags_names = field_energy particle_energy particle_number laser_focus E_Hist_e_all E_Hist_e_ypos E_Hist_H_all E_Hist_H_ypos ypy_e ypy_H
# field_energy
field_energy.type = FieldEnergy
field_energy.intervals = 100
# particle_energy
particle_energy.type = ParticleEnergy
particle_energy.intervals = 100
# particle_number
particle_number.type = ParticleNumber
particle_number.intervals = 100
# laser_focus
laser_focus.type = FieldProbe
laser_focus.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
laser_focus.probe_geometry = Point
laser_focus.x_probe = 0
laser_focus.y_probe = 0
laser_focus.z_probe = 0
# E_Hist_e_all
E_Hist_e_all.type = ParticleHistogram
E_Hist_e_all.species = electrons
E_Hist_e_all.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
E_Hist_e_all.histogram_function(t,x,y,z,ux,uy,uz) = "sqrt(1 + ux**2 + uy**2 + uz**2) -1"
E_Hist_e_all.bin_number = 10240
E_Hist_e_all.bin_min = 0
E_Hist_e_all.bin_max = 2000.0
# E_Hist_e_ypos
E_Hist_e_ypos.type = ParticleHistogram
E_Hist_e_ypos.species = electrons
E_Hist_e_ypos.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
E_Hist_e_ypos.histogram_function(t,x,y,z,ux,uy,uz) = "sqrt(1 + ux**2 + uy**2 + uz**2) -1"
E_Hist_e_ypos.bin_number = 10240
E_Hist_e_ypos.bin_min = 0
E_Hist_e_ypos.bin_max = 2000.0
E_Hist_e_ypos.filter_function(t,x,y,z,ux,uy,uz) = uy>0
# E_Hist_H_all
E_Hist_H_all.type = ParticleHistogram
E_Hist_H_all.species = hydrogen
E_Hist_H_all.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
E_Hist_H_all.histogram_function(t,x,y,z,ux,uy,uz) = "sqrt(1 + ux**2 + uy**2 + uz**2) -1"
E_Hist_H_all.bin_number = 10240
E_Hist_H_all.bin_min = 0
E_Hist_H_all.bin_max = 10
# E_Hist_H_ypos
E_Hist_H_ypos.type = ParticleHistogram
E_Hist_H_ypos.species = hydrogen
E_Hist_H_ypos.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
E_Hist_H_ypos.histogram_function(t,x,y,z,ux,uy,uz) = "sqrt(1 + ux**2 + uy**2 + uz**2) -1"
E_Hist_H_ypos.bin_number = 10240
E_Hist_H_ypos.bin_min = 0
E_Hist_H_ypos.bin_max = 10
E_Hist_H_ypos.filter_function(t,x,y,z,ux,uy,uz) = uy>0
# ypy_e
ypy_e.type = ParticleHistogram2D
ypy_e.species = electrons
ypy_e.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
ypy_e.file_min_digits = 8
ypy_e.histogram_function_abs(t,x,y,z,ux,uy,uz,w) = uy
ypy_e.histogram_function_ord(t,x,y,z,ux,uy,uz,w) = y
ypy_e.value_function(t,x,y,z,ux,uy,uz,w) = "w"
ypy_e.bin_number_abs = 10240
ypy_e.bin_number_ord = 10240
ypy_e.bin_min_abs = -1000.0
ypy_e.bin_min_ord = -3.2e-05
ypy_e.bin_max_abs = 1000.0
ypy_e.bin_max_ord = 0.000128
# ypy_H
ypy_H.type = ParticleHistogram2D
ypy_H.species = hydrogen
ypy_H.intervals = ::5326,4324:20303:53,:4324:533,9650:9650:,20303::533,5656:17640:53,9650:9650:,17640:36282:533
ypy_H.file_min_digits = 8
ypy_H.histogram_function_abs(t,x,y,z,ux,uy,uz,w) = uy
ypy_H.histogram_function_ord(t,x,y,z,ux,uy,uz,w) = y
ypy_H.value_function(t,x,y,z,ux,uy,uz,w) = "w"
ypy_H.bin_number_abs = 10240
ypy_H.bin_number_ord = 10240
ypy_H.bin_min_abs = -10
ypy_H.bin_min_ord = -3.2e-05
ypy_H.bin_max_abs = 10
ypy_H.bin_max_ord = 0.000128

# IO options:
# particles
particles.adios2_engine.type = bp5
particles.adios2_engine.parameters.AggregatorRatio = 1
particles.adios2_engine.parameters.BufferChunkSize = 2147381248
particles.adios2_operator.type = blosc
particles.adios2_operator.parameters.compressor = zstd
particles.adios2_operator.parameters.clevel = 3
particles.adios2_operator.parameters.threshold = 2048
particles.adios2_operator.parameters.nthreads = 10
particles.adios2_operator.parameters.doshuffle = BLOSC_BITSHUFFLE
# fields
fields.adios2_engine.type = bp5
fields.adios2_engine.parameters.AggregatorRatio = 1
fields.adios2_engine.parameters.BufferChunkSize = 2147381248
fields.adios2_operator.type = blosc
fields.adios2_operator.parameters.compressor = zstd
fields.adios2_operator.parameters.clevel = 3
fields.adios2_operator.parameters.threshold = 2048
fields.adios2_operator.parameters.nthreads = 10