Skip to content

GitLab

  • Menu
Projects Groups Snippets
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • airGR airGR
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 60
    • Issues 60
    • List
    • Boards
    • Service Desk
    • Milestones
  • Redmine
    • Redmine
  • Merge requests 7
    • Merge requests 7
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • HYCAR-Hydro
  • airGRairGR
  • Issues
  • #125

Closed
Open
Created May 28, 2021 by Dorchies David@david.dorchiesDeveloper

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.

Edited Sep 14, 2021 by Delaigue Olivier
Assignee
Assign to
Time tracking