This fixes #216 (closed). The initial velocities are now initialized (again) in hydro_convert_quantities by using p->v, before they are used to set the momentum. I also removed the if condition from hydro_end_density, as we do not longer care if the primitive variables are overwritten after the initialization step.