Support spatially varying material parameters

This should work in the end:

2d Richards with varying material parameters. An option here would be to take two sets of parameters as already available in the code and define a function s(x,z) which varies between 0,1 so that the material is e.g, sand if s=1 and clay if s=0. We should probably choose s to be continuous to avoid making the solver more challenging.