man/TransfoParam_CemaNeige.Rd

-\encoding{UTF-8}
-\name{TransfoParam_CemaNeige}
-\alias{TransfoParam_CemaNeige}
-\title{Transformation of the parameters of the CemaNeige module}
-\usage{
-TransfoParam_CemaNeige(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the CemaNeige module (from raw to transformed parameters and vice versa).
-}
-\examples{
-library(airGR)
-
-## transformation Raw->Transformed for the CemaNeige module
-    Xraw  <- matrix(c(+0.19, +1.73,
-                      +0.39, +2.51,
-                      +0.74, +4.06),
-                      ncol = 2, byrow = TRUE)
-    Xtran <- TransfoParam_CemaNeige(ParamIn = Xraw , Direction = "RT")
-
-## transformation Transformed->Raw for the CemaNeige module
-    Xtran <- matrix(c(-6.26, +0.55,
-                      -2.13, +0.92,
-                      +4.86, +1.40),
-                      ncol = 2, byrow = TRUE)
-    Xraw  <- TransfoParam_CemaNeige(ParamIn = Xtran, Direction = "TR")
-}
-\author{
-Laurent Coron (December 2013)
-}
-\seealso{
-\code{\link{TransfoParam}}, \code{\link{TransfoParam_GR4J}}, \code{\link{TransfoParam_GR5J}}, \code{\link{TransfoParam_GR6J}}
-}
-

man/TransfoParam_GR1A.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR1A}
-\alias{TransfoParam_GR1A}
-\title{Transformation of the parameters of the GR1A model}
-\usage{
-TransfoParam_GR1A(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR1A model (from real to transformed parameters and vice versa).
-}
-\author{
-Laurent Coron (March 2015)
-}
-

man/TransfoParam_GR2M.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR2M}
-\alias{TransfoParam_GR2M}
-\title{Transformation of the parameters of the GR2M model}
-\usage{
-TransfoParam_GR2M(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR2M model (from real to transformed parameters and vice versa).
-}
-\author{
-Laurent Coron (March 2015)
-}
-

man/TransfoParam_GR4H.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR4H}
-\alias{TransfoParam_GR4H}
-\title{Transformation of the parameters of the GR4H model}
-\usage{
-TransfoParam_GR4H(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR4H model (from real to transformed parameters and vice versa).
-}
-\author{
-Laurent Coron, Claude Michel, Thibault Mathevet (July 2014)
-}
-

man/TransfoParam_GR4J.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR4J}
-\alias{TransfoParam_GR4J}
-\title{Transformation of the parameters of the GR4J model}
-\usage{
-TransfoParam_GR4J(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR4J model (from real to transformed parameters and vice versa).
-}
-\examples{
-library(airGR)
-
-## transformation Raw->Transformed for the GR4J model
-    Xraw  <- matrix(c(+221.41, -3.63,  +30.00, +1.37,
-                      +347.23, -1.03,  +60.34, +1.76,
-                      +854.06, -0.10, +148.41, +2.34),
-                      ncol = 4, byrow = TRUE)
-    Xtran <- TransfoParam_GR4J(ParamIn = Xraw , Direction = "RT")
-
-## transformation Transformed->Raw for the GR4J model
-    Xtran <- matrix(c(+3.60, -2.00, +3.40, -9.10,
-                      +3.90, -0.90, +4.10, -8.70,
-                      +4.50, -0.10, +5.00, -8.10),
-                      ncol = 4, byrow = TRUE)
-    Xraw  <- TransfoParam_GR4J(ParamIn = Xtran, Direction = "TR")
-}
-\author{
-Laurent Coron (December 2013)
-}
-\seealso{
-\code{\link{TransfoParam}}, \code{\link{TransfoParam_GR5J}}, \code{\link{TransfoParam_GR6J}}, \code{\link{TransfoParam_CemaNeige}}
-}
-

man/TransfoParam_GR5J.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR5J}
-\alias{TransfoParam_GR5J}
-\title{Transformation of the parameters of the GR5J model}
-\usage{
-TransfoParam_GR5J(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR5J model (from real to transformed parameters and vice versa).
-}
-\examples{
-library(airGR)
-
-## transformation Raw->Transformed for the GR5J model
-    Xraw  <- matrix(c(+221.41, -2.65,  +27.11, +1.37, -0.76,
-                      +347.23, -0.64,  +60.34, +1.76, +0.30,
-                      +854.01, -0.10, +148.41, +2.34, +0.52),
-                      ncol = 5, byrow = TRUE)
-    Xtran <- TransfoParam_GR5J(ParamIn = Xraw , Direction = "RT")
-
-## transformation Transformed->Raw for the GR5J model
-    Xtran <- matrix(c(+3.60, -1.70, +3.30, -9.10, -0.70,
-                      +3.90, -0.60, +4.10, -8.70, +0.30,
-                      +4.50, -0.10, +5.00, -8.10, +0.50),
-                      ncol = 5, byrow = TRUE)
-    Xraw  <- TransfoParam_GR5J(ParamIn = Xtran, Direction = "TR")
-}
-\author{
-Laurent Coron (December 2013)
-}
-\seealso{
-\code{\link{TransfoParam}}, \code{\link{TransfoParam_GR4J}}, \code{\link{TransfoParam_GR6J}}, \code{\link{TransfoParam_CemaNeige}}
-}
-

man/TransfoParam_GR6J.Rd

-\encoding{UTF-8}
-\name{TransfoParam_GR6J}
-\alias{TransfoParam_GR6J}
-\title{Transformation of the parameters of the GR6J model}
-\usage{
-TransfoParam_GR6J(ParamIn, Direction)
-}
-\arguments{
-\item{ParamIn}{[numeric] matrix of parameter sets (sets in line, parameter values in column)}
-
-\item{Direction}{[character] direction of the transformation: use "RT" for Raw->Transformed and "TR" for Transformed->Raw}
-}
-\value{
-\emph{ParamOut} [numeric] matrix of parameter sets (sets in line, parameter values in column)
-}
-\description{
-Function which transforms model parameters of the GR6J model (from real to transformed parameters and vice versa).
-}
-\examples{
-library(airGR)
-
-## transformation Raw->Transformed for the GR6J model
-    Xraw  <- matrix(c(+221.41, -1.18,  +27.11, 1.37, -0.18,  +20.09,
-                      +347.23, -0.52,  +60.34, 1.76, +0.02,  +54.60,
-                      +854.06, +0.52, +148.41, 2.34, +0.22, +148.41),
-                      ncol = 6, byrow = TRUE)
-    Xtran <- TransfoParam_GR6J(ParamIn = Xraw , Direction = "RT")
-
-## transformation Transformed->Raw for the GR6J model
-    Xtran <- matrix(c(+3.60, -1.00, +3.30, -9.10, -0.90, +3.00,
-                      +3.90, -0.50, +4.10, -8.70, +0.10, +4.00,
-                      +4.50, +0.50, +5.00, -8.10, +1.10, +5.00),
-                      ncol = 6, byrow = TRUE)
-    Xraw  <- TransfoParam_GR6J(ParamIn = Xtran, Direction = "TR")
-}
-\author{
-Laurent Coron (December 2013)
-}
-\seealso{
-\code{\link{TransfoParam}}, \code{\link{TransfoParam_GR4J}}, \code{\link{TransfoParam_GR5J}}, \code{\link{TransfoParam_CemaNeige}}
-}
-

man/airGR-package.Rd

+\name{airGR-package}
+\alias{airGR}
+\alias{airGR-package}
+\docType{package}
+\encoding{UTF-8}
+\title{\packageTitle{airGR}}
+
+
+\description{
+\packageDescription{airGR}
+Each model core is coded in Fortran to ensure low computational time. The other package functions (i.e. mainly the calibration algorithm and the computation of the efficiency criteria) are coded in R. \cr\cr
+
+
+## --- Functions and objects
+
+The airGR package has been designed to fulfil two major requirements: facilitate the use by non-expert users and allow flexibility regarding the addition of external criteria, models or calibration algorithms. The names of the functions and their arguments were chosen to this end.
+
+The package is mostly based on three families of functions:
+\itemize{
+  \item the functions belonging to the \code{\link{RunModel}} family require three arguments: \emph{InputsModel}, \emph{RunOptions} and \emph{Param};  please refer to help pages \code{\link{CreateInputsModel}} and \code{\link{CreateRunOptions}} for further details and examples;
+  \item the functions belonging to the \code{\link{ErrorCrit}} family require two arguments: \emph{InputsCrit} and \emph{OutputsModel}; please refer to help pages \code{\link{CreateInputsCrit}} and \code{\link{RunModel}} for further details and examples;
+  \item the functions belonging to the \code{\link{Calibration}} family require four arguments: \emph{InputsModel}, \emph{RunOptions}, \emph{InputsCrit} and \emph{CalibOptions};  please refer to help pages \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateInputsCrit}} and \code{\link{CreateCalibOptions}} for further details and examples.
+}
+
+In order to limit the risk of mis-use and increase the flexibility of these main functions, we imposed the structure of their arguments and defined their class. Most users will not need to worry about these imposed structures since functions are provided to prepare these arguments for them: \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateInputsCrit}}, \code{\link{CreateCalibOptions}}. However, advanced users wishing to supplement the package with their own models will need to comply with these imposed structures and refer to the package source codes to get all the specification requirements. \cr\cr
+
+
+## --- Models
+
+Seven hydrological models and one snow melt and accumulation model are implemented in airGR. The hydrological models can be applied either on a lumped way or on a semi-distributed way (on sub-catchments). The snow model can either be used alone or with the daily or hourly hydrological models. Naturally each hydrological model can also be used alone. \cr
+These models can be called within airGR using the following functions:
+\itemize{
+  \item \code{\link{RunModel_GR4H}}: four-parameter hourly lumped hydrological model (Mathevet, 2005)
+  \item \code{\link{RunModel_GR5H}}: five-parameter hourly lumped hydrological model (Ficchi, 2017; Ficchi \emph{et al.}, 2019)
+  \item \code{\link{RunModel_GR4J}}: four-parameter daily lumped hydrological model (Perrin \emph{et al.}, 2003)
+  \item \code{\link{RunModel_GR5J}}: five-parameter daily lumped hydrological model (Le Moine, 2008)
+  \item \code{\link{RunModel_GR6J}}: six-parameter daily lumped hydrological model (Pushpalatha \emph{et al.}, 2011)
+  \item \code{\link{RunModel_GR2M}}: two-parameter monthly lumped hydrological model (Mouelhi, 2003; Mouelhi \emph{et al.}, 2006a)
+  \item \code{\link{RunModel_GR1A}}: one-parameter yearly lumped hydrological model (Mouelhi, 2003; Mouelhi \emph{et al.}, 2006b)
+  \item \code{\link{RunModel_CemaNeige}}: two-parameter degree-day snow melt and accumulation daily model (Valéry \emph{et al.}, 2014; Riboust \emph{et al.}, 2019)
+  \item \code{\link{RunModel_CemaNeigeGR4H}}: combined use of GR4H and CemaNeige
+  \item \code{\link{RunModel_CemaNeigeGR5H}}: combined use of GR5H and CemaNeige
+  \item \code{\link{RunModel_CemaNeigeGR4J}}: combined use of GR4J and CemaNeige
+  \item \code{\link{RunModel_CemaNeigeGR5J}}: combined use of GR5J and CemaNeige
+  \item \code{\link{RunModel_CemaNeigeGR6J}}: combined use of GR6J and CemaNeige
+}
+
+## --- How to get started
+
+To learn how to use the functions from the airGR package, it is recommended to follow the five steps described below:
+\enumerate{
+  \item refer to the help for \code{\link{RunModel_GR4J}} then run the provided example to assess how to make a simulation;
+  \item refer to the help for \code{\link{CreateInputsModel}} to understand how the inputs of a model are prepared/organised;
+  \item refer to the help for \code{\link{CreateRunOptions}} to understand how the run options of a model are parametrised/organised;
+  \item refer to the help for \code{\link{ErrorCrit_NSE}} and \code{\link{CreateInputsCrit}} to understand how the computation of an error criterion is prepared/made;
+  \item refer to the help for \code{\link{Calibration_Michel}}, run the provided example and then refer to the help for \code{\link{CreateCalibOptions}} to understand how a model calibration is prepared/made.
+}
+
+To get started with the package, you can refer to the 'get_started' vignette (\code{vignette("V01_get_started", package = "airGR")}). To know how to use the models on a semi-distributed way, you can refer to the 'sd_model' vignette (\code{vignette("V05_sd_model", package = "airGR")}). For more information, please visit the \href{https://hydrogr.github.io/airGR/}{airGR website}. \cr\cr
+
+
+## --- References
+
+\itemize{
+  \item Ficchi, A. (2017). An adaptive hydrological model for multiple time-steps: Diagnostics and improvements based on fluxes consistency. PhD thesis, UPMC - Irstea Antony, Paris, France.
+  \item Ficchi, A., Perrin, C. and Andréassian, V. (2019). Hydrological modelling at multiple sub-daily time steps: model improvement via flux-matching. Journal of Hydrology, 575, 1308-1327, \doi{10.1016/j.jhydrol.2019.05.084}.
+  \item 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 (in French), UPMC - Cemagref Antony, Paris, France, 324 pp.
+  \item 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, 463 pp.
+  \item Mouelhi, S. (2003). Vers une chaîne cohérente de modèles pluie-débit conceptuels globaux aux pas de temps pluriannuel, annuel, mensuel et journalier, PhD thesis (in French), ENGREF - Cemagref Antony, Paris, France, 323 pp.
+  \item Mouelhi, S., Michel, C., Perrin, C. and Andréassian, V. (2006a). Stepwise development of a two-parameter monthly water balance model. Journal of Hydrology, 318(1-4), 200-214, \doi{10.1016/j.jhydrol.2005.06.014}.
+  \item Mouelhi, S., Michel, C., Perrin, C. and Andréassian, V. (2006b). Linking stream flow to rainfall at the annual time step: the Manabe bucket model revisited. Journal of Hydrology, 328, 283-296, \doi{10.1016/j.jhydrol.2005.12.022}.
+  \item Perrin, C., Michel, C. and Andréassian, V. (2003). Improvement of a parsimonious model for streamflow simulation. Journal of Hydrology, 279(1-4), 275-289, \doi{10.1016/S0022-1694(03)00225-7}.
+  \item Pushpalatha, R., Perrin, C., Le Moine, N., Mathevet, T. and Andréassian, V. (2011). A downward structural sensitivity analysis of hydrological models to improve low-flow simulation. Journal of Hydrology, 411(1-2), 66-76, \doi{10.1016/j.jhydrol.2011.09.034}.
+  \item Riboust, P., Thirel, G., Le Moine N. and Ribstein P. (2019). Revisiting a simple degree-day model for integrating satellite data: Implementation of SWE-SCA hystereses. Journal of Hydrology and Hydromechanics, 67(1), 70–81, \doi{10.1016/j.jhydrol.2014.04.058}.
+  \item Valéry, A., Andréassian, V. and Perrin, C. (2014). "As simple as possible but not simpler": What is useful in a temperature-based snow-accounting routine? Part 2 - Sensitivity analysis of the Cemaneige snow accounting routine on 380 catchments. Journal of Hydrology, 517(0), 1176-1187, \doi{10.1016/j.jhydrol.2014.04.058}.
+}
+
+}
+
+\concept{hydrology}
+\concept{model}
+\concept{efficiency criterion}
+\concept{calibration}
+\concept{GR4J}

