From 057d235cba6709dede7be87ef5af1570bf131c70 Mon Sep 17 00:00:00 2001
From: Dorchies David <david.dorchies@inrae.fr>
Date: Tue, 24 May 2022 21:02:18 +0200
Subject: [PATCH] fix(RunModel_LLR): correctly handle LengthHydro == 0
Refs HYCAR-Hydro/airgr#152, HYCAR-Hydro/airgr#153
---
R/RunModel_LLR.R | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/R/RunModel_LLR.R b/R/RunModel_LLR.R
index 63bd8b3..5ab9884 100644
--- a/R/RunModel_LLR.R
+++ b/R/RunModel_LLR.R
@@ -70,9 +70,11 @@ RunModel_LLR <- function(InputsModel, RunOptions, Param, QcontribDown) {
## propagation time from upstream meshes to outlet
PT <- floor(InputsModel$LengthHydro * 1e3 / TParam / RunOptions$FeatFUN_MOD$TimeStep)
- PK <- sqrt(InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * KParam
+ PK <- ifelse(InputsModel$LengthHydro == 0,
+ 0,
+ sqrt(InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * KParam)
- C0 <- exp(-1/PK)
+ C0 <- ifelse(PK ==0, 0, exp(-1/PK))
C1 <- ((PK/1) * (1-C0))-C0
C2 <- 1 - (PK/1) * (1-C0)
## Lag model computation
--
GitLab