diff --git a/vignettes/01_First_network.Rmd b/vignettes/01_First_network.Rmd
index 2b42ccc993ee18b0305d55b7af51769379b68a3c..1bb92f6457630c6b0b4f2d405206e6e8dc1e5904 100644
--- a/vignettes/01_First_network.Rmd
+++ b/vignettes/01_First_network.Rmd
@@ -71,6 +71,9 @@ l <- lapply(urls, load_ts)
 gits <- Gits(ginet$id[1], l[[ginet$id[1]]], cols = list(date = "Date", Precip = "Ptot", PotEvap = "ETP", Qobs = "Qnat"))
 
 for(id in ginet$id) {
+  l[[id]]$Qnat
+  l[[id]]$Qnat <- l[[id]]$Qnat * 86.4 / girop$area[girop$id == id]
+  l[[id]]$Qnat[l[[id]]$Qnat < 0] <- NA
   gits <- merge(gits, Gits(id, l[[id]], cols = list(date = "Date", Precip = "Ptot", PotEvap = "ETP", Qobs = "Qnat")))
 }
 ```
diff --git a/vignettes/02_First_run.Rmd b/vignettes/02_First_run.Rmd
index 76176318f3142d3777a1eb0002dbf50e81e4bfc2..52a4ebd8674af31c75e9320ee38f2068c2e4e79a 100644
--- a/vignettes/02_First_run.Rmd
+++ b/vignettes/02_First_run.Rmd
@@ -55,7 +55,12 @@ OutputsModels <- RunModelGriwrm(
 ## Plot the result for each basin
 
 ```{r, fig.height = 5, fig.width = 8}
-htmltools::tagList(lapply(names(OutputsModels), function(x) plot(OutputsModels[[x]], main = x)))
+htmltools::tagList(lapply(
+  names(OutputsModels), 
+  function(x) {
+    plot(OutputsModels[[x]], Qobs = gits[[x]]$Qobs[IndPeriod_Run] , main = x)
+  }
+))
 ```