v1.2.16.11 UPDATE: the Error plot of plot.OutputsModel fun now displays Qsim -...

v1.2.16.11 UPDATE: the Error plot of plot.OutputsModel fun now displays Qsim - Qobs instaed of Qsim/Qobs

DESCRIPTION

 Package: airGR
 Type: Package
 Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.2.16.10
+Version: 1.2.16.11
 Date: 2019-06-06
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),

NEWS.rmd

@@ -13,7 +13,8 @@ output:
 
 
 
-6-065-23)
+
+### 1.2.16.11 Release Notes (2019-06-06)
 
 
 #### New features

R/plot.OutputsModel.R

@@ -594,19 +594,22 @@ plot.OutputsModel <- function(x, Qobs = NULL, IndPeriod_Plot = NULL, BasinArea =
   if (BOOLPLOT_Error) {
     kPlot <- kPlot + 1
     mar <- c(3, 5, 1, 5)
-    plotunitError <- "[-]"
     
-    errorQ <- OutputsModel$Qsim / Qobs
+    if (log_scale) {
+      errorQ <- log(OutputsModel$Qsim) - log(Qobs)
+    } else {
+      errorQ <- OutputsModel$Qsim - Qobs
+    }
     par(new = FALSE, mar = mar)
     ylim1 <- range(errorQ[IndPeriod_Plot], na.rm = TRUE)
     
     plot(Xaxis, errorQ[IndPeriod_Plot],
          type = "l", xaxt = "n", yaxt = "n", ylim = ylim1,
          col = "grey50", lwd = lwd * lwdk,
-         xlab = "", ylab = "", log = ifelse(log_scale, "y", ""),
-         panel.first = abline(h = 1, col = "grey", lty = 2), ...)
+         xlab = "", ylab = "",
+         panel.first = abline(h = 0, col = "grey"), ...)
     axis(side = 2, at = pretty(ylim1), labels = pretty(ylim1), cex.axis = cex.axis, ...)
-    mtext(side = 2, paste("flow err.", plotunitError), cex = cex.lab, line = line)
+    mtext(side = 2, paste("flow error", plotunit), cex = cex.lab, line = line)
     if (BOOL_Dates) {
       axis(side = 1, at = Seq1, labels = FALSE, cex.axis = cex.axis, ...)
       axis(side = 1, at = Seq2, labels = Labels2, lwd.ticks = 1.5, cex.axis = cex.axis, ...)