man/airGR.Rd

-\name{airGR}
-\alias{airGR}
-\docType{package}
-\encoding{UTF-8}
-\title{Suite of GR Hydrological Models for Precipitation-Runoff Modelling}
-\description{
-This package brings into R the hydrological modelling tools used at Irstea-Antony (HBAN Research Unit, France), including rainfall-runoff models (GR4H, GR4J, GR5J, GR6J, GR2M, GR1A) and a snowmelt module (CemaNeige). Each model core is coded in FORTRAN to ensure low computational time. The other package functions (i.e. mainly the calibration algorithm and the computation of the efficiency criteria) are coded in R. \cr
-
-
-##### Functions and objects #####
-
-The airGR package has been designed to fulfil two major requirements: facilitate the use by non-expert users and allow flexibility regarding the addition of external criteria, models or calibration algorithms. The names of the functions and their arguments were chosen to this end. 
-
-The package is mostly based on three families of functions: \cr
-- the functions belonging to the \code{\link{RunModel}} family require three arguments: \emph{InputsModel}, \emph{RunOptions} and \emph{Param};  please refer to help pages \code{\link{CreateInputsModel}} and \code{\link{CreateRunOptions}} for further details and examples; \cr
-- the functions belonging to the \code{\link{ErrorCrit}} family require two arguments: \emph{InputsCrit} and \emph{OutputsModel}; please refer to help pages \code{\link{CreateInputsCrit}} and \code{\link{RunModel}} for further details and examples; \cr
-- the functions belonging to the \code{\link{Calibration}} family require four arguments: \emph{InputsModel}, \emph{RunOptions}, \emph{InputsCrit} and \emph{CalibOptions};  please refer to help pages \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateInputsCrit}} and \code{\link{CreateCalibOptions}} for further details and examples.
-
-In order to limit the risk of mis-use and increase the flexibility of these main functions, we imposed the structure of their arguments and defined their class. Most users will not need to worry about these imposed structures since functions are provided to prepare these arguments for them: \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateInputsCrit}}, \code{\link{CreateCalibOptions}}. However, advanced users wishing to supplement the package with their own models will need to comply with these imposed structures and refer to the package source codes to get all the specification requirements. \cr
-
-
-##### Models #####
-
-Six hydrological models and one snow melt and accumulation module are implemented in airGR. The snow module can also be used alone or with the daily hydrological models, and each hydrological model can either be used alone or together with the snow module. \cr
-These models can be called within airGR using the following functions: \cr 
-- \code{\link{RunModel_GR4H}}: four-parameter hourly lumped hydrological model (Mathevet, 2005) \cr
-- \code{\link{RunModel_GR4J}}: four-parameter daily lumped hydrological model (Perrin et al., 2003) \cr
-- \code{\link{RunModel_GR5J}}: five-parameter daily lumped hydrological model (Le Moine, 2008) \cr
-- \code{\link{RunModel_GR6J}}: six-parameter daily lumped hydrological model (Pushpalatha et al., 2011) \cr
-- \code{\link{RunModel_GR2M}}: two-parameter monthly lumped hydrological model (Mouelhi, 2003 ; Mouelhi et al., 2006a) \cr
-- \code{\link{RunModel_GR1A}}: one-parameter yearly lumped hydrological model (Mouelhi, 2003 ; Mouelhi et al., 2006b) \cr
-- \code{\link{RunModel_CemaNeige}}: two-parameter degree-day snow melt and accumulation daily module (Valéry et al., 2014) \cr
-- \code{\link{RunModel_CemaNeigeGR4J}}: combined use of GR4J and CemaNeige \cr
-- \code{\link{RunModel_CemaNeigeGR5J}}: combined use of GR5J and CemaNeige \cr
-- \code{\link{RunModel_CemaNeigeGR6J}}: combined use of GR6J and CemaNeige \cr
-
-
-##### How to get started #####
-
-To learn how to use the functions from the airGR package, it is recommended to follow the five steps described below: \cr
-1. refer to the help for \code{\link{RunModel_GR4H}} then run the provided example to assess how to make a simulation ; \cr
-2. refer to the help for \code{\link{CreateInputsModel}} to understand how the inputs of a model are prepared/organised ; \cr
-3. refer to the help for \code{\link{CreateRunOptions}} to understand how the run options of a model are parametrised/organised ; \cr
-4. refer to the help for \code{\link{ErrorCrit_NSE}} and \code{\link{CreateInputsCrit}} to understand how the computation of an error criterion is prepared/made ; \cr
-5. refer to the help for \code{\link{Calibration_Michel}}, run the provided example and then refer to the help for \code{\link{CreateCalibOptions}} to understand how a model calibration is prepared/made ; \cr
-
-
-##### References #####
-
-- 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 (in French), UPMC - Cemagref Antony, Paris, France, 324 pp. \cr
-- 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, 463 pp. \cr
-- Mouelhi S. (2003), Vers une chaîne cohérente de modèles pluie-débit conceptuels globaux aux pas de temps pluriannuel, annuel, mensuel et journalier, PhD thesis (in French), ENGREF - Cemagref Antony, Paris, France, 323 pp. \cr
-- Mouelhi, S., C. Michel, C. Perrin and V. Andréassian (2006a), Stepwise development of a two-parameter monthly water balance model, Journal of Hydrology, 318(1-4), 200-214, doi:10.1016/j.jhydrol.2005.06.014. \cr
-- Mouelhi, S., C. Michel, C. Perrin. & V. Andreassian (2006b), Linking stream flow to rainfall at the annual time step: the Manabe bucket model revisited, Journal of Hydrology, 328, 283-296, doi:10.1016/j.jhydrol.2005.12.022. \cr
-- Perrin, C., C. Michel and V. Andréassian (2003), Improvement of a parsimonious model for streamflow simulation, Journal of Hydrology, 279(1-4), 275-289, doi:10.1016/S0022-1694(03)00225-7. \cr
-- Pushpalatha, R., C. Perrin, N. Le Moine, T. Mathevet and V. Andréassian (2011), A downward structural sensitivity analysis of hydrological models to improve low-flow simulation, Journal of Hydrology, 411(1-2), 66-76, doi:10.1016/j.jhydrol.2011.09.034. \cr
-- Valéry, A., V. Andréassian and C. Perrin (2014), "As simple as possible but not simpler": What is useful in a temperature-based snow-accounting routine? Part 2 - Sensitivity analysis of the Cemaneige snow accounting routine on 380 catchments, Journal of Hydrology, 517(0): 1176-1187, doi: 1176-1187, doi:10.1016/j.jhydrol.2014.04.058. \cr
-
-}
-
-\keyword{hydrology, model, efficiency criterion, calibration, GR4J}

