import numpy as np

# reference frequency = total plasma frequency
Tsim = 250.
indx = 9; indy = 8
omce = 0.2

vt_c = 0.05; Tc = vt_c**2 # thermal electron temperature
nec = 1.0 # thermal electron density

Debye_length = vt_c
Dx = Debye_length; Dy = Debye_length
Lx = 2.**indx*Dx; Ly = 2.**indy*Dy

a0 = 1e-2 * omce # amplitude of triggering wave

Main(
    geometry = "2d3v",
    interpolation_order = 2,
    
    cell_length = [Dx, Dy],
    number_of_cells = [2.**indx, 2.**indy],

    number_of_patches = [64, 32],

    timestep = 0.025,
    sim_time = Tsim,

    bc_em_type_x = ["silver-muller", "silver-muller"],
    bc_em_type_y = ["silver-muller", "silver-muller"],

    random_seed = 0,
)

ExtField(
    field = "Bx",
    profile = omce,
)

Species(
    species_type = "eonc",
    initPosition_type = "regular",
    initMomentum_type = "maxwell-juettner",
    n_part_per_cell = 100,
    mass = 1.0,
    charge = -1.0,
    nb_density = nec,
    mean_velocity = [0., 0., 0.],
    temperature = [Tc, Tc, Tc],
    bc_part_type_xmin = "refl",
    bc_part_type_xmax = "refl",
    bc_part_type_ymin = "refl",
    bc_part_type_ymax = "refl",
)
Species(
    species_type = "ponc",
    initPosition_type = "regular",
    initMomentum_type = "cold",
    n_part_per_cell = 49,
    mass = 1836.0,
    charge = 1.0,
    nb_density = 1.0,
    mean_velocity = [0., 0., 0.],
    time_frozen = Tsim,
    bc_part_type_xmin = "none",
    bc_part_type_xmax = "none",
    bc_part_type_ymin = "none",
    bc_part_type_ymax = "none",
)



Laser(
    boxSide = "xmin",
    omega = 0.2*omce,
    chirp_profile = tconstant(start=0.),
    time_envelope = tconstant(start=0.),
    space_envelope = [a0, a0], # amplitude included
    phase = [0., np.pi/2.], # right-hand polarized
)
                     
DiagScalar(          
    every = 1       
)                    
                     
DiagFields(          
    every = 50,      
    fields = ["Bx", "By", "Bz", "Ex", "Ey", "Ez", "Rho_eonc"],
)                    
                     
DiagProbe(           
    every = 1,       
    pos = [0., Ly/2.],
    pos_first = [Lx, Ly/2.],
    number = [2**indx],
    fields = ["Bx", "By", "Bz", "Ex", "Ey", "Ez"],
)