diff options
Diffstat (limited to 'lpic/input')
28 files changed, 7270 insertions, 0 deletions
diff --git a/lpic/input/input.2omegap.1 b/lpic/input/input.2omegap.1 new file mode 100644 index 0000000..f695722 --- /dev/null +++ b/lpic/input/input.2omegap.1 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 1.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 50.1317 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 10 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.9 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 50.1317 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 10 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 40 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 8000 # total number of cells +cells_left = 3501 # cells vacuum left +cells_plasma = 1000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 27.5625 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.01 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 51400 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = /tmp/u/rbp/data/data_2omgeap_1 # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 200 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 6 # start time in periods + t_stop = 16 # stop time in periods + x_start = 3000 # left boundary in cells + x_stop = 5000 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 59 # start time in periods + t_stop = 61 # stop time in periods + x_start = 1000 # left boundary in cells + x_stop = 2500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 1 # ex plots? + t_start = 6 # start time in periods + t_stop = 16 # stop time in periods + x_start = 3000 # left boundary in cells + x_stop = 5000 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 200 # stop time in periods + traces = 6 # # of traces at fixed positions x: + t0=2, t1=3350, t2=4000, t3=4450, t4=7990, t5=7998 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.2omegap.2 b/lpic/input/input.2omegap.2 new file mode 100644 index 0000000..298bfe0 --- /dev/null +++ b/lpic/input/input.2omegap.2 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 1.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 55.1501 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 10 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.9 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 55.1501 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 10 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 40 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 8000 # total number of cells +cells_left = 3501 # cells vacuum left +cells_plasma = 1000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 27.5625 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.01 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 51400 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = /tmp/u/rbp/data/data_2omgeap_2 # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 200 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 6 # start time in periods + t_stop = 16 # stop time in periods + x_start = 3000 # left boundary in cells + x_stop = 5000 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 59 # start time in periods + t_stop = 61 # stop time in periods + x_start = 1000 # left boundary in cells + x_stop = 2500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 1 # ex plots? + t_start = 6 # start time in periods + t_stop = 16 # stop time in periods + x_start = 3000 # left boundary in cells + x_stop = 5000 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 200 # stop time in periods + traces = 6 # # of traces at fixed positions x: + t0=2, t1=3350, t2=4000, t3=4450, t4=7990, t5=7998 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.diagnostics b/lpic/input/input.diagnostics new file mode 100644 index 0000000..8d2ec16 --- /dev/null +++ b/lpic/input/input.diagnostics @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 1 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 100 # cells per wavelength, lab frame +cells = 500 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 100 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 6 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 2000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 1 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&ion_phasespace + Q = 1 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&el_velocity + Q = 1 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&ion_velocity + Q = 1 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 1 # ion density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 1 # jx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 1 # jy plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 1 # jz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 1 # ex plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 1 # ey plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 1 # ez plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 1 # bx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 1 # by plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 1 # bz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 1 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=150, t2=200, t3=250, t4=498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.0 b/lpic/input/input.fresnel.0 new file mode 100644 index 0000000..5442ecd --- /dev/null +++ b/lpic/input/input.fresnel.0 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 1700 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 1500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=850, t3=1600, t4=1698 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.25p b/lpic/input/input.fresnel.25p new file mode 100644 index 0000000..860bade --- /dev/null +++ b/lpic/input/input.fresnel.25p @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 24.62 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 2200 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 2000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=1100, t3=2100, t4=2198 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.25p_fine b/lpic/input/input.fresnel.25p_fine new file mode 100644 index 0000000..b4c7acd --- /dev/null +++ b/lpic/input/input.fresnel.25p_fine @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 24.62 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 200 # cells per wavelength, lab frame +cells = 8800 # total number of cells +cells_left = 400 # cells vacuum left +cells_plasma = 8000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=400, t2=4400, t3=8400, t4=8798 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.25s b/lpic/input/input.fresnel.25s new file mode 100644 index 0000000..0dd5b0a --- /dev/null +++ b/lpic/input/input.fresnel.25s @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 24.62 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 2200 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 2000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=1100, t3=2100, t4=2198 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.30p b/lpic/input/input.fresnel.30p new file mode 100644 index 0000000..8816604 --- /dev/null +++ b/lpic/input/input.fresnel.30p @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 30.4503 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 2200 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 2000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=1100, t3=2100, t4=2198 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.30p_fine b/lpic/input/input.fresnel.30p_fine new file mode 100644 index 0000000..9cf19c3 --- /dev/null +++ b/lpic/input/input.fresnel.30p_fine @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 30.4503 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 200 # cells per wavelength, lab frame +cells = 8800 # total number of cells +cells_left = 400 # cells vacuum left +cells_plasma = 8000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=400, t2=4400, t3=8400, t4=8798 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.30s b/lpic/input/input.fresnel.30s new file mode 100644 index 0000000..0f9d962 --- /dev/null +++ b/lpic/input/input.fresnel.30s @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 30.4503 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 2200 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 2000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=1100, t3=2100, t4=2198 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.fresnel.37s b/lpic/input/input.fresnel.37s new file mode 100644 index 0000000..e72fc45 --- /dev/null +++ b/lpic/input/input.fresnel.37s @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude a_0 +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 37.472 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 50 # cells per wavelength, lab frame +cells = 2200 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 2000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0.64 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=100, t2=1100, t3=2100, t4=2198 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.1 b/lpic/input/input.harmonics.1 new file mode 100644 index 0000000..142107a --- /dev/null +++ b/lpic/input/input.harmonics.1 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 7 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.10 b/lpic/input/input.harmonics.10 new file mode 100644 index 0000000..79b6b6e --- /dev/null +++ b/lpic/input/input.harmonics.10 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 6.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 45.233 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 5 # pulse raise/fall time in periods +duration = 30 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 50 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 1000 # cells per wavelength, lab frame +cells = 7000 # total number of cells +cells_left = 3000 # cells vacuum left +cells_plasma = 1000 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 100 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 200 # max. number of MacroParticles per cell +vtherm = 0.02 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 200 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 2500 # left boundary in cells + x_stop = 3500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 2500 # left boundary in cells + x_stop = 3500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=2500, t2=3000, t3=3500, t4=4000, t5=4500, t6=6998 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.2 b/lpic/input/input.harmonics.2 new file mode 100644 index 0000000..96b14df --- /dev/null +++ b/lpic/input/input.harmonics.2 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 4 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data2 # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.3 b/lpic/input/input.harmonics.3 new file mode 100644 index 0000000..4d27f70 --- /dev/null +++ b/lpic/input/input.harmonics.3 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 30.4503 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 4 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.4 b/lpic/input/input.harmonics.4 new file mode 100644 index 0000000..9a674c1 --- /dev/null +++ b/lpic/input/input.harmonics.4 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 30.4503 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 4 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.5 b/lpic/input/input.harmonics.5 new file mode 100644 index 0000000..32a1e9f --- /dev/null +++ b/lpic/input/input.harmonics.5 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 3 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 30.4503 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 13 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.6 b/lpic/input/input.harmonics.6 new file mode 100644 index 0000000..0f61931 --- /dev/null +++ b/lpic/input/input.harmonics.6 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 3 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 30.4503 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 13 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.7 b/lpic/input/input.harmonics.7 new file mode 100644 index 0000000..96fcb9b --- /dev/null +++ b/lpic/input/input.harmonics.7 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 45.233 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 27.5625 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.8 b/lpic/input/input.harmonics.8 new file mode 100644 index 0000000..2550bf6 --- /dev/null +++ b/lpic/input/input.harmonics.8 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 45.233 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 100 # cells in the linear ramp region +n_ion_over_nc = 27.5625 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.9 b/lpic/input/input.harmonics.9 new file mode 100644 index 0000000..094b537 --- /dev/null +++ b/lpic/input/input.harmonics.9 @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 1.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 45.233 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 1000 # cells per wavelength, lab frame +cells = 7000 # total number of cells +cells_left = 3000 # cells vacuum left +cells_plasma = 1000 # occupied cells for cells_ramp=0 +cells_ramp = 200 # cells in the linear ramp region +n_ion_over_nc = 100 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 200 # max. number of MacroParticles per cell +vtherm = 0.02 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 200 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 2500 # left boundary in cells + x_stop = 3500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=2500, t2=3000, t3=3500, t4=4000, t5=4500, t6=6998 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.harmonics.matching b/lpic/input/input.harmonics.matching new file mode 100644 index 0000000..6a437d7 --- /dev/null +++ b/lpic/input/input.harmonics.matching @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.5 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 30 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 500 # cells per wavelength, lab frame +cells = 3500 # total number of cells +cells_left = 1500 # cells vacuum left +cells_plasma = 500 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 4 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 100 # max. number of MacroParticles per cell +vtherm = 0.03 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = /uplas/rbp/local/data/data_harm_matching_alpha30.2 # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 11 # start time in periods + t_stop = 15 # stop time in periods + x_start = 1250 # left boundary in cells + x_stop = 1750 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 7 # # of traces at fixed positions x: + t0=2, t1=1250, t2=1500, t3=1750, t4=2000, t5=2250, t6=3498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.list b/lpic/input/input.list new file mode 100644 index 0000000..21424e5 --- /dev/null +++ b/lpic/input/input.list @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 1 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 100 # cells per wavelength, lab frame +cells = 500 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 100 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 6 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 2000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=150, t2=200, t3=250, t4=498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.mcc b/lpic/input/input.mcc new file mode 100644 index 0000000..96f4476 --- /dev/null +++ b/lpic/input/input.mcc @@ -0,0 +1,277 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpi +// +////////////////////////////////////////////////////////////////////////////////////////// + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.256 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 2 # stop time in periods + +&ionization +------------------------------------------------------------------------------------------ +ion_alpha = 1.32 # atom dependent parameters +ion_u = 15.39 # ionization potential= ion_u * z^ion_alpha +part_z_max = 8 # maximum ion charge +omega_laser = 1.2 # laser frequency in eV +Q_ion_current = 0 # ionization current? +Q_info = 0 # additional info in error file + +&mcc_package +------------------------------------------------------------------------------------------ +xsection_data_file = input/neon_xsection_data + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 10 # cells per wavelength, lab frame +cells = 100 # total number of cells +cells_left = 50 # cells vacuum left +cells_plasma = 1 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 200 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.5 # thermal velocity in units of C + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 8 # q/e +#zmax = +8 # maximum q/e +m = 400000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.1 # thermal velocity in units of C + +&output +------------------------------------------------------------------------------------------ +path = /tmp/u/rbp/data/data_mcc # output path + +&energy + Q = 0 # energy plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 0 # flux plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 0 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 0.5 # time step in periods + +&snapshot_ion + Q = 1 # snapshots? + t_start = 1 # start time in periods + t_stop = 5 # stop time in periods + t_step = 1 # time step in periods + +&snapshot_macro + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 1 # start time in periods + t_stop = 4 # stop time in periods + x_start = 1800 # left boundary in cells + x_stop = 3600 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 1 # start time in periods + t_stop = 4 # stop time in periods + x_start = 1800 # left boundary in cells + x_stop = 3600 # right boundary in cells + +&traces + Q = 0 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=150, t2=200, t3=250, t4=498 + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 0 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.performance b/lpic/input/input.performance new file mode 100644 index 0000000..21424e5 --- /dev/null +++ b/lpic/input/input.performance @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 1 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 100 # cells per wavelength, lab frame +cells = 500 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 100 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 6 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 2000 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=150, t2=200, t3=250, t4=498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.positron b/lpic/input/input.positron new file mode 100644 index 0000000..f2c5d00 --- /dev/null +++ b/lpic/input/input.positron @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.5 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 2 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 2 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 20 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 100 # cells per wavelength, lab frame +cells = 500 # total number of cells +cells_left = 100 # cells vacuum left +cells_plasma = 100 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 6 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.02 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 1 # m/m_e +ppc = 50 # max. number of MacroParticles per cell +vtherm = 0.02 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path for all data + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 0 # snapshots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 1 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 1 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + t_step = 5 # time step in periods + +&de + Q = 1 # electron density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 1 # ion density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 0 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 10 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 5 # # of traces at fixed positions x: + t0=2, t1=150, t2=200, t3=250, t4=498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 3 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.pulse b/lpic/input/input.pulse new file mode 100644 index 0000000..6e7cf87 --- /dev/null +++ b/lpic/input/input.pulse @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 3 # pulse raise/fall time in periods +duration = 6 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 30 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 20 # cells per wavelength, lab frame +cells = 400 # total number of cells +cells_left = 0 # cells vacuum left +cells_plasma = 0 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 0 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 1836 # m/m_e +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 1 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 3 # # of traces at fixed positions x: + t0=1, t1=200, t2=399 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// diff --git a/lpic/input/input.reflection b/lpic/input/input.reflection new file mode 100644 index 0000000..bd538c7 --- /dev/null +++ b/lpic/input/input.reflection @@ -0,0 +1,259 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.01 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] +phase3 = 0 # 3rd harmonic's phase [degree] + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) +amplitude = 0.0 # dimensionless laser field amplitude +amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic +amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic +phase2 = 0 # 2nd harmonic's phase with respect to fundamental +phase3 = 0 # 3rd harmonic's phase + +angle = 0 # in degree +polarization = 1 # s=1, p=2, c=3 +shape = 2 # linear=1, sin=2, sin^2=3 +raise = 0 # pulse raise/fall time in periods +duration = 60 # pulse duration in periods +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 20 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 100 # cells per wavelength, lab frame +cells = 500 # total number of cells +cells_left = 200 # cells vacuum left +cells_plasma = 200 # occupied cells for cells_ramp=0 +cells_ramp = 0 # cells in the linear ramp region +n_ion_over_nc = 5 # maximum density/critical density +box_save = 1 # save configuration? yes=1, no=0 + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 20 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +z = 1 # q/e +m = 50000 # m/m_e +ppc = 10 # max. number of MacroParticles per cell +vtherm = 0.00 # thermal velocity in units of C + + +&output +------------------------------------------------------------------------------------------ +path = ../data # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&flux + Q = 1 # flux plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0 # time step in periods + +&reflex + Q = 1 # reflectivity plot? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 0.5 # time step in periods + +&el_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + t_step = 2 # time step in periods + +&el_velocity + Q = 0 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&ion_velocity + Q = 0 # ion velocity distributions? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + t_step = 10 # time step in periods + +&de + Q = 0 # electron density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&di + Q = 0 # ion density plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ey + Q = 0 # ey plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&ez + Q = 0 # ez plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&by + Q = 0 # by plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&bz + Q = 0 # bz plots? + t_start = 10 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&edens + Q = 1 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 20 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 500 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 100 # stop time in periods + traces = 6 # # of traces at fixed positions x: + t0=2, t1=100, t2=200, t3=250, t4=300, t5=498 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&restart +------------------------------------------------------------------------------------------ +Q = 0 # restart from intermediate stage? +file = restart # start file +Q_save = 1 # save intermediate stages periodically? +file_save = restart # save file + + +¶llel +------------------------------------------------------------------------------------------ +N_domains = 1 # number of parallel processes +Q_reo = 1 # periodic reorganizations? +delta_reo = 1 # laser cycles between reo's + + +////////////////////////////////////////////////////////////////////////////////////////// + + +thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T +------------------------------------------------------------------------------------------ +1e-4 0.005 +3e-4 0.046 +5e-4 0.128 +1e-3 0.511 +4e-3 8.2 +1e-2 51 +2e-2 204 +3e-2 460 +4e-2 818 +4.42e-2 1000 +1e-1 5110 + +choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) + :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl + +////////////////////////////////////////////////////////////////////////////////////////// |