man/plot.OutputsModel.Rd

 \encoding{UTF-8}
-\name{plot.OutputsModel}
+
+
+\name{plot}
 \alias{plot.OutputsModel}
-\alias{plot_OutputsModel}
+\alias{plot}
+\alias{exampleSimPlot}
+\alias{simGR4J}
+\alias{simCNGR4J}
+
+
 \title{Default preview of model outputs}
+
+
+\description{
+Function which creates a screen plot giving an overview of the model outputs.
+}
+
+
 \usage{
 \method{plot}{OutputsModel}(x, Qobs = NULL, IndPeriod_Plot = NULL,
-  BasinArea = NULL, which = "all", log_scale = FALSE,
-  cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1, verbose = TRUE, ...)
+     BasinArea = NULL, which = "synth", log_scale = FALSE,
+     cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1,
+     AxisTS = function(x) axis.POSIXct(side = 1, x = x$DatesR, ...),
+     LayoutMat = NULL,
+     LayoutWidths = rep.int(1, ncol(LayoutMat)),
+     LayoutHeights = rep.int(1, nrow(LayoutMat)),
+     verbose = TRUE, ...)
 }
+
+
 \arguments{
-\item{x}{[object of class \emph{OutputsModel}] list of model outputs (which must at least include DatesR, Precip and Qsim) [POSIXlt, mm, mm]}
+\item{x}{[object of class \emph{OutputsModel}] list of model outputs (which must at least include DatesR, Precip and Qsim) [POSIXlt, mm/time step, mm/time step]}
 
 \item{Qobs}{(optional) [numeric] time series of observed flow (for the same time steps than simulated) [mm/time step]}
 
@@ -17,9 +38,9 @@
 
 \item{BasinArea}{(optional) [numeric] basin area [km2], used to plot flow axes in m3/s}
 
-\item{which}{(optional) [character] choice of plots \cr (e.g. c(\code{"Precip"}, \code{"Temp"}, \code{"SnowPack"}, \code{"Flows"}, \code{"Regime"}, \code{"CumFreq"}, \code{"CorQQ"})), default = \code{"all"}}
+\item{which}{(optional) [character] choice of plots \cr (e.g. c(\code{"Precip"}, \code{"Temp"}, \code{"SnowPack"}, \code{"Flows"}, \code{"Error"}, \code{"Regime"}, \code{"CumFreq"}, \code{"CorQQ"})), default = \code{"synth"}, see details below}
 
-\item{log_scale}{(optional) [boolean] indicating if the flow axis is to be logarithmic, default = \code{FALSE}}
+\item{log_scale}{(optional) [boolean] indicating if the flow and the error time series axis and the flow error time series axis are to be logarithmic, default = \code{FALSE}}
 
 \item{cex.axis}{(optional) [numeric] the magnification to be used for axis annotation relative to the current setting of \code{cex}}
 
@@ -29,29 +50,111 @@
 
 \item{lwd}{(optional) [numeric] the line width (a positive number)}
 
+\item{AxisTS}{(optional) [function] to manage x-axis representing calendar dates and times on time series plots (see \code{\link{axis}} and \code{\link{axis.POSIXct}})}
+
+\item{LayoutMat}{(optional) [numeric] a matrix object specifying the location of the next N figures on the output device. Each value in the matrix must be 0 or a positive integer. If N is the largest positive integer in the matrix, then the integers \eqn{1, \dots, N-1} must also appear at least once in the matrix (see \code{\link{layout}})}
+
+\item{LayoutWidths}{(optional) [numeric] a vector of values for the widths of columns on the device (see \code{\link{layout}})}
+
+\item{LayoutHeights}{(optional) [numeric] a vector of values for the heights of rows on the device (see \code{\link{layout}})}
+
 \item{verbose}{(optional) [boolean] indicating if the function is run in verbose mode or not, default = \code{TRUE}}
 
-\item{...}{other parameters to be passed through to plotting functions}
+\item{...}{(optional) other parameters to be passed through to plotting functions}
 }
