[richards] Implement and use RichardsVelocityOutput
Richards works conceptually with numPhases == 2 but only one velocity should be computed. Therefore we need to overload the velocity output.
Richards works conceptually with numPhases == 2 but only one velocity should be computed. Therefore we need to overload the velocity output.