fix(RunModel_LLR): correctly handle LengthHydro == 0

Refs HYCAR-Hydro/airgr#152, HYCAR-Hydro/airgr#153

fix(RunModel_LLR): correctly handle LengthHydro == 0
Refs HYCAR-Hydro/airgr#152, HYCAR-Hydro/airgr#153
057d235c · Dorchies David · 1b9865fa · 057d235c
Commit 057d235c authored 3 years ago by Dorchies David
Hide whitespace changes
Inline Side-by-side

Showing

with 4 additions and 2 deletions
+4 -2
--- 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