! 3D simulation of plasma wakefield diagnostic using photon acceleration shot in oblique angle.
! There is one laser probe pulses in this case using antenna.
! The driver of the wakefield is an electron beam.
! The probe pulse crosses the wakefield generated by the pump with angle (approx.) 20 degrees and 5 degrees.

! Plasma density: 1.8E19 cm^-3 (n/n0 = 1)
! Laser wavelength: 800 nm (w0/wp = 9.8436826981031498)
! Probe normalised intensity: 0.01

!----------the node configuration for this simulation----------
node_conf 
{
  node_number(1:3) = 32,2,2,
  if_periodic(1:3) = .false.,.true.,.true.,
}
!----------spatial grid----------
grid 
{
  nx_p(1:3) = 6500,400,300,
  coordinates = "cartesian",
}
!----------time step and global data dump timestep number----------
time_step 
{
  dt     =  0.008,
  ndump  =  625, 
}
!----------restart information----------
restart 
{
  ndump_fac = 8,
  if_restart=.false.,
  if_remold=.true.,
}
!----------spatial limits of the simulations----------
! (note that this includes information about
! the motion of the simulation box)
space 
{
  ! for x1 direction, dx = 0.01, to accommodate the laser wavelength
  ! for x2, dx = 0.1, to comply with a rule of thumb kp*dx = 0.1
  xmin(1:3) =  0.000d0,-20.000d0,-15.000d0,
  xmax(1:3) =  65.0d0,20.0d0,15.000d0,
  if_move(1:3)= .true.,.false.,.false.,
  move_u = 9.8436826981031498,
}
!----------time limits ----------
time 
{
  tmin =   0.0, 
  tmax =   120.0,
}
!----------boundary conditions for em-fields ----------
emf_bound 
{
  type(1:2,1) = 30, 30,
  type(1:2,2) = 60, 60,
  type(1:2,3) = 60, 60,
}
!----------diagnostic for electromagnetic fields---------- 
diag_emf 
{
  ndump_fac = 12,
  ndump_fac_lineout = 1,
  reports = "e3", "e3, slice, x2, 200", "e3, slice, x2, 170", "e3, slice, x2, 140", "e3, slice, x2, 110", "e3, slice, x2, 90",
}
!----------number of particle species----------
particles 
{ 
  num_species = 2,
}
!==================================================
! ELECTRONS
!==================================================
species 
{
  name="electrons",
  num_par_max = 30000000,
  rqm=-1.0,
  num_par_x(1:3) = 1,1,1,
  den_min = 1.0d-12,
  interpolation = "quadratic",
  vth(1:3) = 0.0d0, 0.0d0, 0.0d0,
  vfl(1:3)=  0.0d0, 0.0d0, 0.0d0,
  q_real = -1.0,
}
!----------density profile for this species----------
profile 
{
  profile_type = "math func",
  math_func_expr = "1 * if(x1 > 60.0, 1, 0)",
}
!----------boundary conditions for this species----------
spe_bound 
{
  type(1:2,1) =  5, 5,
  type(1:2,2) =  5, 5,
  type(1:2,3) =  5, 5,
}
!----------diagnostic for this species----------
diag_species 
{
  ndump_fac_ave = 12,
  ndump_fac_lineout = 1,
  n_ave(1:3) = 20, 2, 2,
  reports = "charge, savg", "charge, slice, x3, 150",
}

!==================================================
! E-DRIVER
!==================================================
species 
{
  name="edriver",
  num_par_max = 9000000,
  rqm=-1.0,
  num_par_x(1:3) = 1,1,1,
  den_min = 1.0d-12,
  interpolation = "linear",
  vth(1:3) = 0.0d0, 0.0d0, 0.0d0,
  vfl(1:3) = 42277.9d0, -15387.9d0, 0.0d0,
  q_real = -1.0,
}
!----------density profile for this species----------
profile 
{
  !profile_type = "math func",
  !math_func_expr = "exp(-((x1-55)*0.9396926207859084-(x2-10)*0.3420201433256687)^2/4.5 - ((x1-55)*0.3420201433256687+(x2-10)*0.939692607859084)^2/18 - x3^2/18)",
  profile_type(1:3) = "gaussian", "gaussian", "gaussian",
  gauss_center(1:3) = 55.0, 10.0, 0.0,
  gauss_sigma(1:3) = 3.5, 3.5, 3.5,
  gauss_range(1:2,1) = 45.0, 65.0, 
  gauss_range(1:2,2) = 0.0, 20.0,
  gauss_range(1:2,3) = -10.0, 10.0,
  density = 0.33,
}
!----------boundary conditions for this species----------
spe_bound 
{
  type(1:2,1) =  5, 5,
  type(1:2,2) =  5, 5,
  type(1:2,3) =  5, 5,
}
!----------diagnostic for this species----------
diag_species 
{
  !ndump_fac_lineout = 4,
  !reports = "charge, slice, x3, 0",
}
!=====================================
! IONS - not specified to make them immovable
!=====================================

!=====================================
! PROBE PULSE
!=====================================
zpulse
{
  a0 = 0.01,
  omega0 = 9.8436826981031498,
  pol = 90.0,
  direction = 1,
  lon_type = "gaussian",
  lon_duration = 13.3,
  lon_x0 = 25.0,
  lon_range = 50.0,
  per_type = "plane",
}

!-------------smooth for currents------------- 
smooth
{
  smooth_type = "compensated",
}

! ------------- electrical current settings ----------------
current{}
! ------------- electrical charge settings ----------------
charge{}
! ------------- end of osiris input file ---------------
