From 898a459a7f2f68c764884f150042e8395e1ee64d Mon Sep 17 00:00:00 2001
From: Dorchies David <david.dorchies@inrae.fr>
Date: Tue, 2 Aug 2022 11:38:38 +0200
Subject: [PATCH] refactor(RunModel_LRQ): express K0 as a reservoir capacity in
mm/time step
Refs HYCAR-Hydro/airgr#152, HYCAR-Hydro/airgr#153
---
R/RunModel_LRQ.R | 7 ++-----
1 file changed, 2 insertions(+), 5 deletions(-)
diff --git a/R/RunModel_LRQ.R b/R/RunModel_LRQ.R
index 0d512d6..46bcd68 100644
--- a/R/RunModel_LRQ.R
+++ b/R/RunModel_LRQ.R
@@ -66,11 +66,10 @@ RunModel_LRQ <- function(InputsModel, RunOptions, Param, QcontribDown) {
## Parameters set up
TParam <- Param[1L]
- KParam <- Param[2L]
+ KParam <- Param[2L] # KParam en mm/pdt
PT <- InputsModel$LengthHydro * 1e3 / TParam / RunOptions$FeatFUN_MOD$TimeStep
- PK <- KParam * PT
-
+ PK <- sqrt(KParam * InputsModel$BasinAreas[-length(InputsModel$BasinAreas)] * 1e-3) * PT #m3/2 pdt1/2
## Lag set up
@@ -120,8 +119,6 @@ RunModel_LRQ <- function(InputsModel, RunOptions, Param, QcontribDown) {
QsimDown_input <- matrix(QsimDown *
InputsModel$BasinAreas[length(InputsModel$BasinAreas)] * 1e3, ncol = 1)
-
-
for (upstream_basin in seq_len(NbUpBasins)) {
Qupstream_m3 <- c(IniStates[[upstream_basin]],
InputsModel$Qupstream[IndPeriod1, upstream_basin])
--
GitLab