Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
HYCAR-Hydro
airGR
Commits
4458848b
Commit
4458848b
authored
May 22, 2018
by
Delaigue Olivier
Browse files
v1.0.10.5 typo revision of the example of the RunModel_GR4H doc
parent
c3f47687
Changes
3
Hide whitespace changes
Inline
Side-by-side
DESCRIPTION
View file @
4458848b
Package: airGR
Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
Version: 1.0.10.
4
Version: 1.0.10.
5
Date: 2018-05-22
Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl")),
...
...
NEWS.rmd
View file @
4458848b
...
...
@@ -14,7 +14,7 @@ output:
### 1.0.10.
4
Release Notes (2018-05-22)
### 1.0.10.
5
Release Notes (2018-05-22)
____________________________________________________________________________________
...
...
man/RunModel_GR4H.Rd
View file @
4458848b
...
...
@@ -70,27 +70,27 @@ library(airGR)
data(L0123003)
## preparation of the InputsModel object
InputsModel <- CreateInputsModel(FUN_MOD
=
RunModel_GR4H,DatesR
=
BasinObs$DatesR,
Precip
=
BasinObs$P,PotEvap
=
BasinObs$E)
InputsModel <- CreateInputsModel(FUN_MOD
=
RunModel_GR4H,
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/03/2004 00:00"),
which(format(BasinObs$DatesR,format
=
"\%d/\%m/\%Y \%H:\%M")=="31/12/2008 23:00"))
Ind_Run <- seq(which(format(BasinObs$DatesR,
format
=
"\%d/\%m/\%Y \%H:\%M")=="01/03/2004 00:00"),
which(format(BasinObs$DatesR,
format
=
"\%d/\%m/\%Y \%H:\%M")=="31/12/2008 23:00"))
## preparation of the RunOptions object
RunOptions <- CreateRunOptions(FUN_MOD
=
RunModel_GR4H,InputsModel
=
InputsModel,IndPeriod_Run
=
Ind_Run)
RunOptions <- CreateRunOptions(FUN_MOD
=
RunModel_GR4H,
InputsModel
=
InputsModel,
IndPeriod_Run
=
Ind_Run)
## simulation
Param <- c(521.113,-2.918,218.009,4.124)
OutputsModel <- RunModel_GR4H(InputsModel
=
InputsModel,RunOptions
=
RunOptions,Param
=
Param)
Param <- c(521.113,
-2.918,
218.009,
4.124)
OutputsModel <- RunModel_GR4H(InputsModel
=
InputsModel,
RunOptions
=
RunOptions,
Param
=
Param)
## results preview
plot(OutputsModel,Qobs
=
BasinObs$Qmm[Ind_Run])
plot(OutputsModel,
Qobs
=
BasinObs$Qmm[Ind_Run])
## efficiency criterion: Nash-Sutcliffe Efficiency
InputsCrit <- CreateInputsCrit(FUN_CRIT
=
ErrorCrit_NSE,InputsModel
=
InputsModel,
RunOptions
=
RunOptions,Qobs
=
BasinObs$Qmm[Ind_Run])
OutputsCrit <- ErrorCrit_NSE(InputsCrit
=
InputsCrit,OutputsModel
=
OutputsModel)
InputsCrit <- CreateInputsCrit(FUN_CRIT
=
ErrorCrit_NSE,
InputsModel
=
InputsModel,
RunOptions
=
RunOptions,
Qobs
=
BasinObs$Qmm[Ind_Run])
OutputsCrit <- ErrorCrit_NSE(InputsCrit
=
InputsCrit,
OutputsModel
=
OutputsModel)
}
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment