Set water content as initial condition
@hbauser would like to set the water content as initial condition for the Richards solver. This is important for data assimilation, where we usually represent the state in terms of water content.
Add new option
initial.quantity = waterContent.
This requires passing the
RichardsInitialConditionFactory, in order to create a
MatricHeadAdapter which translates water content back into matric head. This adapter is implemented in !134.
How to test the implementation?
Add simple test case to