From 7f8ece8db8ad1e2f8cdaff1cd5146123efc984d3 Mon Sep 17 00:00:00 2001
From: Dorchies David <david.dorchies@inrae.fr>
Date: Tue, 16 Mar 2021 11:02:52 +0100
Subject: [PATCH] fix(RunModel_Lag): bugs on StateEnd
- Add commentaries with debug messages
Fix #103, #104
---
R/RunModel_Lag.R | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/R/RunModel_Lag.R b/R/RunModel_Lag.R
index a832fd9d..b246b1f4 100644
--- a/R/RunModel_Lag.R
+++ b/R/RunModel_Lag.R
@@ -75,6 +75,7 @@ RunModel_Lag <- function(InputsModel, RunOptions, Param) {
rep(0, floor(PT[x] + 1))
})
}
+ #message("Initstates: ",paste(IniStates, collapse = ", "))
for (upstream_basin in seq_len(NbUpBasins)) {
Qupstream <- c(IniStates[[upstream_basin]],
@@ -96,11 +97,14 @@ RunModel_Lag <- function(InputsModel, RunOptions, Param) {
}
# Convert back Qsim to mm
OutputsModel$Qsim <- OutputsModel$Qsim / sum(InputsModel$BasinAreas, na.rm = TRUE) / 1e3
+ #message("Qsim: ",paste(OutputsModel$Qsim, collapse = ", "))
if ("StateEnd" %in% RunOptions$Outputs_Sim) {
OutputsModel$StateEnd$SD <- lapply(seq(NbUpBasins), function(x) {
- Qupstream[(LengthTs - floor(PT[x])):LengthTs]
+ lastTS <- RunOptions$IndPeriod_Run[length(RunOptions$IndPeriod_Run)]
+ InputsModel$Qupstream[(lastTS - floor(PT[x])):lastTS, x]
})
+ #message("StateEnd: ",paste(OutputsModel$StateEnd$SD, collapse = ", "))
}
return(OutputsModel)
--
GitLab