From 25129d3e0f27f65086e035834788cd58ec705795 Mon Sep 17 00:00:00 2001
From: unknown <olivier.delaigue@ANPI1430.antony.irstea.priv>
Date: Tue, 13 Jun 2017 11:05:42 +0200
Subject: [PATCH] v1.0.8.0 new dataset added (Param_Sets_GR4J). It contains the
generalist parameter sets for the GR4J model #4562
---
DESCRIPTION | 2 +-
data/Param_Sets_GR4J.rda | Bin 0 -> 584 bytes
man/Param_Sets_GR4J.Rd | 85 +++++++++++++++++++++++++++++++++++++++
3 files changed, 86 insertions(+), 1 deletion(-)
create mode 100644 data/Param_Sets_GR4J.rda
create mode 100644 man/Param_Sets_GR4J.Rd
diff --git a/DESCRIPTION b/DESCRIPTION
index 08f030d3..2ce59047 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
Package: airGR
Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.0.7.2
+Version: 1.0.8.0
Date: 2017-06-13
Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl")),
diff --git a/data/Param_Sets_GR4J.rda b/data/Param_Sets_GR4J.rda
new file mode 100644
index 0000000000000000000000000000000000000000..7cf3d37dd9f99fbd63962fca7f039bdbfd44872b
GIT binary patch
literal 584
zcmV-O0=NAiiwFP!0000023=H5XcIvcp6qOt#A-#wiy(sF&n0S8ylB~1C~dJKO(T+G
zA#Cf$D8_d4Q=xh&LIp*rUWAGVK@bn3*MdS84}wrZP!Pd`MJ-5sP*RIZHO|cUvd|1n
zzIpF`?|bj<>}H@V(-zOfg%FmIQVL6|qiw05h_FRWVQarzbPL(ES1x6H2HKAaA-8Id
z+MCqgLEgtXwO-TceW4dsYDvB5oCN0XAwGF0rx)d40*eR1AKwk$jlIC@@8Eyh2YWa+
zqt;!@PQvdWfj>V3f7u5X-XMN{74`gE+|MIm`5Nk8J&E`@_u=10ekBjym0rQi%vb$v
z$|V~Z*mEYo`U|8GKcJuLf|1QRZ#nN1;=u;Ce}ETkVcs{+d&RiS`m9r30S4<izo|~J
zo$G&J+fUHS{2JpDbb~%xGan^kP&Wi=ke^e!1<$GTA?*G=!2SH(Pj@pvKz@q(r)~mE
zUw~bFZyx8T`8_$xI-Oj9>>2KL@G<gsVLm#+`UlWIjT{G_J&bq)^HG}n8Gdf!;e6|+
zhxYw^G40gBbAEGOYrf~nwuyO?V-McPJA7w_?WloX5);~SndIWs&e_zGvx+_)z7G-|
zCXwvflx8@~Jcq@+zpVQQ@sJ<kexq7;af$@bWpg}g+8v?@_hGTW1^qZRWhhV6XBxfe
z7QB)=Y&IV!oL1&gcxnqzNj}L;d!<=r<bqpjs5G4MGj7>U<kY`_*F=vMM=vEBg*R$(
Web@DR{i}k1djAK3%8x8y1ONb*1Rb>i
literal 0
HcmV?d00001
diff --git a/man/Param_Sets_GR4J.Rd b/man/Param_Sets_GR4J.Rd
new file mode 100644
index 00000000..9be3fd22
--- /dev/null
+++ b/man/Param_Sets_GR4J.Rd
@@ -0,0 +1,85 @@
+\docType{data}
+\encoding{UTF-8}
+
+\name{Param_Sets_GR4J}
+
+\alias{Param_Sets_GR4J}
+
+
+\title{Generalist parameter sets for the GR4J model}
+
+\format{Compromise short-list containing
+\itemize{
+ \item {GR4J X1} {production store capacity [mm]}
+ \item {GR4J X2} {intercatchment exchange coefficient [mm/d]}
+ \item {GR4J X3} {routing store capacity [mm]}
+ \item {GR4J X4u} {unajusted unit hydrograph time constant [d]}
+}}
+
+\description{
+This dataset of parameter sets can be used for the grid-screening calibration procedure.
+Please pay attention to the fact that the given GR4J X4u variable does not correspond to the actual GR4J X4 parameter.
+As explained in Andréassian \emph{et al.} (2014; section 2.1), the given GR4J X4u value has to be adjusted (rescaled) using the catchment area (S) [km2] as follows: {X4 = X4u * S^0.3}.
+Please, see the example below.
+}
+
+\seealso{
+ \code{\link{RunModel_GR4J}}, \code{\link{Calibration_Michel}}.
+}
+
+\references{
+Andréassian, V., F. Bourgin, L. Oudin, T. Mathevet, C. Perrin, J. Lerat, L. Coron, L. Berthet (2014).
+ Seeking genericity in the selection of parameter sets: impact on hydrological model efficiency.
+ Water Resources Research, 50(10), 8356-8366, doi: 10.1002/2013WR014761.
+}
+
+
+\examples{
+library(airGR)
+
+## loading catchment data
+data(L0123001)
+
+## loading generalist parameter sets
+data(Param_Sets_GR4J)
+str(Param_Sets_GR4J)
+
+## generalist parameter sets with real GR4J X4
+Param_Sets_GR4J$X4 <- Param_Sets_GR4J$X4u * BasinInfo$BasinArea^0.3
+Param_Sets_GR4J$X4u <- NULL
+Param_Sets_GR4J <- as.matrix(Param_Sets_GR4J)
+
+## preparation of the InputsModel object
+InputsModel <- CreateInputsModel(FUN_MOD = RunModel_GR4J, DatesR = BasinObs$DatesR,
+ Precip = BasinObs$P, PotEvap = BasinObs$E)
+
+## run period selection
+Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%d/\%m/\%Y \%H:\%M")=="01/01/1990 00:00"),
+ which(format(BasinObs$DatesR, format = "\%d/\%m/\%Y \%H:\%M")=="31/12/1999 00:00"))
+
+## preparation of RunOptions object
+RunOptions <- CreateRunOptions(FUN_MOD = RunModel_GR4J,
+ InputsModel = InputsModel, IndPeriod_Run = Ind_Run)
+
+## calibration criterion: preparation of the InputsCrit object
+InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel,
+ RunOptions = RunOptions, Qobs = BasinObs$Qmm[Ind_Run])
+
+## preparation of CalibOptions object using the generalist parameter sets
+CalibOptions <- CreateCalibOptions(FUN_MOD = RunModel_GR4J, FUN_CALIB = Calibration_Michel,
+ StartParamList = Param_Sets_GR4J)
+
+## calibration
+OutputsCalib <- Calibration(InputsModel = InputsModel, RunOptions = RunOptions,
+ InputsCrit = InputsCrit, CalibOptions = CalibOptions,
+ FUN_MOD = RunModel_GR4J, FUN_CRIT = ErrorCrit_NSE,
+ FUN_CALIB = Calibration_Michel)
+
+## simulation
+Param <- OutputsCalib$ParamFinalR
+OutputsModel <- RunModel(InputsModel = InputsModel, RunOptions = RunOptions,
+ Param = Param, FUN = RunModel_GR4J)
+
+## results preview
+plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run])
+}
\ No newline at end of file
--
GitLab