+
+
 \value{
-screen plot window
-}
-\description{
-Function which creates a screen plot giving an overview of the model outputs
+Screen plot window.
 }
+
+
 \details{
-Dashboard of results including various graphs (depending on the model):\cr
-  (1) time series of total precipitation\cr
-  (2) time series of temperature (plotted only if CemaNeige is used)\cr
-  (3) time series of snow pack (plotted only if CemaNeige is used)\cr
-  (4) time series of simulated flows (and observed flows if provided)\cr
-  (5) interannual median monthly simulated flow (and observed flows if provided)\cr
-  (6) correlation plot between simulated and observed flows (if observed flows provided)\cr
-  (7) cumulative frequency plot for simulated flows (and observed flows if provided)
+Different types of independent graphs are available (depending on the model, but always drawn in this order):
+\itemize{
+  \item \code{"Precip"}: time series of total precipitation
+  \item \code{"PotEvap"}: time series of potential evapotranspiration
+  \item \code{"ActuEvap"}: time series of simulated actual evapotranspiration (overlaid to \code{"PotEvap"} if already drawn)
+  \item \code{"Temp"}: time series of temperature (plotted only if CemaNeige is used)
+  \item \code{"SnowPack"}: time series of snow water equivalent (plotted only if CemaNeige is used)
+  \item \code{"Flows"}: time series of simulated flows (and observed flows if provided)
+  \item \code{"Error"}: time series of simulated flows minus observed flows (and observed flows if provided)
+  \item \code{"Regime"}: centred 30-day rolling mean applied on interannual average of daily simulated flows (and observed flows if provided)
+  \item \code{"CorQQ"}: correlation plot between simulated and observed flows (only if observed flows provided)
+  \item \code{"CumFreq"}: cumulative frequency plot for simulated flows (and observed flows if provided)
+}
+
+Different dashboards of results including various graphs are available:
+\itemize{
+  \item \code{"perf"}: corresponds to \code{"Error"}, \code{"Regime"}, \code{"CumFreq"} and \code{"CorQQ"}
+  \item \code{"ts"}: corresponds to \code{"Precip"}, \code{"PotEvap"}, \code{"Temp"}, \code{"SnowPack"} and \code{"Flows"}
+  \item \code{"synth"}: corresponds to \code{"Precip"}, \code{"Temp"}, \code{"SnowPack"}, \code{"Flows"}, \code{"Regime"}, \code{"CumFreq"} and \code{"CorQQ"}
+  \item \code{"all"}: corresponds to \code{"Precip"}, \code{"PotEvap"}, \code{"ActuEvap"}, \code{"Temp"}, \code{"SnowPack"}, \code{"Flows"}, \code{"Error"}, \code{"Regime"}, \code{"CumFreq"} and \code{"CorQQ"}
+  }
+
+If several dashboards are selected, or if an independent graph is called with a dashboard, the graphical device will include all the requested graphs without redundancy.
 }
