begin:control
# cells in x direction
nx = 100000
# final time of simulation, allows interaction with target for approx. 40 laser periods
t_end = 500 * femto
# vacuum from x=0 to x=50 microns, vacuum from x=thickness to x=200 microns
x_min = 0.0
x_max = 200 * micron
# first momentum update at 167 fs, so the plasma doesn't expand before laser arrives
particle_tstart = 160 * femto
dlb_threshold = 2500
stdout_frequency = 10
end:control
begin:boundaries
bc_x_min = simple_laser
bc_x_max = simple_outflow
end:boundaries
begin:constant
# ppc number
part = 5
lambda0 = 1 * micron
omega = 2.0*pi*c/lambda0
theta = 0
crit_den = critical(omega)
# solid density al13
den_start = 777*crit_den
# SET THICKNESS
thickness = 50.05*micron
end:constant
begin:laser
# LINEAR POLARISATION SO INTENSITY SHOULD BE 2e23 Wcm^-2, 60 fs
boundary = x_min
intensity_w_cm2 = 2e23
lambda = lambda0 * cos(theta)
pol_angle = 0
phase = 0
t_profile = gauss(time,60.0e-15,36.0e-15)
t_start = 0.0
end:laser
begin:species
#electron
name=electron
charge=-1.0
mass=1.0
npart=nx*part
density=den_start
density=if(x lt 50*micron, 0.0, density(electron))
density=if(x gt thickness, 0.0, density(electron))
# temp chosen so we can resolve Debye length of 2 nm
temp_ev = 60000
identify:electron
end:species
begin:species
name = al13
charge = 13.0
mass = 26.98 * 1836
npart = nx*part
density=density(electron)/13
# low ion background temperature, should be small compared to the electron temp
temp_ev = 10
end:species
begin:species
name = photon
identify:photon
end:species
begin:qed
# radiation reaction included but ignores pair production
# also produces photons but doesn't track their positions
qed_table_location =TABLES
use_qed = T
produce_photons = T
photon_energy_min = 1e6 * ev
use_radiation_reaction = F
end:qed
begin:output
name = fields
file_prefix = fields
# Here starting approx. 2 laser periods early so we can see laser profile just before interaction
time_start = 160 * femto
dt_snapshot = 3.33 * femto
# Properties on grid
grid = always + single
ex = always + single
ey = always + single
ez = always + single
bx = always + single
by = always + single
bz = always + single
# Only need particles on for phase-space plots
particles = always + single
px = always + single
py = always + single
particle_weight = always + single
particle_energy = always + single
number_density = always + species + single + no_sum
particle_probes = never
# Need this on to get the gamma factor
ekbar = always + species + single + no_sum
#grid = always + single
end:output