diff --git a/man/RunModel_CemaNeigeGR4H.Rd b/man/RunModel_CemaNeigeGR4H.Rd
index d2a499da73a6f7faf9e7b8fbff32f9ac4b3ec7df..a5f77e0fb66dc00df797a9af9ab69acdacc86566 100644
--- a/man/RunModel_CemaNeigeGR4H.Rd
+++ b/man/RunModel_CemaNeigeGR4H.Rd
@@ -71,6 +71,8 @@ RunModel_CemaNeigeGR4H(InputsModel, RunOptions, Param)
     \emph{$CemaNeigeLayers[[iLayer]]$Temp        } \tab [numeric] series of air temperature [°C]                     \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Gthreshold  } \tab [numeric] series of melt threshold [mm]                      \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Glocalmax   } \tab [numeric] series of local melt threshold for hysteresis [mm] \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/h] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run: store & unit hydrographs levels [mm], CemaNeige states [mm & °C]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_CemaNeigeGR4J.Rd b/man/RunModel_CemaNeigeGR4J.Rd
index 9ab766451012e0d66c42401884b4f54c795fd080..3cb311a3bd4f1353f1933fd992024841c842a80d 100644
--- a/man/RunModel_CemaNeigeGR4J.Rd
+++ b/man/RunModel_CemaNeigeGR4J.Rd
@@ -71,6 +71,8 @@ RunModel_CemaNeigeGR4J(InputsModel, RunOptions, Param)
     \emph{$CemaNeigeLayers[[iLayer]]$Temp         } \tab [numeric] series of air temperature [°C]                     \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Gthreshold   } \tab [numeric] series of melt threshold [mm]                      \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Glocalmax    } \tab [numeric] series of local melt threshold for hysteresis [mm] \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run: store & unit hydrographs levels [mm], CemaNeige states [mm & °C]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_CemaNeigeGR5H.Rd b/man/RunModel_CemaNeigeGR5H.Rd
index e7c3f8a9826af875b8064ba7884800926204410b..3c45c615623f3b5b538bdce637c9a2ba004e7712 100644
--- a/man/RunModel_CemaNeigeGR5H.Rd
+++ b/man/RunModel_CemaNeigeGR5H.Rd
@@ -74,6 +74,8 @@ RunModel_CemaNeigeGR5H(InputsModel, RunOptions, Param)
     \emph{$CemaNeigeLayers[[iLayer]]$Temp         } \tab [numeric] series of air temperature [°C]                               \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Gthreshold   } \tab [numeric] series of melt threshold [mm]                                \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Glocalmax    } \tab [numeric] series of local melt threshold for hysteresis [mm]           \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/h] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run: store & unit hydrographs levels [mm], CemaNeige states [mm & °C]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_CemaNeigeGR5J.Rd b/man/RunModel_CemaNeigeGR5J.Rd
index 4f2b77c50e6321acbd25936b86828bbd030dfa52..2ca146db3ab7286d8714667dddb5131ecdd74e5b 100644
--- a/man/RunModel_CemaNeigeGR5J.Rd
+++ b/man/RunModel_CemaNeigeGR5J.Rd
@@ -72,6 +72,8 @@ RunModel_CemaNeigeGR5J(InputsModel, RunOptions, Param)
     \emph{$CemaNeigeLayers[[iLayer]]$Temp        } \tab [numeric] series of air temperature [°C]                     \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Gthreshold  } \tab [numeric] series of melt threshold [mm]                      \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Glocalmax   } \tab [numeric] series of local melt threshold for hysteresis [mm] \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run: store & unit hydrographs levels [mm], CemaNeige states [mm & °C]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs
diff --git a/man/RunModel_CemaNeigeGR6J.Rd b/man/RunModel_CemaNeigeGR6J.Rd
index e9a37d4de87e101a6135302f1d0d97e33dfd808c..84c83e82e6e3ee8e11ea97b2aa637d6e84288a03 100644
--- a/man/RunModel_CemaNeigeGR6J.Rd
+++ b/man/RunModel_CemaNeigeGR6J.Rd
@@ -75,6 +75,8 @@ RunModel_CemaNeigeGR6J(InputsModel, RunOptions, Param)
     \emph{$CemaNeigeLayers[[iLayer]]$Temp        } \tab [numeric] series of air temperature [°C]                     \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Gthreshold  } \tab [numeric] series of melt threshold [mm]                      \cr
     \emph{$CemaNeigeLayers[[iLayer]]$Glocalmax   } \tab [numeric] series of local melt threshold for hysteresis [mm] \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run: store & unit hydrographs levels [mm], CemaNeige states [mm & °C]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs
