added parameter READ_VELOCITIES_AT_f0 and implemented a shift of velocities read from the Par_file if needed when viscoelasticity is on