# amr.restart = "/home/labs/malkalab/tata/sim_data/BeamDriven/004_PML_Test/checkpoints/005000/"

# amrex.abort_on_out_of_gpu_memory = 0

warpx.break_signals = "USR2"
warpx.checkpoint_signals = "USR2"


#################################
####### GENERAL PARAMETERS ######
#################################


max_step = 30000       # 180000    # for production, run for longer time, e.g. max_step = 1000


#amr.n_cell = 128 128 1024   # for production, run with finer mesh, e.g. amr.n_cell = 64 64 512
#amr.n_cell = 256 256 8192
amr.n_cell = 128 128 1792  #4608 #13440   # for production, run with finer mesh, e.g. amr.n_cell = 64 64 24576

geometry.dims = 3

amr.max_grid_size_x = 128   # maximum size of each AMReX box, used to decompose the domain
amr.max_grid_size_y = 128
amr.max_grid_size_z = 128

amr.blocking_factor_x = 128  # minimum size of each AMReX box, used to decompose the domain
amr.blocking_factor_y = 128
amr.blocking_factor_z = 128

geometry.prob_lo     = -20.e-6   -20.e-6   -55.e-6    # physical domain
geometry.prob_hi     =  20.e-6    20.e-6    41.e-6

# Mesh refinement
amr.max_level = 0 # Maximum level in hierarchy (1 might be unstable, >1 is not supported)
# warpx.fine_tag_lo = -5.e-6   -5.e-6   -50.e-6
# warpx.fine_tag_hi =  5.e-6    5.e-6   -30.e-6

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

# boundary.field_lo = pml pml pml
# boundary.field_hi = pml pml pml

warpx.pml_ncell = 10
warpx.do_pml_in_domain = 0

warpx.pml_has_particles = 0
warpx.do_pml_j_damping = 0

warpx.do_dive_cleaning = 0
warpx.do_divb_cleaning = 0

#################################
############ NUMERICS ###########
#################################
warpx.verbose = 1

warpx.use_filter = 1
warpx.cfl = 0.999  # if 1., the time step is set to its CFL limit

algo.particle_pusher =  "boris"
algo.maxwell_solver = "psatd"
# algo.maxwell_solver = "ckc"

psatd.nox = 16
psatd.noy = 16
psatd.noz = 16

# Order of particle shape factors
algo.particle_shape = 3
algo.current_deposition = "esirkepov"


##################################
####### MOVING WINDOW ############
##################################

warpx.do_moving_window = 1
warpx.moving_window_dir = z # Only z is supported for the moment
warpx.moving_window_v = 1. # 

warpx.start_moving_window_step = 0
warpx.end_moving_window_step = -1 # -1 to not stop moving window




warpx.do_dynamic_scheduling = 1 # for production, set this to 1 (default)
warpx.serialize_initial_conditions = 0         # for production, set this to 0 (default)



#################################
############ PLASMA #############
#################################


#################################
############ LASER  #############
#################################
lasers.names        = laser1
laser1.profile      = Gaussian
laser1.position     = 0. 0. 40.e-6       # This point is on the laser plane (does not move with the moving window)
laser1.direction    = 0. 0. 1.           # The plane normal direction
laser1.polarization = 0. 1. 0.           # The main polarization vector
laser1.a0        = 2.                    # Maximum amplitude of the laser field (in V/m) a0=1
laser1.profile_waist = 5.e-6            # The waist of the laser (in m)
laser1.profile_duration = 22.e-15        # The duration of the laser (in s)
laser1.profile_t_peak = 40.e-15          # Time at which the laser reaches its peak (in s)
laser1.profile_focal_distance = 0.e-6  # Focal distance from the antenna (in m)
laser1.wavelength = 0.8e-6               # The wavelength of the laser (in m)

laser1.do_continuous_injection = 1

laser1.phi2 = 00.e-30

# Full Diagnostics
diagnostics.diags_names = diag1 checkpoint
diag1.openpmd_backend =    h5
diag1.file_prefix = "/home/labs/malkalab/tata/sim_data/BeamDriven/004_PML_Test/"

diag1.format = "openpmd"
# diag1.intervals = 5000:900000:2500, 100100:100100:1
diag1.intervals = 2500:900000:2500, 100100:100100:1 ,180000:185000:1000 , 105000:125000:1000
diag1.diag_type = Full
diag1.fields_to_plot = Ex Ey Ez Bx By Bz  # jx jy jz

# Corsinging ratio for fields on grid
diag1.coarsening_ratio = 2 2 4
# electrons ions ionsN coarsing ratio for particles



checkpoint.file_prefix = "/home/labs/malkalab/tata/sim_data/BeamDriven/004_PML_Test/checkpoints/"

checkpoint.format = "checkpoint"
checkpoint.intervals = 100:100:1,0:900000:5000
checkpoint.diag_type = Full