+
 \author{
-Laurent Coron (June 2014),
-Olivier Delaigue,
-Guillaume Thirel
+Laurent Coron, Olivier Delaigue, Guillaume Thirel
+}
+
+
+\examples{
+### see examples of RunModel_GR4J or RunModel_CemaNeigeGR4J functions
+### to understand how the example datasets have been prepared
+
+
+## loading examples dataset for GR4J and GR4J + CemaNeige
+data(exampleSimPlot)
+
+
+### Qobs and outputs from GR4J and GR4J + CemaNeige models
+str(simGR4J, max.level = 1)
+str(simCNGR4J, max.level = 1)
+
+
+### default dashboard (which = "synth")
+
+## GR models whithout CemaNeige
+plot(simGR4J$OutputsModel, Qobs = simGR4J$Qobs)
+
+## GR models whith CemaNeige ("Temp" and "SnowPack" added)
+plot(simCNGR4J$OutputsModel, Qobs = simCNGR4J$Qobs)
+
+
+### "Error" and "CorQQ" plots cannot be display whithout Qobs
+plot(simGR4J$OutputsModel, which = "all", Qobs = simGR4J$Qobs)
+plot(simGR4J$OutputsModel, which = "all", Qobs = NULL)
+
+
+### complex plot arrangements
+plot(simGR4J$OutputsModel, Qobs = simGR4J$Qobs,
+     which = c("Flows", "Regime", "CumFreq", "CorQQ"),
+     LayoutMat = matrix(c(1, 2, 3, 1, 4, 4), ncol = 2),
+     LayoutWidths  = c(1.5, 1),
+     LayoutHeights = c(0.5, 1, 1))
+
+
+### customize x-axis on time series plot
+FunAxisTS <- function(x) {
+  axis.POSIXct(side = 1, x = x$DatesR,
+  at = pretty(x$DatesR, n = 10),
+  format = "\%Y-\%m-\%d")
 }
+plot(simGR4J$OutputsModel, Qobs = simGR4J$Qobs, AxisTS = FunAxisTS)
+
+
+### add a main title
 
+## the whole list of settable par's
+opar <- par(no.readonly = TRUE)
+
+## define outer margins and a title inside it
+par(oma = c(0, 0, 3, 0))
+plot(simGR4J$OutputsModel, Qobs = simGR4J$Qobs)
+title(main = "GR4J outputs", outer = TRUE, line = 1.2, cex.main = 1.4)
+
+## reset original par
+par(opar)
+}