From b2a0ca93ec8a4aed8a0f8504b555b3cb02e3fd4c Mon Sep 17 00:00:00 2001
From: Dorchies David <david.dorchies@irstea.fr>
Date: Sat, 23 May 2020 15:38:38 +0200
Subject: [PATCH] feat: Add observed flows to vignette 02_First_run.Rmd
Refs #2
---
vignettes/01_First_network.Rmd | 3 +++
vignettes/02_First_run.Rmd | 7 ++++++-
2 files changed, 9 insertions(+), 1 deletion(-)
diff --git a/vignettes/01_First_network.Rmd b/vignettes/01_First_network.Rmd
index 2b42ccc..1bb92f6 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 7617631..52a4ebd 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)
+ }
+))
```
--
GitLab