Commit 2399efff authored by Dorchies David's avatar Dorchies David
Browse files

change: update scale of claibration maps

parent b3c92b47
Pipeline #33621 passed with stages
in 58 minutes and 3 seconds
......@@ -112,16 +112,10 @@ calcRatio <- function(Qsim, Qobs, FUN) {
pSim / pObs
}
plotRatio <- function(Qsim, Qobs, FUN, title) {
plotRatio <- function(Qsim, Qobs, FUN, title,
breaks = c(0, 0.25, 0.50, 0.75, 0.90,
1.10, 1.25, 1.50, 2, 20)) {
r <- calcRatio(Qsim, Qobs, FUN)
# Define scale centered around 1
decVal <- r[r < 1]
incVal <- r[!(r < 1)]
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)
plot_seine_map(r, breaks, title)
}
```
......@@ -131,7 +125,9 @@ plotRatio <- function(Qsim, Qobs, FUN, title) {
### Ratio débit moyen simulé et observé
```{r, fig.cap="Carte des R-QA des Bassins versants intermédiaires"}
plotRatio(Qsim, Qobs, calcQA, "R-QA")
plotRatio(Qsim, Qobs, calcQA, "R-QA",
breaks = c(0, 0.6, 0.8, 0.90, 0.95,
1.05, 1.10, 1.20, 1.40, 1.6))
```
### VCN30 de période de retour 2 ans
......@@ -195,14 +191,14 @@ indicators_periods <- list(
)
InputsModel <- CreateInputsModel(BasinsObs_with)
saveFlowDB(path = getDataPath(cfg$calibration$path, "Qobs", cfg = cfg),
InputsModel = InputsModel,
saveFlowDB(path = getDataPath(cfg$calibration$path, "Qobs", cfg = cfg),
InputsModel = InputsModel,
OutputsModel = Qobs_m3s,
indicators_periods = indicators_periods,
cfg = cfg)
saveFlowDB(path = getDataPath(cfg$calibration$path, "Qsim", cfg = cfg),
InputsModel = InputsModel,
OutputsModel = OM_with,
saveFlowDB(path = getDataPath(cfg$calibration$path, "Qsim", cfg = cfg),
InputsModel = InputsModel,
OutputsModel = OM_with,
indicators_periods = indicators_periods,
cfg = cfg)
```
......
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