Commit 0afd30a2 authored by Dorchies David's avatar Dorchies David
Browse files

feat(indicateurs): change RCP4.5 to RCP8.5

Refs #12
parent 2399efff
......@@ -168,19 +168,19 @@ Q_indicators <- loadAllIndicators("Q_indicators")
### Données climatiques moyennes mensuelles
```{r, fig.cap="Précipitations moyennes mensuelles du bassin de la Seine à Paris sur la période 1976-2005"}
plot_monthly_mean("rcp4.5", "ref", "H5920010", P_month, list(obs = Pobs_month), "Precipitation (mm)")
plot_monthly_mean("rcp8.5", "ref", "H5920010", P_month, list(obs = Pobs_month), "Precipitation (mm)")
```
```{r, fig.cap="Précipitations moyennes mensuelles du bassin de la Seine à Paris sur la période 2071-2100"}
plot_monthly_mean("rcp4.5", "end", "H5920010", P_month, list(obs = Pobs_month), "Precipitation (mm)")
plot_monthly_mean("rcp8.5", "end", "H5920010", P_month, list(obs = Pobs_month), "Precipitation (mm)")
```
```{r, fig.cap="Températures moyennes mensuelles du bassin de la Seine à Paris sur la période 1976-2005"}
plot_monthly_mean("rcp4.5", "ref", "H5920010", T_month, list(obs = Tobs_month), "Temperature (°C)")
plot_monthly_mean("rcp8.5", "ref", "H5920010", T_month, list(obs = Tobs_month), "Temperature (°C)")
```
```{r, fig.cap="Températures moyennes mensuelles du bassin de la Seine à Paris sur la période 2071-2100"}
plot_monthly_mean("rcp4.5", "end", "H5920010", T_month, list(obs = Tobs_month), "Temperature (°C)")
plot_monthly_mean("rcp8.5", "end", "H5920010", T_month, list(obs = Tobs_month), "Temperature (°C)")
```
......@@ -189,13 +189,13 @@ plot_monthly_mean("rcp4.5", "end", "H5920010", T_month, list(obs = Tobs_month),
Débits moyens mensuels à Paris sur la période de référence:
```{r, fig.asp = 0.7}
plot_monthly_mean("rcp4.5", "ref", "H5920010", Q_month, Qhist_month)
plot_monthly_mean("rcp8.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:
Débits moyens à Paris sur la période 2071-2100 avec le scénario d'émission RCP8.5:
```{r, fig.asp = 0.7}
plot_monthly_mean("rcp4.5", "end", "H5920010", Q_month, Qhist_month)
plot_monthly_mean("rcp8.5", "end", "H5920010", Q_month, Qhist_month)
```
......@@ -249,14 +249,14 @@ tableDeltaStation <- function(station, delta, ind_hist = NULL) {
Exemple pour les pluies moyennes mensuelles pour les 6 premières stations:
```{r}
deltaPM <- calcDelta(P_month, rcp = "rcp4.5", period = "end")
deltaPM <- calcDelta(P_month, rcp = "rcp8.5", period = "end")
lapply(deltaPM, function(x) t(head(x)))
```
Exemple pour la température moyenne mensuelle à Paris:
```{r}
deltaTM <- calcDelta(T_month, rcp = "rcp4.5", period = "end", delta = "+")
deltaTM <- calcDelta(T_month, rcp = "rcp8.5", period = "end", delta = "+")
knitr::kable(tableDeltaStation("H5920010", deltaTM), digits = 1)
```
......@@ -265,27 +265,22 @@ knitr::kable(tableDeltaStation("H5920010", deltaTM), digits = 1)
Exemple sur le débit mensuel moyen à Paris
```{r}
deltaQM <- calcDelta(Q_month, rcp = "rcp4.5", period = "end")
deltaQM <- calcDelta(Q_month, rcp = "rcp8.5", period = "end")
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")
deltaIndicators <- calcDelta(Q_indicators, rcp = "rcp8.5", period = "end")
knitr::kable(tableDeltaStation("H5920010", deltaIndicators, Qhist_indicators), digits = 1)
```
## Cartes d'évolutions
```{r, fig.asp = 1}
plot_map_delta <- function(r, title) {
# Define scale centered around 0
decVal <- r[r < 0]
incVal <- r[!(r < 0)]
decQuant <- quantile(decVal, probs = c(0, 0.1, 0.4, 0.6, 0.9))
incQuant <- quantile(incVal, probs = c(0.1, 0.4, 0.6, 0.9, 1))
breaks <- c(decQuant, incQuant)
breaks[is.na(breaks)] <- 0
plot_map_delta <- function(r, title,
breaks = c(-60, -30, -20, -10, -5,
5, 10, 20, 30, 1.6, 100)) {
plot_seine_map(r, breaks, title)
}
```
......@@ -295,9 +290,9 @@ plot_map_delta(deltaIndicators$med[, "QA"], "QA - évolution RCP4.5 scénario m
```
```{r, fig.asp = 1}
plot_map_delta(deltaIndicators$min[, "QMNA5"], "QMNA5 - évolution RCP4.5 scénario minimum (%)")
plot_map_delta(deltaIndicators$med[, "QMNA5"], "QMNA5 - évolution RCP4.5 scénario minimum (%)")
```
```{r, fig.asp = 1}
plot_map_delta(deltaIndicators$max[, "QJXA10"], "QJXA10 - évolution RCP4.5 scénario maximum (%)")
plot_map_delta(deltaIndicators$med[, "QJXA10"], "QJXA10 - évolution RCP4.5 scénario maximum (%)")
```
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment