Remove warnings about parameter threshold during GR5H calibration
Try this:
library(airGR)
example("RunModel_GR5H")
str(RunOptions) # Line ignored when all lines are copy-paste together. Don't know why
CalibOptions <- CreateCalibOptions(RunModel_GR5H)
OutputsCalib <- Calibration(InputsModel = InputsModel, RunOptions = RunOptions,
InputsCrit = InputsCrit, CalibOptions = CalibOptions,
FUN_MOD = RunModel_GR5H)
You get:
Grid-Screening in progress (0% 20% 40% 60% 80% 100%)
Screening completed (243 runs)
Param = 55.147, -0.400, 228.149, 11.050, 0.126
Crit. NSE[Q] = 0.0650
Steepest-descent local search in progress
Calibration completed (68 iterations, 913 runs)
Param = 451.901, -0.494, 143.170, 3.555, 0.188
Crit. NSE[Q] = 0.9184
There were 24 warnings (use warnings() to see them)
And then:
> warnings()
Warning messages:
1: In FUN_MOD(InputsModel = InputsModel, RunOptions = RunOptions, ... :
Param[4] (X4: unit hydrograph time constant [h]) < 0.50
X4 set to 0.50
2: In FUN_MOD(InputsModel = InputsModel, RunOptions = RunOptions, ... :
Param[4] (X4: unit hydrograph time constant [h]) < 0.50
X4 set to 0.50
...
I suppose that the bounds of the parameter should be correctly handled in the calibration to avoid this kind of warnings.