diff --git a/man/RunModel_GR1A.Rd b/man/RunModel_GR1A.Rd
index 0ebba8426ced4aa870354ba128ca21c9d2839e11..a35613a3e6d8b6a5e13072221a22a2bb33c573b9 100644
--- a/man/RunModel_GR1A.Rd
+++ b/man/RunModel_GR1A.Rd
@@ -37,6 +37,8 @@ RunModel_GR1A(InputsModel, RunOptions, Param)
     \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/y] \cr
     \emph{$Precip  } \tab [numeric] series of input total precipitation [mm/y]          \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge [mm/y]                \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/y] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (NULL) [-]             \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR2M.Rd b/man/RunModel_GR2M.Rd
index 11b3f2ce58a89f8e5339abb80823728bc607f36d..2aff86417d3e32b5345adab8829fee3c28dfa4ec 100644
--- a/man/RunModel_GR2M.Rd
+++ b/man/RunModel_GR2M.Rd
@@ -46,6 +46,8 @@ RunModel_GR2M(InputsModel, RunOptions, Param)
     \emph{$Prod    } \tab [numeric] series of production store level (S) [mm]                   \cr
     \emph{$Rout    } \tab [numeric] series of routing store level (R1) [mm]                     \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge [mm/month] (Q)                \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/month] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                      \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (production store level and routing store level) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR4H.Rd b/man/RunModel_GR4H.Rd
index 7f2b2885e22d37f6f3e5af708cff925f6f39701c..ef0d81ea69032d7700128e0cf2eb9cab034f4669 100644
--- a/man/RunModel_GR4H.Rd
+++ b/man/RunModel_GR4H.Rd
@@ -55,6 +55,8 @@ RunModel_GR4H(InputsModel, RunOptions, Param)
     \emph{$QR      } \tab [numeric] series of routing store outflow (Qr) [mm/h]                         \cr
     \emph{$QD      } \tab [numeric] series of direct flow from UH2 after exchange (Qd) [mm/h]           \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/h]                            \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/h] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR4J.Rd b/man/RunModel_GR4J.Rd
