diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index 139db9c66a3518dbda6c9659506721a009d5ce8e..e90893e844e192f2b05bcf629954209b39161eb0 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -19,10 +19,11 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
 
   ## index time series
   if (!is.null(IndPeriod_Plot)) {
-    if (length(IndPeriod_Plot) == 0) {
+    if (length(IndPeriod_Plot) <= 1 || any(IndPeriod_Plot <= 0)) {
       IndPeriod_Plot <- seq_along(OutputsModel$DatesR)
+      stop('IndPeriod_Plot must be a positive vector of length > 0')
     }
-    IndPeriod_Plot <- seq_along(IndPeriod_Plot)
+    # IndPeriod_Plot <- seq_along(IndPeriod_Plot)
     OutputsModel <- .IndexOutputsModel(OutputsModel, IndPeriod_Plot)
     Qobs <- Qobs[IndPeriod_Plot]
   } else {