Updated solveVelocity.f