Commit 93e07a11 authored by Delaigue Olivier's avatar Delaigue Olivier
Browse files

v1.3.0.8 DOC: RunModel_CemaNeigeGR4H example modified

Showing with 8 additions and 6 deletions
+8 -6
Package: airGR Package: airGR
Type: Package Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
Version: 1.3.0.7 Version: 1.3.0.8
Date: 2019-05-20 Date: 2019-05-20
Authors@R: c( Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")), person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
......
...@@ -14,7 +14,7 @@ output: ...@@ -14,7 +14,7 @@ output:
### 1.3.0.7 Release Notes (2019-05-20) ### 1.3.0.8 Release Notes (2019-05-20)
#### New features #### New features
......
...@@ -89,6 +89,9 @@ library(airGR) ...@@ -89,6 +89,9 @@ library(airGR)
## loading catchment data ## loading catchment data
data(L0123002) data(L0123002)
BasinObs$DatesR <- seq(from = as.POSIXct("1984-01-01", tz = "UTC"),
to = as.POSIXct("1984-01-01", tz = "UTC") + length(BasinObs$DatesR),
by = "hour")
## preparation of the InputsModel object ## preparation of the InputsModel object
InputsModel <- CreateInputsModel(FUN_MOD = RunModel_CemaNeigeGR4H, DatesR = BasinObs$DatesR, InputsModel <- CreateInputsModel(FUN_MOD = RunModel_CemaNeigeGR4H, DatesR = BasinObs$DatesR,
...@@ -97,8 +100,7 @@ InputsModel <- CreateInputsModel(FUN_MOD = RunModel_CemaNeigeGR4H, DatesR = Basi ...@@ -97,8 +100,7 @@ InputsModel <- CreateInputsModel(FUN_MOD = RunModel_CemaNeigeGR4H, DatesR = Basi
HypsoData = BasinInfo$HypsoData, NLayers = 5) HypsoData = BasinInfo$HypsoData, NLayers = 5)
## run period selection ## run period selection
Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-01-01"), Ind_Run <- 2193:5844
which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1999-12-31"))
## ---- original version of CemaNeige ## ---- original version of CemaNeige
...@@ -114,7 +116,7 @@ OutputsModel <- RunModel_CemaNeigeGR4H(InputsModel = InputsModel, ...@@ -114,7 +116,7 @@ OutputsModel <- RunModel_CemaNeigeGR4H(InputsModel = InputsModel,
RunOptions = RunOptions, Param = Param) RunOptions = RunOptions, Param = Param)
## results preview ## results preview
plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run]) #plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run])
## efficiency criterion: Nash-Sutcliffe Efficiency ## efficiency criterion: Nash-Sutcliffe Efficiency
InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel, InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel,
...@@ -135,7 +137,7 @@ OutputsModel <- RunModel_CemaNeigeGR4H(InputsModel = InputsModel, ...@@ -135,7 +137,7 @@ OutputsModel <- RunModel_CemaNeigeGR4H(InputsModel = InputsModel,
RunOptions = RunOptions, Param = Param) RunOptions = RunOptions, Param = Param)
## results preview ## results preview
plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run]) #plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run])
## efficiency criterion: Nash-Sutcliffe Efficiency ## efficiency criterion: Nash-Sutcliffe Efficiency
InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel, InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel,
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment