RunModel_GR4H {airGR} | R Documentation |
Function which performs a single run for the GR4H hourly lumped model.
RunModel_GR4H(InputsModel, RunOptions, Param)
InputsModel |
[object of class InputsModel] see | |||||||||
RunOptions |
[object of class RunOptions] see | |||||||||
Param |
[numeric] vector of 4 parameters
|
For further details on the model, see the references section.
For further details on the argument structures and initialisation options, see CreateRunOptions
.
[list] list containing the function outputs organised as follows:
$DatesR | [POSIXlt] series of dates |
$PotEvap | [numeric] series of input potential evapotranspiration [mm/h] |
$Precip | [numeric] series of input total precipitation [mm/h] |
$Prod | [numeric] series of production store level [mm] |
$AE | [numeric] series of actual evapotranspiration [mm/h] |
$Perc | [numeric] series of percolation (PERC) [mm/h] |
$PR | [numeric] series of PR=Pn-Ps+Perc [mm/h] |
$Q9 | [numeric] series of UH1 outflow (Q9) [mm/h] |
$Q1 | [numeric] series of UH2 outflow (Q1) [mm/h] |
$Rout | [numeric] series of routing store level [mm] |
$Exch | [numeric] series of potential semi-exchange between catchments [mm/h] |
$AExch | [numeric] series of actual exchange between catchments (1+2) [mm/h] |
$QR | [numeric] series of routing store outflow (QR) [mm/h] |
$QD | [numeric] series of direct flow from UH2 after exchange (QD) [mm/h] |
$Qsim | [numeric] series of simulated discharge [mm/h] |
$StateEnd | [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm], see CreateIniStates for more details |
(refer to the provided references or to the package source code for further details on these model outputs)
Laurent Coron, Charles Perrin, Thibaut Mathevet
Mathevet, T. (2005).
Quels modèles pluie-débit globaux pour le pas de temps horaire ?
Développement empirique et comparaison de modèles sur un large échantillon de bassins versants.
PhD thesis (in French), ENGREF - Cemagref (Antony), Paris, France.
Le Moine, N. (2008).
Le bassin versant de surface vu par le souterrain :
une voie d'amélioration des performances et du réalisme des modèles pluie-débit ?
PhD thesis (french), UPMC, Paris, France.
RunModel_GR4J
,
CreateInputsModel
, CreateRunOptions
, CreateIniStates
.
library(airGR) ## load of catchment data data(L0123003) ## preparation of the InputsModel object 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 = "%Y-%m-%d %H:%M")=="2004-03-01 00:00"), which(format(BasinObs$DatesR, format = "%Y-%m-%d %H:%M")=="2008-12-31 23:00")) ## preparation of the RunOptions object RunOptions <- CreateRunOptions(FUN_MOD = RunModel_GR4H, InputsModel = InputsModel, IndPeriod_Run = Ind_Run) ## simulation Param <- c(X1 = 521.113, X2 = -2.918, X3 = 218.009, X4 = 4.124) OutputsModel <- RunModel_GR4H(InputsModel = InputsModel, RunOptions = RunOptions, Param = Param) ## results preview plot(OutputsModel, Qobs = BasinObs$Qmm[Ind_Run]) ## efficiency criterion: Nash-Sutcliffe Efficiency InputsCrit <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsModel, RunOptions = RunOptions, Obs = BasinObs$Qmm[Ind_Run]) OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)