Moved some of the global copies ou of the loop

I have also changed the time stepper to the velocity verlet method. It
is about 10% faster than the Beeman method and requires one less array.
I am thinking that it might be necessary to have such a high order
integrator for the newtonian motion.
I can't remember what has happened