diff --git a/DESCRIPTION b/DESCRIPTION
index b17135006ef4a0b624958f4b16cb4ed44312b29f..245757480e146c6556c02931a6c96070d34340ab 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: airGR
 Type: Package
 Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.2.16.12
+Version: 1.2.16.13
 Date: 2019-06-06
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
diff --git a/NEWS.rmd b/NEWS.rmd
index 8e30c857aff76dc949eb64cd1279cfb2855d7063..95705802315712c42f9307a3211cf367b2af564b 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
 
 
 
-### 1.2.16.12 Release Notes (2019-06-06)
+### 1.2.16.13 Release Notes (2019-06-06)
 
 
 #### New features
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index d36c3d6f3c0372444efb284eaaa6a2cd14ede152..e6178e6f4df59d09410f2366cc63b00a90b38dcd 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -605,7 +605,7 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
     
     plot(Xaxis, errorQ[IndPeriod_Plot],
          type = "l", xaxt = "n", yaxt = "n", ylim = ylim1,
-         col = "grey50", lwd = lwd * lwdk,
+         col = "grey40", lwd = lwd * lwdk,
          xlab = "", ylab = "",
          panel.first = abline(h = 0, col = "royalblue"), ...)
     axis(side = 2, at = pretty(ylim1), labels = pretty(ylim1), cex.axis = cex.axis, ...)
@@ -875,7 +875,7 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
       ylim <- log(range(c(Qobs[IndPeriod_Plot][SelectQobsNotZero & SelectQsimNotZero], OutputsModel$Qsim[IndPeriod_Plot][SelectQobsNotZero & SelectQsimNotZero]), na.rm = TRUE))
       plot(log(Qobs[IndPeriod_Plot][SelectQobsNotZero & SelectQsimNotZero]),
            log(OutputsModel$Qsim[IndPeriod_Plot][SelectQobsNotZero & SelectQsimNotZero]),
-           type = "p", pch = 1, cex = 0.9, col = "grey50", lwd = lwd,
+           type = "p", pch = 1, cex = 0.9, col = "grey40", lwd = lwd,
            xlim = ylim, ylim = ylim, xaxt = "n", yaxt = "n", xlab = "", ylab = "", ...)
       abline(a = 0, b = 1, col = "royalblue", lwd = lwd)
       axis(side = 1, at = seqDATA1, labels = seqDATA2, cex = cex.leg, cex.axis = cex.axis, ...)