diff --git a/DESCRIPTION b/DESCRIPTION
index 11cdda52dff05c0e8848016c825ddda5291bd757..576b8e983a54e3e20d9d71ff071b22b5b6ae2ecc 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.9.15
+Version: 1.2.9.16
 Date: 2019-03-18
 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 214188665d81eca27ba499885a21a2e245a63278..e0c999a897f7ec349d7f93b3e4cfd06802dc3bf7 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -13,7 +13,7 @@ output:
 
 
 
-### 1.2.9.15 Release Notes (2019-03-18) 
+### 1.2.9.16 Release Notes (2019-03-18) 
 
 
 
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index 647583083757ddfd988e2a6d325dcef3edc0f966..fddf412e2db1b3328d63150fa3cd0ce7344f8e87 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -108,6 +108,9 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
       NameTS <- "month"
       plotunit <- "[mm/month]"
       formatAxis <- "%m/%Y"
+      if (format(OutputsModel$DatesR[1L], format = "%d") != "01") {
+        OutputsModel$DatesR <- as.POSIXlt(format(OutputsModel$DatesR, format = "%Y-%m-01"), tz = "UTC", format = "%Y-%m-%d")
+      }
     }
     if (inherits(OutputsModel, "yearly") &
         TimeStep %in% (c(365, 366) * 24 * 60 * 60)) {