]>
Logging level of the Richards solver.
info
trace, debug, info, warning, error, critical
Path to the directory where most of the outputs are stored.
DORiE will attempt to create it if it does not exist.
./
path
Base file name for VTK output.
string
Plot vertex based (``true``) or cell-centered (``false``)
data into VTK files. Prefer vertex over cell data for full-precision
output. System tests and plotting functions (``dorie plot``) require
cell-centered data.
true, false
true
Policy to write the data.
endOfRichardsStep
endOfRichardsStep, none
Plot VTK files with virtually refined grids. VTK only
supports bilinear triangulations and displays higher-order solutions
inappropriately. Use level 0 for order 1, and level N for order N.
For level > 0, the printed grid does not resemble the actual grid.
This parameter defaults to 0 if not given in the config file. Notice
that subsampling significantly increases the output file size!
int
0
Defines whether VTK files should be written as ASCII (``true``)
or binary (``false``). ASCII is easier to parse in case you want to write
your own post-processing, but takes a lot more space on your hard drive.
true, false
false
Path to the boundary condition file.
path
Type of spatial segmentation of the boundaries specified in the BC file
rectangularGrid
rectangularGrid
Choose type of boundary segmentation: rectangularGrid
Whether to interpolate between the boundary conditions
at different times linearly (``true``) or not at all (``false``). May require
different boundary condition files.
true, false
false
The data type representing the initial condition. Either an
HDF datafile (``data``), or analytic equations (``analytic``).
data, analytic
analytic
Choose initial condition type: data, analytic
The physical quantity represented by the initial condition
data.
matricHead
matricHead
Choose quantity represented: matricHead
Equation for the initial condition
equation [x,y,z,h]
-h
Path to the initial condition data file
(``data`` type only). DORiE currently only supports H5 files with
file extension ``.h5``.
path
Dataset to use as initial condition (``data`` type only).
string
Interpolation type used for the data (``data`` type only).
nearest, linear
linear
Starting time in seconds.
float
0
Ending time in seconds.
float
1E6
Minimum time step that is allowed before DORiE stops running.
with an error, in seconds.
float
0.1
Value of the first time step in seconds.
float
10
Largest allowed time step in seconds. Use this to control the
density of your output.
float
1E5
Minimum number of Newton iterations of the solver per
time step. At maxTimestep, the Newton solver is not allowed to calculate more
than this number of iterations.
int
1
Maximum number of Newton iterations of the solver per
time step. At minTimestep, the Newton solver is not allowed to calculate more
than this number of iterations.
int
12
Factor the current time step is multiplied with when increasing
the time step.
float > 1
1.5
Factor the current time step is multiplied with when decreasing
the time step.
float < 1
0.5
YAML file containing the parameter definitions.
path
richards_param.yml
Polynomial order of the DG method used. Setting this value
to 0 (zero) selects the finite volume (FV) solver (only compatible to
structured rectangular grids).
1
0, 1, 2, 3
Select '0' for the finite volume (FV) solver
Upwinding method for skeleton terms. Upwinding typically
increases numeric stability while reducing accuracy.
**semiUpwind:** Apply upwinding to conductivity factor (only).
**fullUpwind:** Apply upwinding to conductivity.
**Not recommended for DG**.
none, semiUpwind, fullUpwind
none
Choose upwinding type: 'none', 'semiUpwind', 'fullUpwind'
DG discretization method for skeleton terms.
**SIPG:** Symmetric Interior Penalty
**NIPG:** Non-Symmetric Interior Penalty
**OBB:** Oden, Babuska, Baumann: no penalty term
**IIP:** Incomplete Interior Penalty: no symmetry term
SIPG, NIPG, OBB, IIP
SIPG
Apply harmonic weighting to skeleton term contributions
in DG.
true, false
true
Penalty factor to be used in the Discontinuous Galerkin scheme
float
10
Apply the flux reconstruction method to the solved matric
head and obtain conservative gradients. It always computes
(internally) the local lifting independent whether the ``lifting``
keyword is active or not.
true
true, false
Check that flux reconstruction engine is creating conforming
normal fluxes up to ``jumpTol``. ProTip: Setting warnings together
with a very low tolerance will let you track the changes on the
quality of the flux reconstruction.
none
none, warn, error
Whenever ``checkJumps`` is activated, it check that flux
reconstruction engine is creating conforming normal fluxes up to ``jumpTol``.
float > 0
1E-10
Force the Newton solver to calculate at least one iteration,
even if the initial defect is below ``AbsoluteLimit``. This ensures
that dynamics cannot be 'skipped' at very small time steps.
true, false
false
Absolute error tolerance of the Newton solver.
Reduce this value to increase precision.
1E-11 < float < 1E-8
1E-10
Relative error tolerance of the Newton solver.
Reduce this value to increase precision.
float
1E-4
If the ratio between last and current calculation defect exceeds this value,
the linear operator matrix is reassembled to ensure convergence.
float
5E-2
Maximum iteration count of linear searches performed to deduce
the optimal damping factor for reducing the defect.
int
10
Required defect reduction of the linear solver. The Newton solver calculates
the required linear reduction for second order Newton convergence automatically and
chooses the smaller value of both.
float
1E-3