RunModel_GR4H.Rd 4.9 KB
 Delaigue Olivier committed Oct 28, 2015 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 % Generated by roxygen2 (4.1.1): do not edit by hand % Please edit documentation in R/RunModel_GR4H.R \encoding{UTF-8} \name{RunModel_GR4H} \alias{RunModel_GR4H} \title{Run with the GR4H hydrological model} \usage{ RunModel_GR4H(InputsModel, RunOptions, Param) } \arguments{ \item{InputsModel}{[object of class \emph{InputsModel}] see \code{\link{CreateInputsModel}} for details} \item{RunOptions}{[object of class \emph{RunOptions}] see \code{\link{CreateRunOptions}} for details} \item{Param}{[numeric] vector of 4 parameters \tabular{ll}{ GR4H X1 \tab production store capacity [mm] \cr GR4H X2 \tab groundwater exchange coefficient [mm/h] \cr GR4H X3 \tab routing store capacity [mm] \cr GR4H X4 \tab unit hydrograph time constant [h] \cr }} } \value{ [list] list containing the function outputs organised as follows: \tabular{ll}{ \emph{$DatesR } \tab [POSIXlt] series of dates \cr \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/h] \cr \emph{$Precip } \tab [numeric] series of input total precipitation [mm/h] \cr \emph{$Prod } \tab [numeric] series of production store level (X(2)) [mm] \cr \emph{$AE } \tab [numeric] series of actual evapotranspiration [mm/h] \cr \emph{$Perc } \tab [numeric] series of percolation (PERC) [mm/h] \cr \emph{$PR } \tab [numeric] series of PR=PN-PS+PERC [mm/h] \cr \emph{$Q9 } \tab [numeric] series of HU1 outflow (Q9) [mm/h] \cr \emph{$Q1 } \tab [numeric] series of HU2 outflow (Q1) [mm/h] \cr \emph{$Rout } \tab [numeric] series of routing store level (X(1)) [mm] \cr \emph{$Exch } \tab [numeric] series of potential semi-exchange between catchments [mm/h] \cr \emph{$AExch } \tab [numeric] series of actual exchange between catchments (1+2) [mm/h] \cr \emph{$QR } \tab [numeric] series of routing store outflow (QR) [mm/h] \cr \emph{$QD } \tab [numeric] series of direct flow from HU2 after exchange (QD) [mm/h] \cr \emph{$Qsim } \tab [numeric] series of Qsim [mm/h] \cr \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, HU1 levels, HU2 levels) [mm] \cr } (refer to the provided references or to the package source code for further details on these model outputs) } \description{ Function which performs a single run for the GR4H hourly lumped model. } \details{ For further details on the model, see the references section. For further details on the argument structures and initialisation options, see \code{\link{CreateRunOptions}}. } \examples{ ## load of catchment data require(airGR) 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="\%d/\%m/\%Y \%H:\%M")=="01/03/2004 00:00"), which(format(BasinObs$DatesR,format="\%d/\%m/\%Y \%H:\%M")=="31/12/2008 00:00")) ## preparation of the RunOptions object 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) ## results preview plot_OutputsModel(OutputsModel=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) cat(paste(" Crit ",OutputsCrit$CritName," ",round(OutputsCrit\$CritValue,4),"\\n",sep="")) } \author{ Laurent Coron (July 2014) } \references{ 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. } \seealso{ \code{\link{RunModel_GR4J}}, \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}. }