diff --git a/DESCRIPTION b/DESCRIPTION
index 8bd17dd685ca22212bf6aa71ced429e244b766d7..5b92168bc0c9d0cfc6812a6d1c2f99dbf0bb761c 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.0.9.6
-Date: 2017-06-26
+Version: 1.0.9.7
+Date: 2017-07-11
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl")),
   person("Charles", "Perrin", role = c("aut", "ths")),
diff --git a/man/RunModel.Rd b/man/RunModel.Rd
index 88ec5aec3e8a1ad086a91a8c587a82353a49292e..16c22d31df3161693d4aa8bb293f1e2ee951486e 100644
--- a/man/RunModel.Rd
+++ b/man/RunModel.Rd
@@ -52,9 +52,9 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron (June 2014)
+Laurent Coron
 }
 \seealso{
-\code{\link{RunModel_GR4J}}, \code{\link{RunModel_CemaNeigeGR4J}}, \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+\code{\link{RunModel_GR4J}}, \code{\link{RunModel_CemaNeigeGR4J}}, \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_CemaNeige.Rd b/man/RunModel_CemaNeige.Rd
index f77a57944803eba28683f2f41466d75dad11581e..a8394a82b5c7ac54d4da9b0e0f687c75a70163a4 100644
--- a/man/RunModel_CemaNeige.Rd
+++ b/man/RunModel_CemaNeige.Rd
@@ -70,7 +70,7 @@ plot(OutputsModel)
 
 }
 \author{
-Laurent Coron (January 2014)
+Laurent Coron
 }
 \references{
 Valéry, A., V. Andréassian and C. Perrin (2014), 
@@ -81,6 +81,6 @@ Valéry, A., V. Andréassian and C. Perrin (2014),
       Part 2 - Sensitivity analysis of the Cemaneige snow accounting routine on 380 catchments, Journal of Hydrology, doi:10.1016/j.jhydrol.2014.04.058.
 }
 \seealso{
-\code{\link{RunModel_CemaNeigeGR4J}}, \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+\code{\link{RunModel_CemaNeigeGR4J}}, \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_CemaNeigeGR4J.Rd b/man/RunModel_CemaNeigeGR4J.Rd
index 63993f9eeb2f48687b9658eca26e32defa39ba8e..a0e5fb86c226a9170e60698f4085fec8e4a89f61 100644
--- a/man/RunModel_CemaNeigeGR4J.Rd
+++ b/man/RunModel_CemaNeigeGR4J.Rd
@@ -28,10 +28,10 @@ CemaNeige X2 \tab degree-day melt coefficient [mm/°C/d]                       \
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/d]                          \cr
          \emph{$Prod    }          \tab [numeric] series of production store level [mm]                        \cr
          \emph{$Pn      }          \tab [numeric] series of net rainfall [mm/d]                         			  \cr
