feat(06-evolution): add plot of monthly flows

Refs #12

bookdown/06-evolution_indicateurs_hydrologiques.Rmd

@@ -204,6 +204,47 @@ Qhist_month5 <- lapply(historiQ, loadHistoriQ, indicator = "Q_monthly_5years")
 Qhist_indicators <- lapply(historiQ, loadHistoriQ, indicator = "Q_indicators")
 ```
 
+```{r}
+plotDataMonth <- function(rcp, period, station, drias, hist) {
+    scenarios <- names(drias)
+    names(scenarios) <- scenarios
+    drias <- lapply(scenarios, function(x) {
+        if (grepl(rcp, x, fixed = TRUE)) {
+            
+            drias[[x]][[period]][station, ]
+        } else {
+            NULL
+        }
+    })
+    drias <- drias[!sapply(drias,is.null)]
+    df <- as.data.frame(do.call(cbind, drias))
+    df$Mois <- seq.int(12)
+    dfHist <- as.data.frame(do.call(cbind, lapply(hist, function(x) x[station, ])))
+    dfHist$Mois <- seq.int(12)
+    library(ggplot2)
+    ggplot(tidyr::gather(df, key = "scenario", value = "Q", -Mois)) +
+        geom_line(aes(x = Mois, y = Q, color = scenario), size = 1) +
+        geom_line(data = tidyr::gather(dfHist, key = "historique", value = "Q", -Mois), 
+                  aes(x = Mois, y = Q, linetype = historique), 
+                  size = 1) +
+        scale_x_continuous(breaks = seq.int(12), 
+                           labels =  strsplit("JFMAMJJASOND", "")[[1]])
+        
+}
+```
+
+Débits moyens mensuels à Paris sur la période de référence:
+
+```{r}
+plotDataMonth("rcp4.5", "ref", "H5920010", Q_month, Qhist_month)
+```
+
+Débits moyens à Paris sur la période 2071-2100 avec le scénario d'émission RCP4.5:
+
+```{r}
+plotDataMonth("rcp4.5", "end", "H5920010", Q_month, Qhist_month)
+```
+
 
 Les tableaux de synthèse fournissent la valeur minimale, médiane et maximale des évolutions parmi les couples scénario/modèle climatiques.
 
@@ -256,7 +297,7 @@ Exemple pour la température moyenne mensuelle à Paris:
 
 ```{r}
 deltaTM <- calcDelta(T_month, rcp = "rcp4.5", period = "end", delta = "+")
-knitr::kable(tableDeltaStation("H8012010", deltaTM), digits = 1)
+knitr::kable(tableDeltaStation("H5920010", deltaTM), digits = 1)
 ```
 
 **N.B.:** les données climatiques correspondent ici au données moyennes du sous-bassin versant et pas le bassin versant entier. Il faudrait agréger les données des bassins amont pour avoir une moyenne du bassin pour chaque station.
@@ -265,14 +306,14 @@ Exemple sur le débit mensuel moyen à Paris
 
 ```{r}
 deltaQM <-  calcDelta(Q_month, rcp = "rcp4.5", period = "end")
-knitr::kable(tableDeltaStation("H8012010", deltaQM), digits = 1)
+knitr::kable(tableDeltaStation("H5920010", deltaQM, Qhist_month), digits = 1)
 ```
 
 Exemple sur tous les indicateurs hydrologiques à Paris:
 
 ```{r}
 deltaIndicators <- calcDelta(Q_indicators, rcp = "rcp4.5", period = "end")
-knitr::kable(tableDeltaStation("H8012010", deltaIndicators), digits = 1)
+knitr::kable(tableDeltaStation("H5920010", deltaIndicators, Qhist_indicators), digits = 1)
 ```
 ## Cartes d'évolutions