diff --git a/DESCRIPTION b/DESCRIPTION
index 4974dc87ede8dee97fdefacc6edbdc1125b6298a..c9ea9e23c990fa98cd8aee0d197a99851fc5cd57 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.5
+Version: 1.2.16.6
 Date: 2019-05-24
 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 02290d96b3c7f7948286db0629dd495b69db8ca9..4390a999f7dca8501a20275dcdbe145ff645866c 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
 
 
 
-### 1.2.16.5 Release Notes (2019-05-23)
+### 1.2.16.6 Release Notes (2019-05-23)
 
 
 #### New features
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index 7493ec7b760d748f1edc8d8700c14db29012b98d..08d46ffad2dfb9d63ea3b5a342161ee034b75a1f 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -1,6 +1,7 @@
 plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea = NULL, which = "synth", log_scale = FALSE,
-                              cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1, verbose = TRUE,
-                              LayoutMat = NULL, LayoutWidths = rep.int(1, ncol(LayoutMat)), LayoutHeights = rep.int(1, nrow(LayoutMat)), ...) {
+                              cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1,
+                              LayoutMat = NULL, LayoutWidths = rep.int(1, ncol(LayoutMat)), LayoutHeights = rep.int(1, nrow(LayoutMat)),
+                              verbose = TRUE, ...) {
   
   
   ## save default graphical parameters and resetting on exit
diff --git a/man/plot.OutputsModel.Rd b/man/plot.OutputsModel.Rd
index f8abac14ed38d2b3c09935097acc1fc7490feee2..e3e3fa85246c4a3cf4e769def8ec56ecf1686943 100644
--- a/man/plot.OutputsModel.Rd
+++ b/man/plot.OutputsModel.Rd
@@ -12,10 +12,11 @@
 \usage{
 \method{plot}{OutputsModel}(x, Qobs = NULL, IndPeriod_Plot = NULL,
   BasinArea = NULL, which = "synth", log_scale = FALSE,
-  cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1, verbose = TRUE,
+  cex.axis = 1, cex.lab = 0.9, cex.leg = 0.9, lwd = 1,
   LayoutMat = NULL,
   LayoutWidths = rep.int(1, ncol(LayoutMat)),
-  LayoutHeights = rep.int(1, nrow(LayoutMat)), ...)
+  LayoutHeights = rep.int(1, nrow(LayoutMat)),
+  verbose = TRUE, ...)
 }
 
 
@@ -40,6 +41,12 @@
 
 \item{lwd}{(optional) [numeric] the line width (a positive number)}
 
+\item{LayoutMat}{(optional) [numeric] a matrix object specifying the location of the next N figures on the output device. Each value in the matrix must be 0 or a positive integer. If N is the largest positive integer in the matrix, then the integers {1, \dots, N-1} must also appear at least once in the matrix (see \code{\link{layout}})}
+
+\item{LayoutWidths}{(optional) [numeric] a vector of values for the widths of columns on the device (see \code{\link{layout}})}
+
+\item{LayoutHeights}{(optional) [numeric] a vector of values for the heights of rows on the device (see \code{\link{layout}})}
+
 \item{verbose}{(optional) [boolean] indicating if the function is run in verbose mode or not, default = \code{TRUE}}
 
 \item{...}{other parameters to be passed through to plotting functions}