-         \emph{$Ps      }          \tab [numeric] series of the part of Ps filling the production store [mm/d]        \cr
+         \emph{$Ps      }          \tab [numeric] series of the part of Pn filling the production store [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 [mm/d]                                      \cr
+         \emph{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [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 [mm]                           \cr
@@ -41,7 +41,7 @@ CemaNeige X2 \tab degree-day melt coefficient [mm/°C/d]                       \
          \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 Qsim [mm/d]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/d]                                               \cr
          \emph{$CemaNeigeLayers}   \tab [list] list of CemaNeige outputs (1 list per layer)                          \cr
          \emph{$CemaNeigeLayers[[iLayer]]$Pliq         }   \tab [numeric] series of liquid precip. [mm/d]                          \cr
          \emph{$CemaNeigeLayers[[iLayer]]$Psol         }   \tab [numeric] series of solid precip. [mm/d]                           \cr
@@ -97,7 +97,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Audrey Valéry, Claude Michel, Charles Perrin, Vazken Andréassian (December 2013)
+Laurent Coron, Audrey Valéry, Claude Michel, Charles Perrin, Vazken Andréassian
 }
 \references{
 Perrin, C., C. Michel and V. Andréassian (2003), 
@@ -112,6 +112,6 @@ Perrin, C., C. Michel and V. Andréassian (2003),
 }
 \seealso{
 \code{\link{RunModel_CemaNeigeGR5J}}, \code{\link{RunModel_CemaNeigeGR6J}}, \code{\link{RunModel_GR4J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_CemaNeigeGR5J.Rd b/man/RunModel_CemaNeigeGR5J.Rd
index 6ccb4e2d38a0e546a96a34e71ba8abd39ff020c3..46866ea5cedb24dc9182ea8526874f468c3dfd98 100644
--- a/man/RunModel_CemaNeigeGR5J.Rd
+++ b/man/RunModel_CemaNeigeGR5J.Rd
@@ -29,10 +29,10 @@ CemaNeige X2 \tab degree-day melt coefficient [mm/°C/d]                       \
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/d]                          \cr
          \emph{$Prod    }          \tab [numeric] series of production store level [mm]                        \cr
          \emph{$Pn      }          \tab [numeric] series of net rainfall [mm/d]                         			  \cr
-         \emph{$Ps      }          \tab [numeric] series of the part of Ps filling the production store [mm/d]        \cr
+         \emph{$Ps      }          \tab [numeric] series of the part of Pn filling the production store [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 [mm/d]                                      \cr
+         \emph{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [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 [mm]                           \cr
@@ -42,7 +42,7 @@ CemaNeige X2 \tab degree-day melt coefficient [mm/°C/d]                       \
          \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 Qsim [mm/d]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/d]                                               \cr
          \emph{$CemaNeigeLayers}   \tab [list] list of CemaNeige outputs (1 list per layer)                           \cr
          \emph{$CemaNeigeLayers[[iLayer]]$Pliq         }   \tab [numeric] series of liquid precip. [mm/d]                          \cr
          \emph{$CemaNeigeLayers[[iLayer]]$Psol         }   \tab [numeric] series of solid precip. [mm/d]                           \cr
@@ -98,7 +98,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Audrey Valéry, Claude Michel, Nicolas Le Moine, Charles Perrin, Vazken Andréassian (December 2013)
+Laurent Coron, Audrey Valéry, Claude Michel, Nicolas Le Moine, Charles Perrin, Vazken Andréassian
 }
 \references{
 Le Moine, N. (2008), Le bassin versant de surface vu par le souterrain : une voie d'amélioration des performances 
@@ -115,6 +115,6 @@ Le Moine, N. (2008), Le bassin versant de surface vu par le souterrain : une voi
 }
 \seealso{
 \code{\link{RunModel_CemaNeigeGR4J}}, \code{\link{RunModel_CemaNeigeGR6J}}, \code{\link{RunModel_GR5J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_GR1A.Rd b/man/RunModel_GR1A.Rd
index eb76db2e5bae2b94212f5e048c898aae1f20947c..f34bbb8c4b87740461acb7ea7c85523059d854b7 100644
--- a/man/RunModel_GR1A.Rd
+++ b/man/RunModel_GR1A.Rd
@@ -21,7 +21,7 @@ GR1A X1      \tab model parameter [-] \cr
          \emph{$DatesR  }          \tab [POSIXlt] series of dates                                                     \cr
          \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 Qsim [mm/y]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/y]                                               \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)
@@ -76,11 +76,14 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Claude Michel (March 2015)
+Laurent Coron, Claude Michel
 }
 \references{
 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, France. \cr
 }
+\seealso{
+\code\code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+}
 
diff --git a/man/RunModel_GR2M.Rd b/man/RunModel_GR2M.Rd
index f697c347f3dc18eab33324da254d57c5b957ff01..3199ab2716bc4cc8057c43acb8d0c2b2edd1394a 100644
--- a/man/RunModel_GR2M.Rd
+++ b/man/RunModel_GR2M.Rd
@@ -22,7 +22,7 @@ GR2M X2      \tab groundwater exchange coefficient [-]    \cr
          \emph{$DatesR  }          \tab [POSIXlt] series of dates                                                          \cr
          \emph{$PotEvap }          \tab [numeric] series of input potential evapotranspiration [mm/month]                  \cr
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/month]                           \cr
-         \emph{$Qsim    }          \tab [numeric] series of Qsim [mm/month]                                                \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/month]                                                \cr
          \emph{$StateEnd}          \tab [numeric] states at the end of the run (production store level and routing store level) [mm], \cr\tab 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)
@@ -75,7 +75,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Claude Michel, Safouane Mouelhi (March 2015)
+Laurent Coron, Claude Michel, Safouane Mouelhi
 }
 \references{
 Mouelhi S. (2003), 
@@ -85,4 +85,6 @@ Mouelhi S. (2003),
       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.
 }
