lcodePy-team
diff --git a/‎diag-test.ipynb
+9-10 b/‎diag-test.ipynb
+9-10
diff --git a/‎lcode2dPy/push_solver.py
+1-1 b/‎lcode2dPy/push_solver.py
+1-1
@@ -95,7 +95,7 @@ def step_dt(self, plasma_particles, plasma_fields, beam_source, beam_drain, t, d
 
             # Every xi step diagnostics
             if diagnostics:
-                diagnostics.dxi(t + self.config.getfloat('time-step'), layer_idx, plasma_particles, plasma_fields, rho_beam, beam_slice_to_move)
+                diagnostics.dxi(t + self.config.getfloat('time-step'), layer_idx, plasma_particles, plasma_fields, rho_beam, stable_slice)
             else:
                 if layer_idx % 100 == 0:
                     print('xi={xi:.6f} Ez={Ez:e} N={N}'.format(xi=layer_idx * -self.xi_step_p, Ez=plasma_fields.E_z[0], N=steps))