added parameter READ_VELOCITIES_AT_f0 and implemented a shift of velocit...

...ies read from the Par_file if needed when viscoelasticity is on