CATHY breaks due to positive surface pressure heads after DA update

BenjMy / pycathy_wrapper

A Python wrapper for CATHY (Catchment Hydrology) - [Camporese et al., 2010]

https://benjmy.github.io/pycathy_wrapper/

Boost Software License 1.0

6 stars 1 forks source link

CATHY breaks due to positive surface pressure heads after DA update #9

Closed BenjMy closed 2 years ago

BenjMy commented 2 years ago

Problem description

The simulation consists in an infiltration in a confined domain (no flux on all the side boundaries). During DA update, after 3hours of infiltration, the state update generates a positive small value of pressure head (1e-3). CATHY breaks (the initial conditions being wrong) and the DA stops.

Simulation parameters

Pressure head initial: -5
Constant and uniform atmbc: 5e-7
no input of ponding head initial conditions
analysis method: Data Assimilation

BenjMy commented 2 years ago

Running the same, but decreasing atmbc to 5e-8 solved the issue.

[ ] Implement a user warning for positive surface pressure head during DA update if no pounding is set.