Miracle error #874
Description
That is what I call devilry.
Here you see two images of Ey field taken 1/2 laser period apart at 40fs and 40.66fs (laser period of 0.4 micron laser is 1.3333 fs) . This is just the top part of the simulation box. Laser falls from the top. The color values of strips should flip on opposite and where was red (positive value) should be blue (negative one) and vice versa but they are exactly the same!
How this miracle might happen ? The fields solver is not running every time step?
# Focusing is set below but it is not used
begin:control
nx = 1000
ny = 46
nz = 46
t_end = 42 * femto
x_min = 0. * micron
x_max = 10.* micron
y_min = -0.23 * micron
y_max = 0.23 * micron
z_min = -0.23 * micron
z_max = 0.23 * micron
stdout_frequency = 100
nprocx=1
end:control
begin:boundaries
bc_x_min = simple_laser
bc_x_max = open
bc_y_min = periodic
bc_y_max = periodic
bc_z_min = periodic
bc_z_max = periodic
end:boundaries
begin:constant
lambda0 = 0.4 * micron
laser_fwhm2 = 8. * micron
laser_dur_fwhm = 50.0e-15
laser_angle = 0 # pi / 4
laser_period_s1 = 1.33333e-15
w_0 = 2.5479 * micron
x_spot = 9. * micron
x_R = pi * w_0^2 / lambda0 # Rayleigh range
RC = x_spot * (1.0 + (x_R/x_spot)^2) # Radius of curvature on x_min
w_bnd = w_0 * sqrt( 1.0 + (x_spot/x_R)^2) # Spot size at x_min
gouy = atan(x_spot/x_R) # Gouy phase shift at x_min
omega = 2 * pi * c / lambda0
n_crit = critical(omega)
den = 10. * n_crit
Hcharge = 1
sub_top = 5.* micron
sub_dep = 0.5*micron
r = sqrt(y^2 + z^2)
end:constant
begin:species
name = Hydrogen
charge = Hcharge
mass = 1836.2
number_density = if( (x gt sub_top) and (x lt sub_top+sub_dep), den, 0. )
temp_ev = 10.
npart_per_cell = 50
end:species
begin:species
name = electron
charge = -1.
mass = 1.
temp_ev = 10.
number_density = number_density(Hydrogen)
npart_per_cell = 50
end:species
begin:laser
boundary = x_min
lambda = lambda0
phase = 0
intensity_w_cm2 = 1e18
profile = gauss ( r, 0, laser_fwhm2 / (sqrt(2.*loge(2.0))))
t_profile = gauss ( time, 100.*femto, laser_dur_fwhm / (sqrt(2.*loge(2.0))))
end:laser
begin:output
name = normal0
file_prefix = m
time_start = 40.e-15
time_stop = 41.333333e-15
# time_stop = time_start + laser_period_s1 # 41.333333e-15
dt_snapshot = laser_period_s1 / 2.
grid = always
ex = always
ey = always
ez = always
bx = always
by = always
bz = always
end:output