diff options
Diffstat (limited to 'postpost/inp.lpi')
| -rw-r--r-- | postpost/inp.lpi | 397 |
1 files changed, 397 insertions, 0 deletions
diff --git a/postpost/inp.lpi b/postpost/inp.lpi new file mode 100644 index 0000000..4426727 --- /dev/null +++ b/postpost/inp.lpi @@ -0,0 +1,397 @@ +////////////////////////////////////////////////////////////////////////////////////////// +// +// input for lpic +// +////////////////////////////////////////////////////////////////////////////////////////// + + +&pulse_front +------------------------------------------------------------------------------------------ +Q = 1 # switch ON (Q=1) OFF (Q=0) +angle = 45 # in degree +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.0001 # time step in periods + + +&pulse_front_part_0 +------------------------------------------------------------------------------------------ +amplitude = 0 # dimensionless laser field amplitude +polarization = 2 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in periods +delay = 0 # propagation delay in cycles +phase = 0 # propagation delay in cycles +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 1 # frequency in fundamental frequencies + + +&pulse_front_part_1 +------------------------------------------------------------------------------------------ +amplitude = 5 # dimensionless laser field amplitude, 2nd harmonic +polarization = 2 # s=1, p=2, c=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # 2nd harmonic's phase with respect to fundamental [degree] +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 2 # frequency in fundamental frequencies +delay = 0 # delay of 2nd harmonic envelope with respect to main envelope + + +&pulse_front_part_2 +------------------------------------------------------------------------------------------ +amplitude = 0.0 # dimensionless laser field amplitude, 3rd harmonic +polarization = 2 # s=1, p=2, c=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # 3rd harmonic's phase [degree] +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 3 # frequency in fundamental frequencies +delay = 0 # delay of 3rd harmonic envelope with respect to main envelope + + +&pulse_front_part_3 +------------------------------------------------------------------------------------------ +amplitude = 0.1 # dimensionless laser field amplitude, nth harmonic +polarization = 1 # s=1, p=2, c=3 +raise = 0.31954069 # pulse raise/fall time in periods +duration = 0.63908137 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # nth harmonic's phase [degree] +chirp = -4445.9584 # linear chirp (dw/dt) in 2 pi T^-2 +delay = 15 # delay of nth harmonic envelope with respect to main envelope +frequency = 64.516129 # frequency of nth harmonic in fundamental frequencies + + +&pulse_rear +------------------------------------------------------------------------------------------ +Q = 0 # switch ON (Q=1) OFF (Q=0) + +angle = 0 # in degree +pulse_save = 1 # save pulse shape? yes=1, no=0 +pulse_save_step = 0.02 # time step in periods + + +&pulse_rear_part_0 +------------------------------------------------------------------------------------------ +amplitude = 0.0 # dimensionless laser field amplitude +polarization = 1 # s=1, p=2, c=3 +shape = 3 # linear=1, sin=2, sin^2=3 +raise = 10 # pulse raise/fall time in periods +duration = 20 # pulse duration in periods +delay = 0 # propagation delay in cycles +phase = 0 # propagation delay in cycles +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 1 # frequency in fundamental frequencies + + +&pulse_rear_part_1 +------------------------------------------------------------------------------------------ +amplitude = 0.0 # dimensionless laser field amplitude, 2nd harmonic +polarization = 2 # s=1, p=2, c=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # 2nd harmonic's phase with respect to fundamental +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 2 # frequency in fundamental frequencies +delay = 0 # delay of 2nd harmonic envelope with respect to main envelope + + +&pulse_rear_part_2 +------------------------------------------------------------------------------------------ +amplitude = 0.0 # dimensionless laser field amplitude, 3rd harmonic +polarization = 2 # s=1, p=2, c=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # 3rd harmonic's phase +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +frequency = 3 # frequency in fundamental frequencies +delay = 0 # delay of 3rd harmonic envelope with respect to main envelope + + +&pulse_rear_part_3 +------------------------------------------------------------------------------------------ +amplitude = 0.0 # dimensionless laser field amplitude, nth harmonic +polarization = 2 # s=1, p=2, c=3 +raise = 15 # pulse raise/fall time in periods +duration = 30 # pulse duration in fundamental periods +shape = 3 # linear=1, sin=2, sin^2=3 +phase = 0 # nth harmonic's phase [degree] +chirp = 0 # linear chirp (dw/dt) in 2 pi T^-2 +delay = 0 # delay of nth harmonic envelope with respect to main envelope +frequency = 14.0 # frequency of nth harmonic in fundamental frequencies + + +&propagate +------------------------------------------------------------------------------------------ +prop_start = 0 # start time in periods +prop_stop = 40 # stop time in periods + + +&box +------------------------------------------------------------------------------------------ +cells_per_wl = 1000 # cells per wavelength, lab frame +cells = 16200 # total number of cells +cells_left = 3500 # cells vacuum left +cells_plasma = 12500 # occupied cells for cells_ramp=0 +cells_rampgluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the righter exp-ramp +n_ion_over_nc = 196 # maximum density/critical density +ux0 = 0.0 # initial momentum in x-direction +box_save = 1 # save configuration? yes=1, no=0 + + +&lramp +------------------------------------------------------------------------------------------ +length = 0 # linear ramp length (form=0) or righter scale length (form=1 or 2) +form = 1 # 0 = linear; 1 = exponential; 2=doubleexponential +cutoff = -50 # max ramplength (<0 <=> disabled) +gluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the righter exp-ramp +middens = 0.5 # density of the constant portion in the middle (form=3) in fractions of ni +length2 = 0 # lefter scale length (form=2) + + +&rramp +------------------------------------------------------------------------------------------ +length = 0 # cells in the linear ramp region (form=0) or scalelength +form = 0 # 0 = linear; 1 = exponential +length2 = 100 # righter scale length (form=2) +cutoff = -1 # max ramplength (<0 <=> disabled) +gluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the lefter exp-ramp +middens = 0.5 # density of the constant portion in the middle (form=3) in fractions of ni + + +&electrons +------------------------------------------------------------------------------------------ +fix = 0 # 0->not fixed 1->fixed +ppc = 200 # max. number of MacroParticles per cell +vtherm = 0.000 # thermal velocity in units of C + + +&ions +------------------------------------------------------------------------------------------ +fix = 1 # 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_Thomsonstreaking_1_1 # output path + +&energy + Q = 1 # energy plot? + t_start = 0 # start time in periods + t_stop = 60 # 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 = 60 # 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 = 60 # stop time in periods + t_step = 1 # time step in periods + +&snapshot + Q = 1 # snapshots? + t_start = 0 # start time in periods + t_stop = 60 # 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 = 60 # stop time in periods + t_step = 0.05 # time step in periods + +&ion_phasespace + Q = 0 # phasespace plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + t_step = 0.05 # time step in periods + +&el_velocity + Q = 1 # electron velocity distributions? + t_start = 0 # start time in periods + t_stop = 60 # 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 = 60 # 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 = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&di + Q = 1 # ion density plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&jx + Q = 0 # jx plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&jy + Q = 0 # jy plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&jz + Q = 0 # jz plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&ex + Q = 0 # ex plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&ey + Q = 1 # ey plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&ez + Q = 1 # ez plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&bx + Q = 0 # bx plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&by + Q = 1 # by plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&bz + Q = 1 # bz plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&etx +------------------------------------------------------------------------------------------ + Q = 0 # thermal x-energy plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + + +&ety +------------------------------------------------------------------------------------------ + Q = 0 # thermal y-energy plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + + +&etz +------------------------------------------------------------------------------------------ + Q = 0 # thermal z-energy plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + + +&edens + Q = 1 # field energy density plots? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + x_start = 0 # left boundary in cells + x_stop = 100000 # right boundary in cells + +&traces + Q = 1 # traces? + t_start = 0 # start time in periods + t_stop = 60 # stop time in periods + traces = 22 # of traces at fixed positions x: + t0 =2, t1=3400, t2=3410,t3=3420, t4=3430, t5=3440, t6=3450, t7=3460; t8=3470, t9=3480, t10=3490, t11=3500, t12=3510, t13=3520, t14=3530, t15=3540, t16=3550, t17=3560, t18=3570, t19=3580, t20=3590, t21=3600 + + +////////////////////////////////////////////////////////////////////////////////////////// + + +&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 + +////////////////////////////////////////////////////////////////////////////////////////// + +#source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 p-pol transm. 1 minp-pol transm. 1 maxp-pol transm. 1 +#source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 +#source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 +#source p-pol refl. 1 |