Commit fb5320c2 authored by Dorchies David's avatar Dorchies David
Browse files

feat: add RCP8.5 scenarios

Refs #15
parent 0194d181
Pipeline #31848 passed with stages
in 52 minutes and 16 seconds
# Calage du modèle {#calage}
```{r setup, echo=FALSE, results='hide'}
```{r setup, include=FALSE}
library(seinebasin2)
cfg <- loadConfig()
```
......
# Débits désinfluencés {#desinfluencement}
```{r setup, echo=FALSE, results='hide'}
```{r setup, include=FALSE}
library(seinebasin2)
cfg <- loadConfig()
```
......
# Les scénarios climatiques DRIAS 2020 {#drias2020}
```{r setup, echo=FALSE, results='hide'}
```{r setup, include=FALSE}
library(seinebasin2)
library(sp)
library(dplyr)
......@@ -73,22 +73,29 @@ drias_data_files <- listDataFiles(file.path(cfg$hydroclim$path, "drias"), cfg =
scenarioDriasFiles <- sapply(cfg$hydroclim$drias$scenarios, getDrias2020filenames, drias_data_files)
```
La liste des scénarios sélectionnés pour l'étude est la suivante : `r paste(cfg$hydroclim$drias$scenarios, collapse = ", ")`. Pour chaque scenario, les données sont présentes pour une période de référence (1950-2005) et une projection RCP4.5 (2006-2100).
La liste des scénarios sélectionnés pour l'étude est la suivante : `r paste(cfg$hydroclim$drias$scenarios, collapse = ", ")`. Pour chaque scenario, les données sont présentes pour une période de référence (1950-2005) et deux projections (2006-2100) pour les scénarios d'émission RCP4.5 et RCP 8.5.
### Aggrégation des données DRIAS 2020
```{r, eval=!cfg$data$cloud}
driasPath <- getDataPath(cfg$hydroclim$path)
lapply(colnames(scenarioDriasFiles), function(scenario) {
BasinsObs <- createBasinsObs(scenarioDriasFiles[, scenario],
saveBasinsObs <- function(rcp, scenario) {
files <- scenarioDriasFiles[, scenario]
files <- files[grep(paste0("Historical|", rcp), files)]
BasinsObs <- createBasinsObs(files,
dfMailles = dfMailles, cfg = cfg)
path <- file.path(
driasPath,
paste0(
paste("BasinObs", cfg$hydroclim$drias$rcp[2], gsub("/", "_", scenario), sep="_"),
paste("BasinObs", rcp, gsub("/", "_", scenario), sep="_"),
".RDS"
)
)
saveRDS(BasinsObs, path)
})
}
mapply(saveBasinsObs,
rcp = cfg$hydroclim$drias$rcp[-1],
scenario = rep(colnames(scenarioDriasFiles), 2))
```
......@@ -4,7 +4,7 @@ Dans ce chapitre, nous simulerons les débits naturels (c.-à-d. sans prise en c
## Simulation et enregistrement des chroniques de débits naturels
```{r setup, echo=FALSE, results='hide'}
```{r setup, include=FALSE}
library(seinebasin2)
cfg <- loadConfig()
```
......@@ -32,7 +32,7 @@ Il suffit ensuite d'effectuer l'opération pour tous les scénarios&nbsp;:
```{r}
rcps <- cfg$hydroclim$drias$rcp[-1]
scenarios <- gsub("/", "_", cfg$hydroclim$drias$scenarios)
scenarios <- rep(scenarios, length(rcps))
scenarios <- rep(scenarios, each = length(rcps))
```
```{r, eval=!cfg$data$cloud, message=FALSE}
......@@ -118,7 +118,14 @@ plotVCN <- function(CdSiteHydro, df) {
xlab("Année") + ylab(bquote('Débit naturel'~(m^3/s))) +
labs(color = "Source", shape = "Source", size = "Source")
}
plotVCN("H5920010", dfVCN10)
```
```{r, fig.cap="VCN10 annuel et moyenne mobile sur 5 ans des débits naturels historiques des bases Hydratec et Inrae et des projections Drias 2020 historiques et futures du scénario d'émission RCP 4.5 du modèle Inrae sur la Seine à Paris"}
plotVCN("H5920010", dfVCN10[grep("Hydratec|Inrae|rcp4\\.5", dfVCN10$source), ])
```
```{r, fig.cap="VCN10 annuel et moyenne mobile sur 5 ans des débits naturels historiques des bases Hydratec et Inrae et des projections Drias 2020 historiques et futures du scénario d'émission RCP 4.5 du modèle Inrae sur la Seine à Paris"}
plotVCN("H5920010", dfVCN10[grep("Hydratec|Inrae|rcp8\\.5", dfVCN10$source), ])
```
Enregistrement des graphiques pour toutes les stations sur le cloud&nbsp;:
......
# Evolution des indicateurs hydrologiques
```{r setup, echo=FALSE, results='hide'}
```{r setup, include=FALSE}
library(seinebasin2)
cfg <- loadConfig()
```
......
......@@ -16,8 +16,12 @@ for (scenario in cfg$hydroclim$drias$scenarios) {
url_download <- file.path(url, fileName)
message("Download ", url_download, "...", appendLF = FALSE)
path_write <- getDataPath(cfg$hydroclim$path, "drias", fileName, cfg = cfg)
httr::GET(url_download, httr::write_disk(path_write, overwrite = TRUE), httr::progress())
message(" Done")
if (!file.exists(path_write)) {
httr::GET(url_download, httr::write_disk(path_write, overwrite = TRUE), httr::progress())
message(" Done")
} else {
message(" Already done")
}
}
}
}
......@@ -19,6 +19,7 @@ default:
rcp:
- historical
- rcp4.5
- rcp8.5
bbox:
x:
min: 532000
......
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