Set water content as initial condition
Description
@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.
Proposal
Add new option initial.quantity = waterContent
.
This requires passing the FlowParameters
to 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 test-initial-condition
.
Related issues
See #127.