diff --git a/R/CreateInputsModel.GRiwrm.R b/R/CreateInputsModel.GRiwrm.R
index e36d26bea6eb412ce1fb6c0e44b683e6ced6c946..e9f560c18c8343f7e99b32c3828bebd9e4785f36 100644
--- a/R/CreateInputsModel.GRiwrm.R
+++ b/R/CreateInputsModel.GRiwrm.R
@@ -215,7 +215,8 @@ CreateInputsModel.GRiwrm <- function(x, DatesR,
                                  HypsoData = getInputBV(HypsoData, id),
                                  NLayers = getInputBV(NLayers, id, 5),
                                  Qobs = Qobs,
-                                 Qmin = getInputBV(Qmin, id)
+                                 Qmin = getInputBV(Qmin, id),
+                                 IsHyst = IsHyst
                                  )
   }
   attr(InputsModel, "TimeStep") <- getModelTimeStep(InputsModel)
@@ -250,7 +251,7 @@ CreateEmptyGRiwrmInputsModel <- function(griwrm) {
 #'
 #' @return \emph{InputsModel} object for one.
 #' @noRd
-CreateOneGRiwrmInputsModel <- function(id, griwrm, ..., Qobs, Qmin) {
+CreateOneGRiwrmInputsModel <- function(id, griwrm, ..., Qobs, Qmin, IsHyst) {
   np <- getNodeProperties(id, griwrm)
 
   if (np$Diversion) {