index 0536238b05b88bb04ebea47d91f88e52dae3dcdf..7c97e6c209cc91a1aef4f04dee86479dac21527f 100644
--- a/man/RunModel_GR4J.Rd
+++ b/man/RunModel_GR4J.Rd
@@ -36,25 +36,27 @@ RunModel_GR4J(InputsModel, RunOptions, Param)
 \value{
 [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/d] (E)                  \cr
-    \emph{$Precip  } \tab [numeric] series of input total precipitation (P) [mm/d]                           \cr
-    \emph{$Prod    } \tab [numeric] series of production store level (S) [mm]                                \cr
-    \emph{$Pn      } \tab [numeric] series of net rainfall (Pn) [mm/d]                                       \cr
-    \emph{$Ps      } \tab [numeric] series of the part of Pn filling the production store (Ps) [mm/d]        \cr
-    \emph{$AE      } \tab [numeric] series of actual evapotranspiration [mm/d]                               \cr
-    \emph{$Perc    } \tab [numeric] series of percolation (Perc) [mm/d]                                      \cr
-    \emph{$PR      } \tab [numeric] series of Pr=Pn-Ps+Perc (Pr) [mm/d]                                      \cr
-    \emph{$Q9      } \tab [numeric] series of UH1 outflow (Q9) [mm/d]                                        \cr
-    \emph{$Q1      } \tab [numeric] series of UH2 outflow (Q1) [mm/d]                                        \cr
-    \emph{$Rout    } \tab [numeric] series of routing store level (R1) [mm]                                  \cr
-    \emph{$Exch    } \tab [numeric] series of potential semi-exchange between catchments [mm/d]              \cr
-    \emph{$AExch1  } \tab [numeric] series of actual exchange between catchments for branch 1 [mm/d]         \cr
-    \emph{$AExch2  } \tab [numeric] series of actual exchange between catchments for branch 2 [mm/d]         \cr
-    \emph{$AExch   } \tab [numeric] series of actual exchange between catchments (1+2) [mm/d]                \cr
-    \emph{$QR      } \tab [numeric] series of routing store outflow (Qr) [mm/d]                              \cr
-    \emph{$QD      } \tab [numeric] series of direct flow from UH2 after exchange (Qd) [mm/d]                \cr
-    \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/d]                                 \cr
+    \emph{$DatesR  } \tab [POSIXlt] series of dates                                                   \cr
+    \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/d] (E)           \cr
+    \emph{$Precip  } \tab [numeric] series of input total precipitation (P) [mm/d]                    \cr
+    \emph{$Prod    } \tab [numeric] series of production store level (S) [mm]                         \cr
+    \emph{$Pn      } \tab [numeric] series of net rainfall (Pn) [mm/d]                                \cr
+    \emph{$Ps      } \tab [numeric] series of the part of Pn filling the production store (Ps) [mm/d] \cr
+    \emph{$AE      } \tab [numeric] series of actual evapotranspiration [mm/d]                        \cr
+    \emph{$Perc    } \tab [numeric] series of percolation (Perc) [mm/d]                               \cr
+    \emph{$PR      } \tab [numeric] series of Pr=Pn-Ps+Perc (Pr) [mm/d]                               \cr
+    \emph{$Q9      } \tab [numeric] series of UH1 outflow (Q9) [mm/d]                                 \cr
+    \emph{$Q1      } \tab [numeric] series of UH2 outflow (Q1) [mm/d]                                 \cr
+    \emph{$Rout    } \tab [numeric] series of routing store level (R1) [mm]                           \cr
+    \emph{$Exch    } \tab [numeric] series of potential semi-exchange between catchments [mm/d]       \cr
+    \emph{$AExch1  } \tab [numeric] series of actual exchange between catchments for branch 1 [mm/d]  \cr
+    \emph{$AExch2  } \tab [numeric] series of actual exchange between catchments for branch 2 [mm/d]  \cr
+    \emph{$AExch   } \tab [numeric] series of actual exchange between catchments (1+2) [mm/d]         \cr
+    \emph{$QR      } \tab [numeric] series of routing store outflow (Qr) [mm/d]                       \cr
+    \emph{$QD      } \tab [numeric] series of direct flow from UH2 after exchange (Qd) [mm/d]         \cr
+    \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/d]                          \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR5H.Rd b/man/RunModel_GR5H.Rd
index 997566b79e7931a97801fcfcdaef71421e549362..f1f782e8c7eedbb5fea3922e6d2863080fdf6e8e 100644
--- a/man/RunModel_GR5H.Rd
+++ b/man/RunModel_GR5H.Rd
@@ -59,6 +59,8 @@ RunModel_GR5H(InputsModel, RunOptions, Param)
     \emph{$QR      } \tab [numeric] series of routing store outflow (Qr) [mm/h]                                                 \cr
     \emph{$QD      } \tab [numeric] series of direct flow from UH after exchange (Qd) [mm/h]                                    \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/h]                                                    \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/h]                  \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                                   \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH levels) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR5J.Rd b/man/RunModel_GR5J.Rd
index 227d86182fc22d8395352d37d636906873585f0e..755207efc7cbc422d1895f472cec0917081c965d 100644
--- a/man/RunModel_GR5J.Rd
+++ b/man/RunModel_GR5J.Rd
@@ -56,6 +56,8 @@ RunModel_GR5J(InputsModel, RunOptions, Param)
     \emph{$QR      } \tab [numeric] series of routing store outflow (Qr) [mm/d]                         \cr
     \emph{$QD      } \tab [numeric] series of direct flow from UH after exchange (Qd) [mm/d]            \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/d]                            \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH levels) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.
diff --git a/man/RunModel_GR6J.Rd b/man/RunModel_GR6J.Rd
index 7a1849b5ec58407b9c9ee01a645307c65dd4e2e7..fc510e7a028f89f3282920cc87b2ade29dbde5af 100644
--- a/man/RunModel_GR6J.Rd
+++ b/man/RunModel_GR6J.Rd
@@ -59,6 +59,8 @@ RunModel_GR6J(InputsModel, RunOptions, Param)
     \emph{$Exp     } \tab [numeric] series of exponential store level (negative) (R2) [mm]              \cr
     \emph{$QD      } \tab [numeric] series of direct flow from UH2 after exchange (Qd) [mm/d]           \cr
     \emph{$Qsim    } \tab [numeric] series of simulated discharge (Q) [mm/d]                            \cr
+    \emph{RunOptions$WarmUpQsim} \tab [numeric] series of simulated discharge (Q) on the warm-up period [mm/d] \cr
+    \emph{RunOptions$Param     } \tab [numeric] parameter set parameter set used by the model                  \cr
     \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm]. See \code{\link{CreateIniStates}} for more details \cr
   }
 Refer to the provided references or to the package source code for further details on these model outputs.