Check Richards analytic derivatives

The new benchmark cases merged in !2432 (merged) do not convergence (failing Newton convergence) with the analytic Jacobian implementation. There is a test using the analytic Jacobian (Richards lens test) so most parts should be tested.

Possible differences in the benchmark case: extrusion factor, ...?