View source: R/reservoir_strategic.R
reservoir_strategic | R Documentation |
db_fill()
Function generates the reservoir parameter file from a pre-processed reservoir vector map and optional supplemental parameter file.
reservoir_strategic(res_vect = NULL, subbasin = NULL,
res_file = NULL, dir_out = "./", reservoir_file = "reservoir.txt",
overwrite = F, silent = F)
res_vect |
Name of reservoir vector map in GRASS location. Should be point
instead of polygon feature (i.e. reservoir outlet locations; consider function
|
subbasin |
Subbasin raster map in GRASS location. Can be created with
|
res_file |
tab-delimited file containing reservoir properties (fields see details) and key |
dir_out |
Character string specifying output directory (will be created if it
does not yet exist). The result file can later be imported to the database with |
reservoir_file |
Output: File of parameters for the strategic reservoirs
assigned to subbasins. To be filled into a database using |
overwrite |
|
silent |
|
For each reservoir that should be modelled explicitly within WASA the
following information need to be collected and written into the vector file's
attribute table or res_file
. Column order is not important. Additional
columns can be given but will be ignored:
res_id
Unique numeric reservoir identifier (if res_file
is given, it also needs to be
defined in the vector file's attribute table!).
name
OPTIONAL: name of the reservoir. Will be filled with <NA>
if not given.
minlevel
Initial minimum level in the reservoir [m]. Value varies because of sediment
accumulation. Default: 0.
maxlevel
Maximum water level in the reservoir [m]. Default: Estimated using Molle's equation (alpha = 2.7, k = 1500).
vol0
Initial volume of the reservoir [10^3 m^3]. Value varies because of sediment
accumulation. Set to '-999' if information is not available. Default: Estimated using Molle's equation (alpha = 2.7, k = 1500).
storecap
Initial storage capacity of the reservoir [10^3 m^3]. Value varies because of
sediment accumulation. DEfault: vol0
damflow
Target outflow discharge of the reservoir (90 % reliability) [m^3/s]. Default: 9.99
damq_frac
Fraction of Q90 released from the reservoir in regular years [-]. Default: 1
withdrawal
Water withdrawal discharge from the reservoir to supply the water use sectors
[m^3/s]. Outflow discharge through the dam is not considered. Default: 0
damyear
Year of construction of the dam (YYYY). Default: 1900
maxdamarea
Initial maximum area of the reservoir [ha]. Value varies because of sediment
accumulation. Default from GIS
damdead
Initial dead volume of the reservoir [10^3 m^3]. Value varies because of
sediment accumulation. Default: 0.01 * storecap
damalert
Initial alert volume of the reservoir [10^3 m^3]. Value varies because of
sediment accumulation. Default: 0.01 * storecap
dama, damb
Parameters of the area-volume relationship in the reservoir:
area = dama * Vol^damb [-]. Values of reservoir area and volume are
expressed in m^2 and m^3, respectively. Attention: These are not the Molle parameters! The correct values can be derived from them as shown in the default.
Default: dama = (1/k * (alpha * k)^(alpha/(alpha-1)))^((alpha-1)/alpha),
damb = (alpha-1)/alpha
q_outlet
Maximum outflow discharge released through the bottom outlets of the
reservoir [m^3/s]. Default: 0.999
fvol_botm
Fraction of storage capacity that indicates the minimum storage volume for
sediment release through the bottom outlets of the reservoir [-]. Default: 0.01
fvol_over
Fraction of storage capacity that indicates the minimum storage volume for
water release through the spillway of the reservoir [-]. Default: 1
damc, damd
Parameters of the spillway rating curve of the reservoir: Qout = damc * Hv^damd
[-]. Values of water height over the spillway and overflow discharges are
expressed in m and m^3/s, respectively.
Default: 1.6 * 30, 1.5
elevbottom
Bottom outlet elevation of the reservoir [m]. Currently ignored in WASA. Default: 99
The output file reservoir_file
contains the additional column pid
which is the corresponding subbasin ID determined from input subbasin
.
If you applied reservoir_outlet
using the outlet locations
of strategic reservoirs as drainage points, it might be necessary to use
locations of the function's output points_processed
instead of the
true reservoir outlet locations as otherwise the reservoirs might get assigned to
the wrong subbasins!
Tobias Pilz tpilz@uni-potsdam.de
lumpR package introduction with literature study and sensitivity analysis:
Pilz, T.; Francke, T.; Bronstert, A. (2017): lumpR 2.0.0: an R package facilitating
landscape discretisation for hillslope-based hydrological models.
Geosci. Model Dev., 10, 3001-3023, doi: 10.5194/gmd-10-3001-2017
WASA model in general:
Guentner, A. (2002): Large-scale hydrological modelling in the semi-arid
North-East of Brazil. PIK Report 77, Potsdam Institute for Climate
Impact Research, Potsdam, Germany.
Reservoir module of the WASA model:
Mamede, G. L. (2008): Reservoir sedimentation in dryland catchments: Modeling
and management. PhD Thesis, University of Potsdam, Germany.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.