Add steady state as initial condition
When discussing initial conditions, it came the idea to include steady states as initial condition for initialization of the simulation (See #122 (comment 24093)).
Solve the Richards/Transport equation in a steady state and its result feed a initial condition for the transient case.
How to test the implementation?
We could test it against the ODE solvers in python as we do for the transient case.
See #122 (closed)