begin:constant
ndens=5.0e22
tm=2.5e4
tb=2.5e4
mr=12853.1
vte=sqrt(kb*tm/me)
wpe=sqrt(ndens*(qe)^2/(me*epsilon0))
dl=vte/wpe
tend=4400.0/wpe
end:constant

begin:control
  n_zeros=3
  nx=7992
  ny=24
  nz=24
  
  nprocx=12
  nprocy=4
  nprocz=4

dlb_threshold=0.9
  
  #maximum number of iterations
  #set to -1 to run until finished
  nsteps=-1

  #final time of simulation
  t_end=tend

  #size of domain
  x_min=0.0
  x_max=dl*nx*270.0
  
  y_min=0.0
  y_max=dl*ny*200.0
  
  z_min=0.0
  z_max=dl*nz*200.0
  
end:control


begin:boundaries

  bc_x_min = periodic
  bc_x_max = periodic
  bc_y_min = periodic
  bc_y_max = periodic
  bc_z_min = periodic
  bc_z_max = periodic
    
end:boundaries


begin:species
  name=ele
  charge=-1.0
  mass=1.0
  npart=nx*ny*nz*16.0
density=ndens*(1.0+(x)/x_max*9.0)*(tanh(x/(0.005*x_max))+tanh(-(x-x_max)/(0.005*x_max))-1.0)
temp=tm
end:species

begin:species
  name=ion
  charge=1.0
  mass=12853.1
  npart=nx*ny*nz*16.0
density=ndens*(1.0+(x)/x_max*9.0)*(tanh(x/(0.005*x_max))+tanh(-(x-x_max)/(0.005*x_max))-1.0)
temp=tm
end:species

begin:species
  name=ebt
  charge=-1.0
  mass=1.0
  npart=nx*ny*nz*16.0
density=ndens*exp(-((x-(10.0*(c/wpe)))^2)/((2.0*(2.0*c/wpe)^2)) )*exp(-((y-y_max/2.0)^2)/(2.0*(c/wpe)^2))*exp(-((z-z_max/2.0)^2)/(2.0*(c/wpe)^2))
temp=tb
drift_x=1.087587e-17
end:species

begin:species
  name=ebd
  charge=-1.0
  mass=1.0
  npart=nx*ny*nz*16.0
density=ndens*2.5*exp(-((x-(15.7*(c/wpe)))^2)/((2.0*(c/wpe)^2)) )*exp(-((y-y_max/2.0)^2)/(2.0*(c/wpe)^2))*exp(-((z-z_max/2.0)^2)/(2.0*(c/wpe)^2))
temp=tb
drift_x=1.087587e-17
end:species

begin:output

  #number of timesteps between output dumps
  dt_snapshot=(tend/10.0)
  restart_dump_every=-1
  
  #Properties on grid
  grid=always
  ex=always
  ey=always
  ez=always
  bx=always
  by=always
  bz=always
  jx=always
  jy=always
  jz=always
  ekbar=always
  total_energy_sum=always
  charge_density=always+species
  number_density=always+species
  temperature=always+species

  #extended io
  distribution_functions=always
end:output

begin:dist_fn
  name=dfpx
  ndims=1
  dumpmask=always

  direction1=dir_px

  range1=(-1.4e-17,1.4e-17)
  resolution1=100000000

  include_species:ele
  include_species:ion
  include_species:ebt
  include_species:ebd
  end:dist_fn

begin:fields
ex=0.0
ey=0.0
ez=0.0
bx=0.0
by=0.0
bz=0.0
end:fields