diff --git a/DESCRIPTION b/DESCRIPTION
index b769dd27f55f31c28cd0004384bd59b87fe68405..e096b0f588bcc63a36fe0adbbc84d2dc3d18e094 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.8.5
-Date: 2017-06-21
+Version: 1.0.8.6
+Date: 2017-06-22
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl")),
   person("Charles", "Perrin", role = c("aut", "ths")),
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index d664ec414be42617eed4956174cabb5827104e02..4a49ccfc2f9c9c09b1b6547e243b5d382fb71152 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -436,6 +436,7 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
                                 ylim <- range(log(OutputsModel$Qsim[IndPeriod_Plot][SelectNotZero]), na.rm = TRUE); }
     if ( BOOL_Qobs &  BOOL_Qsim) { SelectNotZero <- SelectQobsNotZero & SelectQsimNotZero;
                                 ylim <- range(log(c(Qobs[IndPeriod_Plot][SelectNotZero], OutputsModel$Qsim[IndPeriod_Plot][SelectNotZero])), na.rm = TRUE); }
+    SelectNotZero <- ifelse(is.na(SelectNotZero), FALSE, SelectNotZero)
     if (any(SelectNotZero)) {
       plot(0, 0, type = "n",
            xlim = xlim, ylim = ylim,