Loading libraries, network and time series data

Run vignette("01_First_network", package = "airGRiwrm") and vignette("02_First_run", package = "airGRiwrm") before this one in order to create the Rdata files loaded below:

library(airGRiwrm)
load("_cache/V01.RData")
load("_cache/V02.RData")
library(seinebasin)
data(QNAT)

InputsCrit object

InputsCrit <- CreateInputsCrit(
  InputsModel = InputsModel,
  FUN_CRIT = airGR::ErrorCrit_KGE2,
  RunOptions = RunOptions, Obs = Qnat[IndPeriod_Run,]
)
str(InputsCrit)
## List of 25
##  $ TRANN_01:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0943 0.0832 0.1443 0.1387 0.1498 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ STDIZ_04:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.11 0.155 0.409 0.184 0.258 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ BAR-S_06:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.118 0.111 0.111 0.118 0.111 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ CHAUM_07:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.599 0.239 0.239 0.12 0.2 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ CUSSY_08:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0697 0.0348 0.0348 0.0348 0.0348 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ STGER_09:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.579 0.322 0.386 0.429 0.472 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ GUILL_10:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0177 0.0177 0.0177 0.0177 0.0177 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ AISY-_11:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.064 0.064 0.064 0.064 0.064 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ EPISY_14:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.179 0.165 0.159 0.159 0.159 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ MONTR_18:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.263 0.263 0.263 0.263 0.263 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ LOUVE_19:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.15 0.187 0.206 0.337 0.281 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ LASSI_20:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0591 0.0591 0.0986 0.0986 0.0986 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ VITRY_25:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.086 0.0737 0.0737 0.0737 0.0737 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ MERY-_22:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0886 0.0975 0.0886 0.0886 0.082 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ GURGY_02:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.21 0.219 0.285 0.271 0.256 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ CHABL_12:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0929 0.0929 0.0155 0.0232 0.0155 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ BRIEN_03:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.113 0.104 0.104 0.113 0.104 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ ARCIS_24:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0697 0.0529 0.0649 0.0601 0.0961 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ CHALO_21:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.121 0.157 0.135 0.159 0.264 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ COURL_23:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.205 0.196 0.2 0.218 0.209 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ NOISI_17:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.139 0.146 0.163 0.17 0.176 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ NOGEN_13:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.0762 0.0612 0.0612 0.0706 0.0612 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ MONTE_15:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.122 0.129 0.132 0.127 0.137 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ ALFOR_16:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.194 0.196 0.198 0.198 0.201 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  $ PARIS_05:List of 8
##   ..$ FUN_CRIT: chr "ErrorCrit_KGE2"
##   ..$ Obs     : num [1:18628] 0.176 0.179 0.183 0.187 0.191 ...
##   ..$ VarObs  : chr "Q"
##   ..$ BoolCrit: logi [1:18628] TRUE TRUE TRUE TRUE TRUE TRUE ...
##   ..$ idLayer : logi NA
##   ..$ transfo : chr ""
##   ..$ epsilon : NULL
##   ..$ Weights : NULL
##   ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit"
##  - attr(*, "class")= chr [1:2] "list" "GRiwrmInputsCrit"

GRiwrmCalibOptions object

CalibOptions <- CreateCalibOptions(InputsModel)
str(CalibOptions)
## List of 25
##  $ TRANN_01:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ STDIZ_04:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ BAR-S_06:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ CHAUM_07:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ CUSSY_08:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ STGER_09:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ GUILL_10:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ AISY-_11:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ EPISY_14:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ MONTR_18:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ LOUVE_19:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ LASSI_20:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ VITRY_25:List of 4
##   ..$ FixedParam       : logi [1:4] NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:4] 4.59e-05 2.18e+04 -1.09e+04 1.09e+04 4.59e-05 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:4] 169.017 247.151 432.681 -2.376 -0.649 ...
##   ..- attr(*, "class")= chr [1:3] "CalibOptions" "GR4J" "HBAN"
##  $ MERY-_22:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ GURGY_02:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ CHABL_12:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ BRIEN_03:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ ARCIS_24:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ CHALO_21:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ COURL_23:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ NOISI_17:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ NOGEN_13:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ MONTE_15:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ ALFOR_16:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"
##  $ PARIS_05:List of 4
##   ..$ FixedParam       : logi [1:5] NA NA NA NA NA
##   ..$ SearchRanges     : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ...
##   ..$ FUN_TRANSFO      :function (ParamIn, Direction)  
##   ..$ StartParamDistrib: num [1:3, 1:5] 11.2 12.5 15 169 247.2 ...
##   ..- attr(*, "class")= chr [1:4] "CalibOptions" "GR4J" "SD" "HBAN"

