diff --git a/DESCRIPTION b/DESCRIPTION
index 52f2699b6f48e6e74d0b482cf32f1625bffa5896..e097255aa8cef88c989aff0049987d064bf1ea43 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: airGR
 Type: Package
 Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.6.9.31
-Date: 2021-01-25
+Version: 1.6.9.32
+Date: 2021-01-26
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
   person("Olivier", "Delaigue", role = c("aut", "cre"), comment = c(ORCID = "0000-0002-7668-8468"), email = "airGR@inrae.fr"),
diff --git a/man/CreateCalibOptions.Rd b/man/CreateCalibOptions.Rd
index b418995a04c37002f14ad8bc3bfb6332167f725f..c98ad748b70f67f09ada0082c21502b5358963ac 100644
--- a/man/CreateCalibOptions.Rd
+++ b/man/CreateCalibOptions.Rd
@@ -80,7 +80,7 @@ CreateCalibOptions(FUN_MOD, FUN_CALIB = Calibration_Michel,
 Users wanting to use \code{FUN_MOD}, \code{FUN_CALIB} or \code{FUN_TRANSFO} functions that are not included in
 the package must create their own \code{CalibOptions} object accordingly. \cr
 
-## ---- CemaNeige version
+## --- CemaNeige version
 
 If \code{IsHyst = FALSE}, the original CemaNeige version from Valéry et al. (2014) is used. \cr
 If \code{IsHyst = TRUE}, the CemaNeige version from Riboust et al. (2019) is used. Compared to the original version, this version of CemaNeige needs two more parameters and it includes a representation of the hysteretic relationship between the Snow Cover Area (SCA) and the Snow Water Equivalent (SWE) in the catchment. The hysteresis included in airGR is the Modified Linear hysteresis (LH*); it is represented on panel b) of Fig. 3 in Riboust et al. (2019). Riboust et al. (2019) advise to use the LH* version of CemaNeige with parameters calibrated using an objective function combining 75 \% of KGE calculated on discharge simulated from a rainfall-runoff model compared to observed discharge and 5 \% of KGE calculated on SCA on 5 CemaNeige elevation bands compared to satellite (e.g. MODIS) SCA (see Eq. (18), Table 3 and Fig. 6). Riboust et al. (2019)'s tests were realized with GR4J as the chosen rainfall-runoff model. \cr
