Revert to 1st order time integration.
Now using the velocity bounds given in the initialisation file to build the ion histogram