ion_velocities = 'random' in cp.x now works as expected.
added new routine randvel to do the initialization, extracted from the old one with some modification to do the equivalent of randpos for the velocities
the old routine maybe can be modified (not sure what it wanted to do in the middle of the cpr loop with tcap = .true. )