diff --git a/man/CreateInputsCrit.Rd b/man/CreateInputsCrit.Rd
index 848a35fd1d0a01e5eb3b3c4ed64008671e080d00..eed50f03672fc974f63264c022db962ef7a30433 100644
--- a/man/CreateInputsCrit.Rd
+++ b/man/CreateInputsCrit.Rd
@@ -73,11 +73,11 @@ To calculate composite or multiple criteria, it is necessary to use the \code{Er
 \details{
 Users wanting to use \code{FUN_CRIT} functions that are not included in the package must create their own InputsCrit object accordingly. \cr \cr
 
-## ---- Period of calculation
+## --- Period of calculation
 
 Criteria can be calculated over discontinuous periods (i.e. only over winter periods, or when observed discharge is below a certain threshold). To do so, please indicate in \code{Bool_Crit} which indices must be used for calcullation. Discontinuous periods are allowed in the \code{Bool_Crit} argument.
 
-## ---- Transformations
+## --- Transformations
 
 Transformations are simple functions applied to the observed and simulated variables used in order to change their distribution. Transformations are often used in hydrology for modifying the weight put on errors made for high flows or low flows. The following transformations are available: \cr \cr
 \itemize{
@@ -93,11 +93,11 @@ We do not advise computing KGE or KGE' with log-transformation as it might be wr
 In order to make sure that KGE and KGE2 remain dimensionless and are not impacted by zero values, the Box-Cox transformation (\code{transfo = "boxcox"}) uses the formulation given in Equation 10 of Santos et al. (2018). Lambda is set to 0.25 accordingly. \cr \cr
 The syntax of the power transformation allows a numeric or a string of characters. For example for a squared transformation, the following can be used: \code{transfo = 2}, \code{transfo = "2"} or \code{transfo = "^2"}. Negative values are allowed. Fraction values are not allowed (e.g., \code{"-1/2"} must instead be written \code{"-0.5"}).\cr \cr
 
-## ---- The epsilon value
+## --- The epsilon value
 
 The epsilon value is useful when \code{"log"} or \code{"inv"} transformations are used (to avoid calculation of the inverse or of the logarithm of zero). If an epsilon value is provided, then it is added to the observed and simulated variable time series at each time step and before the application of a transformation. The epsilon value has no effect when the \code{"boxcox"} transformation is used. The impact of this value and a recommendation about the epsilon value to use (usually one hundredth of average observation) are discussed in Pushpalatha et al. (2012) for NSE and in Santos et al. (2018) for KGE and KGE'. \cr \cr
 
-## ---- Single, multiple or composite criteria calculation
+## --- Single, multiple or composite criteria calculation
 
 Users can set the following arguments as atomic or list: \code{FUN_CRIT}, \code{Obs}, \code{VarObs}, \code{BoolCrit}, \code{transfo}, \code{Weights}. If the list format is chosen, all the lists must have the same length. \cr
 Calculation of a single criterion (e.g. NSE computed on discharge) is prepared by providing to \code{CreateInputsCrit} arguments atomics only. \cr
diff --git a/man/CreateRunOptions.Rd b/man/CreateRunOptions.Rd
index 6ed188c4670de2f1944cf2dc00ae4af3b979ce2b..20859d4c03a58c843d4c3b7c103d9d0dda56b0db 100644
--- a/man/CreateRunOptions.Rd
+++ b/man/CreateRunOptions.Rd
@@ -72,11 +72,11 @@ CreateRunOptions(FUN_MOD, InputsModel,
 Users wanting to use \code{FUN_MOD} functions that are not included in 
 the package must create their own \code{RunOptions} object accordingly.
 
-## ---- IndPeriod_WarmUp and IndPeriod_Run
+## --- IndPeriod_WarmUp and IndPeriod_Run
 
 Since the hydrological models included in airGR are continuous models, meaning that internal states of the models are propagated to the next time step, \code{IndPeriod_WarmUp} and \code{IndPeriod_Run} must be continuous periods, represented by continuous indices values; no gaps are allowed. To calculate criteria or to calibrate a model over discontinuous periods, please see the \code{Bool_Crit} argument of the \code{\link{CreateInputsCrit}} function.
 
-## ---- Initialisation options
+## --- Initialisation options
 
 The model initialisation options can either be set to a default configuration or be defined by the user.
 
@@ -111,7 +111,7 @@ However, it is also possible to perform a long-term initialisation if other indi
 }
 }
 
-## ---- CemaNeige version
+## --- CemaNeige version
 
 If \code{IsHyst = FALSE}, the original CemaNeige version from Valéry et al. (2014) is used.  \cr
 If \code{IsHyst = TRUE}, the CemaNeige version from Riboust et al. (2019) is used. Compared to the original version, this version of CemaNeige needs two more parameters and it includes a representation of the hysteretic relationship between the Snow Cover Area (SCA) and the Snow Water Equivalent (SWE) in the catchment. The hysteresis included in airGR is the Modified Linear hysteresis (LH*); it is represented on panel b) of Fig. 3 in Riboust et al. (2019). Riboust et al. (2019) advise to use the LH* version of CemaNeige with parameters calibrated using an objective function combining 75 \% of KGE calculated on discharge simulated from a rainfall-runoff model compared to observed discharge and 5 \% of KGE calculated on SCA on 5 CemaNeige elevation bands compared to satellite (e.g. MODIS) SCA (see Eq. (18), Table 3 and Fig. 6). Riboust et al. (2019)'s tests were realized with GR4J as the chosen rainfall-runoff model. \cr
diff --git a/man/Param_Sets_GR4J.Rd b/man/Param_Sets_GR4J.Rd
index 7a3abb7cb111b98e0b79055443cf256aa96a8111..5798b2952e164a8e9ea759dca97d6c8bca928d24 100644
--- a/man/Param_Sets_GR4J.Rd
+++ b/man/Param_Sets_GR4J.Rd
@@ -56,7 +56,7 @@ Param_Sets_GR4J <- as.matrix(Param_Sets_GR4J)
 InputsModel <- CreateInputsModel(FUN_MOD = RunModel_GR4J, DatesR = BasinObs$DatesR,
                                  Precip = BasinObs$P, PotEvap = BasinObs$E)
 
-## ---- calibration step
+## --- calibration step
 
 ## short calibration period selection (< 6 months)
 Ind_Cal <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-01-01"),
@@ -82,7 +82,7 @@ OutputsCrit_Loop <- apply(Param_Sets_GR4J, 1, function(Param) {
 Param_Best <- unlist(Param_Sets_GR4J[which.max(OutputsCrit_Loop), ])
 
 
-## ---- validation step
+## --- validation step
 
 ## validation period selection
 Ind_Val <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-03-01"),
diff --git a/man/RunModel_CemaNeige.Rd b/man/RunModel_CemaNeige.Rd
index 2fc511d5641a82f178833e331779eb2b9cb2ff2a..685cb2a6d5eacb7961b52a0a4fdb3645e30cfb85 100644
--- a/man/RunModel_CemaNeige.Rd
+++ b/man/RunModel_CemaNeige.Rd
@@ -79,7 +79,7 @@ Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-01-0
                which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1999-12-31"))
 
 
-## ---- original version of CemaNeige
+## --- original version of CemaNeige
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeige, InputsModel = InputsModel,
@@ -94,7 +94,7 @@ OutputsModel <- RunModel_CemaNeige(InputsModel = InputsModel,
 plot(OutputsModel)
 
 
-## ---- version of CemaNeige with the Linear Hysteresis
+## --- version of CemaNeige with the Linear Hysteresis
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeige, InputsModel = InputsModel,
diff --git a/man/RunModel_CemaNeigeGR4H.Rd b/man/RunModel_CemaNeigeGR4H.Rd
index 7dbe52223adbac99b1419c6f246fbe48f1cded56..d2a499da73a6f7faf9e7b8fbff32f9ac4b3ec7df 100644
--- a/man/RunModel_CemaNeigeGR4H.Rd
+++ b/man/RunModel_CemaNeigeGR4H.Rd
@@ -108,7 +108,7 @@ Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d \%H:\%M")=="2
                which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d \%H:\%M")=="2008-12-31 23:00"))
 
 
-## ---- original version of CemaNeige
+## --- original version of CemaNeige
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeigeGR4H, InputsModel = InputsModel,
diff --git a/man/RunModel_CemaNeigeGR4J.Rd b/man/RunModel_CemaNeigeGR4J.Rd
index c451c6a3f2be63b5ba5c29736b287358fe363b01..9ab766451012e0d66c42401884b4f54c795fd080 100644
--- a/man/RunModel_CemaNeigeGR4J.Rd
+++ b/man/RunModel_CemaNeigeGR4J.Rd
@@ -104,7 +104,7 @@ Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-01-0
                which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1999-12-31"))
 
 
-## ---- original version of CemaNeige
+## --- original version of CemaNeige
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeigeGR4J, InputsModel = InputsModel,
@@ -125,7 +125,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 
 
-## ---- version of CemaNeige with the Linear Hysteresis
+## --- version of CemaNeige with the Linear Hysteresis
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeigeGR4J, InputsModel = InputsModel,
diff --git a/man/RunModel_CemaNeigeGR5H.Rd b/man/RunModel_CemaNeigeGR5H.Rd
index 5b9f07e5e0008631fe6a381e8f23bcf6faff6861..e7c3f8a9826af875b8064ba7884800926204410b 100644
--- a/man/RunModel_CemaNeigeGR5H.Rd
+++ b/man/RunModel_CemaNeigeGR5H.Rd
@@ -109,7 +109,7 @@ Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d \%H:\%M")=="2
                which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d \%H:\%M")=="2008-12-31 23:00"))
 
 
-## ---- original version of CemaNeige
+## --- original version of CemaNeige
 
 ## Imax computation
 Imax <- Imax(InputsModel = InputsModel, IndPeriod_Run = Ind_Run,
diff --git a/man/RunModel_CemaNeigeGR6J.Rd b/man/RunModel_CemaNeigeGR6J.Rd
index 65e9635d50cf00bc89d0421c6830654f867ddddb..e9a37d4de87e101a6135302f1d0d97e33dfd808c 100644
--- a/man/RunModel_CemaNeigeGR6J.Rd
+++ b/man/RunModel_CemaNeigeGR6J.Rd
@@ -108,7 +108,7 @@ Ind_Run <- seq(which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1990-01-0
                which(format(BasinObs$DatesR, format = "\%Y-\%m-\%d")=="1999-12-31"))
 
 
-## ---- original version of CemaNeige
+## --- original version of CemaNeige
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeigeGR6J, InputsModel = InputsModel,
@@ -129,7 +129,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 
 
-## ---- version of CemaNeige with the Linear Hysteresis
+## --- version of CemaNeige with the Linear Hysteresis
 
 ## preparation of the RunOptions object
 RunOptions <- CreateRunOptions(FUN_MOD = RunModel_CemaNeigeGR6J, InputsModel = InputsModel,
diff --git a/man/TransfoParam.Rd b/man/TransfoParam.Rd
index de58b74b8908342e621d973d77e5e7eaef589440..2029ed45c330e662c20d34a163f33e9a3febfe49 100644
--- a/man/TransfoParam.Rd
+++ b/man/TransfoParam.Rd
@@ -55,7 +55,7 @@ TransfoParam_CemaNeigeHyst(ParamIn, Direction)
 \examples{
 library(airGR)
 
-## ---- generic function
+## --- generic function
 
 ## transformation Raw -> Transformed for the GR4J model
 Xraw  <- matrix(c(+221.41, -3.63,  +30.00, +1.37,
@@ -72,7 +72,7 @@ Xtran <- matrix(c(+3.60, -2.00, +3.40, -9.10,
 Xraw  <- TransfoParam(ParamIn = Xtran, Direction = "TR", FUN_TRANSFO = TransfoParam_GR4J)
 
 
-## ---- specific function
+## --- specific function
 
 ## transformation Raw -> Transformed for the GR4J model
 Xraw  <- matrix(c(+221.41, -3.63,  +30.00, +1.37,