Calibration

OutputsCalib <- Calibration(InputsModel, RunOptions, InputsCrit, CalibOptions)
## Calibration.GRiwrmInputsModel: Treating sub-basin TRANN_01...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  247.151,   -0.020,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8828
## Steepest-descent local search in progress
##   Calibration completed (47 iterations, 454 runs)
##       Param =  202.172,   -0.124,   76.792,    5.505
##       Crit. KGE2[Q]      = 0.9549
## Calibration.GRiwrmInputsModel: Treating sub-basin STDIZ_04...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  169.017,   -0.020,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.9015
## Steepest-descent local search in progress
##   Calibration completed (34 iterations, 342 runs)
##       Param =  161.253,   -0.238,   69.848,    3.746
##       Crit. KGE2[Q]      = 0.9472
## Calibration.GRiwrmInputsModel: Treating sub-basin BAR-S_06...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  247.151,   -0.020,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8826
## Steepest-descent local search in progress
##   Calibration completed (27 iterations, 281 runs)
##       Param =  247.151,   -0.347,   91.836,    5.302
##       Crit. KGE2[Q]      = 0.9588
## Calibration.GRiwrmInputsModel: Treating sub-basin CHAUM_07...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  432.681,   -0.020,   83.096,    1.417
##       Crit. KGE2[Q]      = 0.9019
## Steepest-descent local search in progress
##   Calibration completed (32 iterations, 331 runs)
##       Param =  269.272,    0.421,  166.665,    1.257
##       Crit. KGE2[Q]      = 0.9434
## Calibration.GRiwrmInputsModel: Treating sub-basin CUSSY_08...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  247.151,   -0.649,   42.098,    2.384
##       Crit. KGE2[Q]      = 0.8985
## Steepest-descent local search in progress
##   Calibration completed (26 iterations, 278 runs)
##       Param =  214.586,   -0.974,   59.761,    2.087
##       Crit. KGE2[Q]      = 0.9066
## Calibration.GRiwrmInputsModel: Treating sub-basin STGER_09...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  432.681,   -0.020,   83.096,    1.417
##       Crit. KGE2[Q]      = 0.9225
## Steepest-descent local search in progress
##   Calibration completed (28 iterations, 297 runs)
##       Param =  325.478,   -0.313,  142.419,    1.145
##       Crit. KGE2[Q]      = 0.9376
## Calibration.GRiwrmInputsModel: Treating sub-basin GUILL_10...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  169.017,   -2.376,   42.098,    2.384
##       Crit. KGE2[Q]      = 0.8562
## Steepest-descent local search in progress
##   Calibration completed (32 iterations, 331 runs)
##       Param =  183.842,   -1.733,   24.004,    2.586
##       Crit. KGE2[Q]      = 0.8879
## Calibration.GRiwrmInputsModel: Treating sub-basin AISY-_11...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  169.017,   -2.376,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8725
## Steepest-descent local search in progress
##   Calibration completed (41 iterations, 411 runs)
##       Param =  172.375,   -1.653,   49.462,    2.668
##       Crit. KGE2[Q]      = 0.9175
## Calibration.GRiwrmInputsModel: Treating sub-basin EPISY_14...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  432.681,   -0.649,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8520
## Steepest-descent local search in progress
##   Calibration completed (36 iterations, 369 runs)
##       Param =  610.285,   -0.593,   41.294,    3.754
##       Crit. KGE2[Q]      = 0.9395
## Calibration.GRiwrmInputsModel: Treating sub-basin MONTR_18...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  247.151,   -0.649,   42.098,    2.384
##       Crit. KGE2[Q]      = 0.8528
## Steepest-descent local search in progress
##   Calibration completed (18 iterations, 208 runs)
##       Param =  292.949,   -0.578,   40.854,    2.169
##       Crit. KGE2[Q]      = 0.8634
## Calibration.GRiwrmInputsModel: Treating sub-basin LOUVE_19...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  169.017,   -0.649,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8814
## Steepest-descent local search in progress
##   Calibration completed (17 iterations, 200 runs)
##       Param =  162.390,   -1.099,   83.096,    3.438
##       Crit. KGE2[Q]      = 0.9138
## Calibration.GRiwrmInputsModel: Treating sub-basin LASSI_20...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  247.151,   -0.649,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8510
## Steepest-descent local search in progress
##   Calibration completed (32 iterations, 322 runs)
##       Param =  223.632,   -1.160,   75.944,    4.355
##       Crit. KGE2[Q]      = 0.9075
## Calibration.GRiwrmInputsModel: Treating sub-basin VITRY_25...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (81 runs)
##       Param =  432.681,   -0.649,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8652
## Steepest-descent local search in progress
##   Calibration completed (83 iterations, 785 runs)
##       Param =  278.332,   -1.382,   97.715,    5.076
##       Crit. KGE2[Q]      = 0.9550
## Calibration.GRiwrmInputsModel: Treating sub-basin MERY-_22...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   42.098,    1.417
##       Crit. KGE2[Q]      = 0.8785
## Steepest-descent local search in progress
##   Calibration completed (77 iterations, 1000 runs)
##       Param =    0.413,  620.534,   -2.454,   27.938,   16.854
##       Crit. KGE2[Q]      = 0.9606
## Calibration.GRiwrmInputsModel: Treating sub-basin GURGY_02...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.9454
## Steepest-descent local search in progress
##   Calibration completed (85 iterations, 1111 runs)
##       Param =    2.132,  412.218,   -2.931,   97.170,    3.419
##       Crit. KGE2[Q]      = 0.9623
## Calibration.GRiwrmInputsModel: Treating sub-basin CHABL_12...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  247.151,   -2.376,   42.098,    2.384
##       Crit. KGE2[Q]      = 0.8633
## Steepest-descent local search in progress
##   Calibration completed (60 iterations, 826 runs)
##       Param =    0.830,  198.343,   -2.761,   46.063,    5.526
##       Crit. KGE2[Q]      = 0.9386
## Calibration.GRiwrmInputsModel: Treating sub-basin BRIEN_03...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  169.017,   -0.649,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.8976
## Steepest-descent local search in progress
##   Calibration completed (71 iterations, 946 runs)
##       Param =    1.381,  204.434,   -0.277,   49.597,    5.416
##       Crit. KGE2[Q]      = 0.9402
## Calibration.GRiwrmInputsModel: Treating sub-basin ARCIS_24...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   83.096,    2.384
##       Crit. KGE2[Q]      = 0.9094
## Steepest-descent local search in progress
##   Calibration completed (86 iterations, 1091 runs)
##       Param =    0.435,  431.915,   -3.245,   61.559,   14.581
##       Crit. KGE2[Q]      = 0.9632
## Calibration.GRiwrmInputsModel: Treating sub-basin CHALO_21...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   83.096,    1.417
##       Crit. KGE2[Q]      = 0.8970
## Steepest-descent local search in progress
##   Calibration completed (60 iterations, 830 runs)
##       Param =    0.433,  527.943,   -3.949,  151.618,   11.002
##       Crit. KGE2[Q]      = 0.9666
## Calibration.GRiwrmInputsModel: Treating sub-basin COURL_23...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   83.096,    1.417
##       Crit. KGE2[Q]      = 0.9154
## Steepest-descent local search in progress
##   Calibration completed (109 iterations, 1357 runs)
##       Param =    1.235,  633.742,   -3.667,  244.182,    2.730
##       Crit. KGE2[Q]      = 0.9657
## Calibration.GRiwrmInputsModel: Treating sub-basin NOISI_17...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   42.098,    1.417
##       Crit. KGE2[Q]      = 0.8131
## Steepest-descent local search in progress
##   Calibration completed (87 iterations, 1103 runs)
##       Param =    0.622, 1393.644,   -4.706,   54.015,    5.612
##       Crit. KGE2[Q]      = 0.9602
## Calibration.GRiwrmInputsModel: Treating sub-basin NOGEN_13...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   20.697,    1.417
##       Crit. KGE2[Q]      = 0.9172
## Steepest-descent local search in progress
##   Calibration completed (164 iterations, 1903 runs)
##       Param =    0.235, 1277.068,   -1.667,   13.408,    7.013
##       Crit. KGE2[Q]      = 0.9666
## Calibration.GRiwrmInputsModel: Treating sub-basin MONTE_15...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   20.697,    1.417
##       Crit. KGE2[Q]      = 0.9555
## Steepest-descent local search in progress
##   Calibration completed (86 iterations, 1087 runs)
##       Param =    0.340, 3983.834,   -4.643,    9.025,    1.574
##       Crit. KGE2[Q]      = 0.9766
## Calibration.GRiwrmInputsModel: Treating sub-basin ALFOR_16...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   20.697,    1.417
##       Crit. KGE2[Q]      = 0.9569
## Steepest-descent local search in progress
##   Calibration completed (37 iterations, 593 runs)
##       Param =    1.670, 1085.721,   -2.376,   21.542,    2.540
##       Crit. KGE2[Q]      = 0.9784
## Calibration.GRiwrmInputsModel: Treating sub-basin PARIS_05...
## Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
##   Screening completed (243 runs)
##       Param =   11.250,  432.681,   -2.376,   20.697,    2.384
##       Crit. KGE2[Q]      = 0.9709
## Steepest-descent local search in progress
##   Calibration completed (64 iterations, 869 runs)
##       Param =    1.337, 3070.182, -2199.053,   21.542,    3.644
##       Crit. KGE2[Q]      = 0.9748