-
+\seealso{
+\code{\code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
+}
diff --git a/man/RunModel_GR4H.Rd b/man/RunModel_GR4H.Rd
index 07a9e7dc05e2082742e98ac8bde0d098652a7803..16287a3318b07f8974fc85a3faf180e8282f4b09 100644
--- a/man/RunModel_GR4H.Rd
+++ b/man/RunModel_GR4H.Rd
@@ -27,7 +27,7 @@ GR4H X4      \tab unit hydrograph time constant [h]
          \emph{$Prod    }          \tab [numeric] series of production store level [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{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [mm/h]                                      \cr
          \emph{$Q9      }          \tab [numeric] series of UH1 outflow (Q9) [mm/h]                                   \cr
          \emph{$Q1      }          \tab [numeric] series of UH2 outflow (Q1) [mm/h]                                   \cr
          \emph{$Rout    }          \tab [numeric] series of routing store level [mm]                           \cr
@@ -35,7 +35,7 @@ GR4H X4      \tab unit hydrograph time constant [h]
          \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 UH2 after exchange (QD) [mm/h]           \cr
-         \emph{$Qsim    }          \tab [numeric] series of Qsim [mm/h]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/h]                                               \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)
@@ -88,6 +88,6 @@ Le Moine, N. (2008), Le bassin versant de surface vu par le souterrain : une voi
 }
 \seealso{
 \code{\link{RunModel_GR4J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_GR4J.Rd b/man/RunModel_GR4J.Rd
index 924b112ff522e444372e2b2cfa892935f786a38d..da94a292662585ee9b6d86a1a00e1abc26584195 100644
--- a/man/RunModel_GR4J.Rd
+++ b/man/RunModel_GR4J.Rd
@@ -26,10 +26,10 @@ GR4J X4      \tab unit hydrograph time constant [d]
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/d]                          \cr
          \emph{$Prod    }          \tab [numeric] series of production store level [mm]                 	          \cr
          \emph{$Pn      }          \tab [numeric] series of net rainfall [mm/d]                         			  \cr
-         \emph{$Ps      }          \tab [numeric] series of the part of Ps filling the production store [mm/d]        \cr
+         \emph{$Ps      }          \tab [numeric] series of the part of Pn filling the production store [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 [mm/d]                                      \cr
+         \emph{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [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 [mm]                     		      \cr
@@ -39,7 +39,7 @@ GR4J X4      \tab unit hydrograph time constant [d]
          \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 Qsim [mm/d]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/d]                                               \cr
          \emph{$StateEnd}          \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm], \cr\tab 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)
@@ -82,7 +82,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Claude Michel, Charles Perrin (December 2013)
+Laurent Coron, Claude Michel, Charles Perrin
 }
 \references{
 Perrin, C., C. Michel and V. Andréassian (2003), 
@@ -91,6 +91,6 @@ Perrin, C., C. Michel and V. Andréassian (2003),
 }
 \seealso{
 \code{\link{RunModel_GR5J}}, \code{\link{RunModel_GR6J}}, \code{\link{RunModel_CemaNeigeGR4J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_GR5J.Rd b/man/RunModel_GR5J.Rd
index 94f3f4863906929672db4f92c3498515ed02c3ae..e27c1813379cddd2658fa428839aff0bdd380682 100644
--- a/man/RunModel_GR5J.Rd
+++ b/man/RunModel_GR5J.Rd
@@ -27,10 +27,10 @@ GR5J X5      \tab intercatchment exchange threshold [-]                     \cr
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/d]                          \cr
          \emph{$Prod    }          \tab [numeric] series of production store level [mm]                        \cr
          \emph{$Pn      }          \tab [numeric] series of net rainfall [mm/d]                         			  \cr
-         \emph{$Ps      }          \tab [numeric] series of the part of Ps filling the production store [mm/d]        \cr
+         \emph{$Ps      }          \tab [numeric] series of the part of Pn filling the production store [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 [mm/d]                                      \cr
+         \emph{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [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 [mm]                           \cr
@@ -40,7 +40,7 @@ GR5J X5      \tab intercatchment exchange threshold [-]                     \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 Qsim [mm/d]                                               \cr
+         \emph{$Qsim    }          \tab [numeric] series of simulated discharge [mm/d]                                               \cr
          \emph{$StateEnd}          \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm], \cr\tab 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)
@@ -84,7 +84,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Claude Michel, Nicolas Le Moine (December 2013)
+Laurent Coron, Claude Michel, Nicolas Le Moine
 }
 \references{
 Le Moine, N. (2008), Le bassin versant de surface vu par le souterrain : une voie d'amélioration des performances
@@ -95,6 +95,6 @@ Pushpalatha, R., C. Perrin, N. Le Moine, T. Mathevet, and V. Andréassian (2011)
 }
 \seealso{
 \code{\link{RunModel_GR4J}}, \code{\link{RunModel_GR6J}}, \code{\link{RunModel_CemaNeigeGR5J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }
 
diff --git a/man/RunModel_GR6J.Rd b/man/RunModel_GR6J.Rd
index 60c30a1d73f06678b2c2315fec372d73b129ac9b..c36bedd5298bb99b46360312e78dda6a88e0aa1f 100644
--- a/man/RunModel_GR6J.Rd
+++ b/man/RunModel_GR6J.Rd
@@ -28,10 +28,10 @@ GR6J X6      \tab coefficient for emptying exponential store [mm]
          \emph{$Precip  }          \tab [numeric] series of input total precipitation [mm/d]                          \cr
          \emph{$Prod    }          \tab [numeric] series of production store level [mm]                        \cr
          \emph{$Pn      }          \tab [numeric] series of net rainfall [mm/d]                         			  \cr
-         \emph{$Ps      }          \tab [numeric] series of the part of Ps filling the production store [mm/d]        \cr
+         \emph{$Ps      }          \tab [numeric] series of the part of Pn filling the production store [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 [mm/d]                                      \cr
+         \emph{$PR      }          \tab [numeric] series of PR=Pn-Ps+Perc [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 [mm]                           \cr
@@ -86,7 +86,7 @@ InputsCrit  <- CreateInputsCrit(FUN_CRIT = ErrorCrit_NSE, InputsModel = InputsMo
 OutputsCrit <- ErrorCrit_NSE(InputsCrit = InputsCrit, OutputsModel = OutputsModel)
 }
 \author{
-Laurent Coron, Claude Michel, Charles Perrin, Raji Pushpalatha, Nicolas Le Moine (December 2013)
+Laurent Coron, Claude Michel, Charles Perrin, Raji Pushpalatha, Nicolas Le Moine
 }
 \references{
 Pushpalatha, R., C. Perrin, N. Le Moine, T. Mathevet and V. Andréassian (2011), 
@@ -95,6 +95,6 @@ Pushpalatha, R., C. Perrin, N. Le Moine, T. Mathevet and V. Andréassian (2011),
 }
 \seealso{
 \code{\link{RunModel_GR4J}}, \code{\link{RunModel_GR5J}}, \code{\link{RunModel_CemaNeigeGR6J}},
-         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}.
+         \code{\link{CreateInputsModel}}, \code{\link{CreateRunOptions}}, \code{\link{CreateIniStates}}.
 }