Code works for simulation argon. Possible bug for small timestep if particle moves more than one cell outside domain