Run model with Michel calibration

ParamMichel <- sapply(griwrm$id, function(x) {OutputsCalib[[x]]$Param})

OutputsModels <- RunModel(
  InputsModel,
  RunOptions = RunOptions,
  Param = ParamMichel
)
## RunModel.GRiwrmInputsModel: Treating sub-basin TRANN_01...
## RunModel.GRiwrmInputsModel: Treating sub-basin STDIZ_04...
## RunModel.GRiwrmInputsModel: Treating sub-basin BAR-S_06...
## RunModel.GRiwrmInputsModel: Treating sub-basin CHAUM_07...
## RunModel.GRiwrmInputsModel: Treating sub-basin CUSSY_08...
## RunModel.GRiwrmInputsModel: Treating sub-basin STGER_09...
## RunModel.GRiwrmInputsModel: Treating sub-basin GUILL_10...
## RunModel.GRiwrmInputsModel: Treating sub-basin AISY-_11...
## RunModel.GRiwrmInputsModel: Treating sub-basin EPISY_14...
## RunModel.GRiwrmInputsModel: Treating sub-basin MONTR_18...
## RunModel.GRiwrmInputsModel: Treating sub-basin LOUVE_19...
## RunModel.GRiwrmInputsModel: Treating sub-basin LASSI_20...
## RunModel.GRiwrmInputsModel: Treating sub-basin VITRY_25...
## RunModel.GRiwrmInputsModel: Treating sub-basin MERY-_22...
## RunModel.GRiwrmInputsModel: Treating sub-basin GURGY_02...
## RunModel.GRiwrmInputsModel: Treating sub-basin CHABL_12...
## RunModel.GRiwrmInputsModel: Treating sub-basin BRIEN_03...
## RunModel.GRiwrmInputsModel: Treating sub-basin ARCIS_24...
## RunModel.GRiwrmInputsModel: Treating sub-basin CHALO_21...
## RunModel.GRiwrmInputsModel: Treating sub-basin COURL_23...
## RunModel.GRiwrmInputsModel: Treating sub-basin NOISI_17...
## RunModel.GRiwrmInputsModel: Treating sub-basin NOGEN_13...
## RunModel.GRiwrmInputsModel: Treating sub-basin MONTE_15...
## RunModel.GRiwrmInputsModel: Treating sub-basin ALFOR_16...
## RunModel.GRiwrmInputsModel: Treating sub-basin PARIS_05...

Save calibration data for next vignettes

save(ParamMichel, file = "_cache/V03.RData")

Plot the result for each basin

plot(OutputsModels, Qobs = Qnat[IndPeriod_Run,])

## [[1]]
## NULL
## 
## [[2]]
## NULL
## 
## [[3]]
## NULL
## 
## [[4]]
## NULL
## 
## [[5]]
## NULL
## 
## [[6]]
## NULL
## 
## [[7]]
## NULL
## 
## [[8]]
## NULL
## 
## [[9]]
## NULL
## 
## [[10]]
## NULL
## 
## [[11]]
## NULL
## 
## [[12]]
## NULL
## 
## [[13]]
## NULL
## 
## [[14]]
## NULL
## 
## [[15]]
## NULL
## 
## [[16]]
## NULL
## 
## [[17]]
## NULL
## 
## [[18]]
## NULL
## 
## [[19]]
## NULL
## 
## [[20]]
## NULL
## 
## [[21]]
## NULL
## 
## [[22]]
## NULL
## 
## [[23]]
## NULL
## 
## [[24]]
## NULL
## 
## [[25]]
## NULL