phydrus.model.Model.add_waterflow¶
-
Model.add_waterflow(model=0, maxit=10, tolth=0.001, tolh=1, ha=1e-06, hb=10000.0, linitw=False, top_bc=0, bot_bc=0, hseep=0, rtop=None, rbot=None, rroot=None, gw_level=None, aqh=None, bqh=None, hysteresis=0, ikappa=- 1)[source]¶ Method to add a water_flow module to the model.
- Parameters
model (int, optional) – Soil hydraulic properties model: 0 = van Genuchten”s [1980] model with 6 parameters. 1 = modified van Genuchten”s model with 10 parameters [Vogel and Císlerová, 1988]. 2 = Brooks and Corey”s [1964] model with 6 parameters. 3 = van Genuchten”s [1980] model with air-entry value of -2 cm and with 6 parameters. 4 = Kosugi’s [1996] model with 6 parameters. 5 = dual-porosity model of Durner [1994] with 9 parameters. 6 = dual-porosity system with transfer proportional to the effective saturation (9 parameters). 7 = dual-porosity system with transfer proportional to the pressure head (11 parameters). 9 = dual-permeability system with transfer proportional to the pressure head (17 parameters) !!! model>3 options are not available with the major ion chemistry. module.
maxit (int, optional) – Maximum number of iterations allowed during any time step.
tolth (float, optional) – Absolute water content tolerance for nodes in the unsaturated part of the flow region [-]. TolTh represents the maximum desired absolute change in the value of the water content, θ, between two successive iterations during a particular time step.
tolh (float, optional) – Absolute pressure head tolerance for nodes in the saturated part of the flow region [L] (its recommended value is 0.1 cm). TolH represents the maximum desired absolute change in the value of the pressure head, h, between two successive iterations during a particular time step.
ha (float, optional) – Absolute value of the upper limit [L] of the pressure head interval below which a table of hydraulic properties will be generated internally for each material.
hb (float, optional) – Absolute value of the lower limit [L] of the pressure head interval for which a table of hydraulic properties will be generated internally for each material.
linitw (bool, optional) – Set to True if the initial condition is given in terms of the water content. Set to False if given in terms of pressure head.
top_bc (int, optional) – Upper Boundary Condition: 0 = Constant Pressure Head. 1 = Constant Flux. 2 = Atmospheric Boundary Condition with Surface Layer. 3 = Atmospheric Boundary Condition with Surface Run Off. 4 = Variable Pressure Head. 5 = Variable Pressure Head/Flux.
bot_bc (int, optional) – Lower Boundary Condition: 0 = Constant Pressure Head. 1 = Constant Flux. 2 = Variable Pressure Head. 3 = Variable Flux. 4 = Free Drainage. 5 = Deep Drainage. 6 = Seepage Face. 7 = Horizontal Drains.
hseep (float, optional) – Pressure head (i.e., 0) that initiates flow over the seepage face bottom boundary.
rtop (float, optional) – Prescribed top flux [LT-1] (in case of a Dirichlet BC set this variable equal to zero).
rbot (float, optional) – Prescribed bottom flux [LT-1] (in case of a Dirichlet BC set this variable equal to zero).
rroot (float, optional) – Prescribed potential transpiration rate [LT-1] (if no transpiration occurs or if transpiration is variable in time set this variable equal to zero).
gw_level (float, optional) – Reference position of the groundwater table (e.g., the x-coordinate of the soil surface).
aqh (float, optional) – Value of the parameter Aqh [LT-1] in the q(GWL)-relationship.
bqh (float, optional) – Value of the parameter Bqh [L-1] in the q(GWL)-relationship.
hysteresis (int, optional) – Hysteresis in the soil hydraulic properties: 0 = No hysteresis. 1 = Hysteresis in the retention curve only. 2 = Hysteresis in both the retention and hydraulic conductivity functions. 3 = Hysteresis using Robert Lenhard’s model [Lenhard et al., 1991; Lenhard and Parker, 1992]. (Not available with major ion chemistry module.)
ikappa (int, optional) – Set to -1 if the initial condition is to be calculated from the main drying branch. Set to 1 if the initial condition is to be calculated from the main wetting branch.