diff --git a/DESCRIPTION b/DESCRIPTION
index 0f9697c7183df34a7e7640c69844ce48c08161b2..625e2de4f78fcd78855ad0766d6261c44d90a60b 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.3.2.49
+Version: 1.3.2.50
 Date: 2019-11-08
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
diff --git a/NEWS.md b/NEWS.md
index 2abfdb8e536bb99294217a08a342f3145c4f9c32..505d232fbe1d44d947e9126119b3fd3737c860f8 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -2,7 +2,7 @@
 
 
 
-### 1.3.2.49 Release Notes (2019-11-08)
+### 1.3.2.50 Release Notes (2019-11-08)
 
 
 #### Bug fixes
diff --git a/R/plot.OutputsModel.R b/R/plot.OutputsModel.R
index deca5dc582c0d3b5491c10bb10021581941a9f3a..9168980fd5b6400ad20ae962a34cecded0410386 100644
--- a/R/plot.OutputsModel.R
+++ b/R/plot.OutputsModel.R
@@ -487,10 +487,11 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
     mar <- c(3, 5, 1, 5)
     par(new = FALSE, mar = mar)
     
-    DATA2 <- Qobs
-    DATA2[!SelectQobsNotZero] <- mean(Qobs, na.rm = TRUE) / 10000
-    DATA2 <- log(DATA2)
-    
+    if (BOOL_Qobs) {    
+      DATA2 <- Qobs
+      DATA2[!SelectQobsNotZero] <- mean(Qobs, na.rm = TRUE) / 10000
+      DATA2 <- log(DATA2)
+    }
     DATA3 <- OutputsModel$Qsim
     DATA3[!SelectQsimNotZero] <- mean(OutputsModel$Qsim, na.rm = TRUE) / 10000
     DATA3 <- log(DATA3)