Arrays for outputs, refl_events and traj_events, are given length (2*nevents+1), but function only assigns values until position nevents+1
Arrays for outputs, refl_events and traj_events, are given length (2*nevents+1), but function only assigns values until position nevents+1