diff --git a/DESCRIPTION b/DESCRIPTION
index 9a32e72be6e53b2c2566e4115925fa83d350ac32..b17135006ef4a0b624958f4b16cb4ed44312b29f 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.11
+Version: 1.2.16.12
 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 c0f9e0b3afe833239dcb2b7bd541ba461d1a6d9d..8e30c857aff76dc949eb64cd1279cfb2855d7063 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
 
 
 
-### 1.2.16.11 Release Notes (2019-06-06)
+### 1.2.16.12 Release Notes (2019-06-06)
 
 
 #### New features
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index e87e934f66cfb671cd53040d2b9c3d050fa9fec2..d36c3d6f3c0372444efb284eaaa6a2cd14ede152 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -607,7 +607,7 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
          type = "l", xaxt = "n", yaxt = "n", ylim = ylim1,
          col = "grey50", lwd = lwd * lwdk,
          xlab = "", ylab = "",
-         panel.first = abline(h = 0, col = "grey"), ...)
+         panel.first = abline(h = 0, col = "royalblue"), ...)
     axis(side = 2, at = pretty(ylim1), labels = pretty(ylim1), cex.axis = cex.axis, ...)
     mtext(side = 2, paste("flow error", plotunit), cex = cex.lab, line = line)
     if (BOOL_Dates) {
@@ -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 = par("fg"), lwd = lwd,
+           type = "p", pch = 1, cex = 0.9, col = "grey50", 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, ...)