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

feat(saveFlowDB): add hydrological indicators

parent 3832fc4f
......@@ -13,16 +13,19 @@ saveFlowDB <- function(path, InputsModel, OutputsModel, cfg = loadConfig()) {
dir.create(path, showWarnings = FALSE, recursive = TRUE)
# Gauging station list
BS_reseau <- attr(InputsModel, "GRiwrm")
BS_reseau <- BS_reseau[, c("id", "lambert2.x", "lambert2.y", "area", "nom", "down", "length", "model")]
names(BS_reseau) <- c("CdSiteHydro", "lambert2.x", "lambert2.y", "area", "nom", "id_aval", "distance_aval", "model")
readr::write_tsv(BS_reseau, file.path(path, "stations.tsv"))
# Simulated flows
Qsim <- attr(OutputsModel, "Qm3s")
Qsim$DatesR <- format(Qsim$DatesR, format = "%Y-%m-%d")
readr::write_tsv(Qsim, file.path(path, "ts_Qsim.tsv")
)
# Simulated contributing flows (GR outputs)
listQGR <- lapply(names(OutputsModel), function(id) {
if (is.null(OutputsModel[[id]]$QsimDown)) {
return(OutputsModel[[id]]$Qsim_m3)
......@@ -37,8 +40,30 @@ saveFlowDB <- function(path, InputsModel, OutputsModel, cfg = loadConfig()) {
readr::write_tsv(Qcontrib, file.path(path, "ts_Qcontrib.tsv"))
# Global hydrological indicators
indicators <- c("QA",
"VCN10_2", "VCN30_2", "QMNA2",
"VCN10_5", "VCN30_5", "QMNA5",
"VCN10_10", "VCN30_10", "QMNA10",
"QJXA2", "QJXA10", "QJXA20")
lInd <- lapply(
indicators,
calcQcaract,
Q = attr(OutputsModel, "Qm3s")
)
mInd <- do.call(cbind, lInd)
colnames(mInd) <- indicators
dfInd <- cbind(data.frame(Id = rownames(mInd)), mInd)
readr::write_tsv(dfInd, file.path(path, "Q_indicators.tsv"))
# README
file.copy(system.file("LISEZMOI.md", package = "seinebasin2"),
file.path(path, "LISEZMOI.md"),
overwrite = TRUE)
}
}
calcQcaract <- function(indicator, Q) {
FUN <- match.fun(paste0("calc", indicator))
FUN(Q)
}
......@@ -9,10 +9,11 @@ Il contient les fichiers suivants:
- stations.tsv
- ts_Qsim.tsv
- ts_Qcontrib.tsv
- Q_indicators.tsv
Ce sont des fichiers texte avec séparateur tabulation (Tab-Separated-Values).
## Description du fichier `stations.tsv`
## Description du fichier "stations.tsv"
Ce fichier contient la liste des stations hydrométriques utilisées dans le modèle 'seinebasin2' avec les colonnes suivantes:
......@@ -26,7 +27,7 @@ Ce fichier contient la liste des stations hydrométriques utilisées dans le mod
Les sites issus de la banque hydro ont un identifiant commençant par la lettre `H`. Les sites supplémentaires (`A_PMARNE`, `TRANN_01`, `CHAUM_07`) proviennent de stations gérées par l'EPTB Seine Grands Lacs.
## Description du fichier `ts_Qsim.tsv`
## Description du fichier "ts_Qsim.tsv"
Ce fichier contient les séries chronologiques des débits moyens journaliers modélisés sur chacune des stations hydrométriques.
......@@ -34,7 +35,7 @@ La première colonne `DatesR` contient la date de la données moyenne journaliè
Les colonnes suivantes portent le nom de la station hydrométique correspondant à la colonne `CdSiteHydro` du fichier `stations.tsv` et elles contiennent le débit moyen journalier modélisé en m<sup>3</sup>/s.
## Description du fichier `ts_Qcontrib.tsv`
## Description du fichier "ts_Qcontrib.tsv"
Ce fichier contient les séries chronologiques de la contribution moyenne journalière modélisée sur chaque bassin versant intermédiaire. Les bassins versants intermédiaires sont délimités par le bassin versant dont l'exutoire est la station hydrométrique donnant son nom au bassin versant auquels on soustrait tous les bassins versants identifiés à l'amont.
......@@ -42,3 +43,17 @@ La première colonne `DatesR` contient la date de la données moyenne journaliè
Les colonnes suivantes portent le nom de la station hydrométique correspondant à la colonne `CdSiteHydro` du fichier `stations.tsv` et elles contiennent la contribution moyenne journalière du bassin versant intermédiaire modélisée en m<sup>3</sup>/s.
## Description du fichier "Q_indicators.tsv"
Ce fichier contient les indicateurs hydrologiques suivants en m³/s:
- `QA`: débit annuel moyen
- `VCN10_2`, `VCN30_2`, `VCN10_5`, `VCN_30_5`, `VCN10_10`, `VCN30_10`: avec `VCN[k]_[Y]` le débit moyen minimum sur une durée de `k` jours de période de retour `Y` ans
- `QMNA2`, `QMNA5`, `QMNA10`: avec `QMNA[Y]` le débit moyen minimum mensuel de période de retour `Y` ans
- `QJXA2`, `QJXA10`, `QJXA20`: avec `QJXA[Y]` le débit moyen journalier maximum annuel de période de retour `Y` ans
Pour chaque indicateur, le fichier contient:
- l'indicateur calculé sur les débits observés (colonnes `[indicateur].obs`)
- l'indicateur calculé sur les débits simulés sur la période de débits observés (colonnes `[indicateur].sim`)
- le rapport entre l'indicateur calculé sur les débits simulés et celui sur les débits observés (colonnes `R-[indicateur]`)
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