Include growing plant roots as a sink in the Richards source-term function
Description
The goal of my project is to simulate plant growth depending on atmospheric forcing (precipitation or evaporation) using the Feddes root-water uptake model.
Proposal
Include a root hydraulic head depending root function into the Flow-source function. It is also possible to include more simple root configurations for future use.
What I want for my project are three connected functionalities of which only one calculates the sink term. It is called Root. The first of the other two uses the water uptake to calculate biomass production (used for root and shoot growth) called Biomass while the second uses the upper boundary condition and the shoot biomass to calculate the potential transpiration for the next time-step and is called Shoot.
Main questions:
- how to read in initial parameters and conditions
- how to access the boundary condition at the top boundary (for calculation of T_p)
- how to access h_m
variables to save every time-step
(possibly as arrays to enable multiple plants in the future)
- Biomass B
- potential Transpiration T_p
- Water-uptake W
- hm (is already being saved)
- (Sink-value (only for possible multiple plants))
Flow diagram of the model idea
graph TB;
subgraph model
C[Shoot]-->|Tp|A;
B[Richards]-->|hm|A;
B-->|water-uptake|D;
D[Biomass production]-->|Biomass|A
D-->|Biomass|C
A[Root]==>|Sink value|B;
end
subgraph external parameters
S[Soil parameter]-.->B
P[Plant initial conditions and parameter]-.->A
P-.->C
P-.->D
Bo[Boundary conditions]-.->B
Bo-.->|only top|C
Bo-.->|grid extensions|A
end
Disclaimer:
This is my first Issue so please inform me about what I information I should add or clarify.