Commit 5da5ccf0 authored by Dorchies David's avatar Dorchies David
Browse files

lint: RunModel_LLR and RunModel_LBLR

Refs HYCAR-Hydro/airgr#152, HYCAR-Hydro/airgr#153
parent 6e4e1a1a
No related merge requests found
Showing with 9 additions and 10 deletions
+9 -10
...@@ -64,21 +64,20 @@ RunModel_LBLR <- function(InputsModel, RunOptions, Param, QcontribDown) { ...@@ -64,21 +64,20 @@ RunModel_LBLR <- function(InputsModel, RunOptions, Param, QcontribDown) {
} }
## Parameters set up ## Parameters set up
## Param[1L] et Param[2L] : m/s ## Param[1L] and Param[2L] : m/s
TParam <- max(InputsModel$LengthHydro) * 1e3 / Param[1L] / RunOptions$FeatFUN_MOD$TimeStep #T en pdt TParam <- max(InputsModel$LengthHydro) * 1e3 / Param[1L] / RunOptions$FeatFUN_MOD$TimeStep # Tau0 en pdt
KParam <- max(InputsModel$LengthHydro) * 1e3 / Param[2L] / RunOptions$FeatFUN_MOD$TimeStep # T en pdt KParam <- max(InputsModel$LengthHydro) * 1e3 / Param[2L] / RunOptions$FeatFUN_MOD$TimeStep # K0 en pdt
PK <- ifelse(InputsModel$LengthHydro == 0, 0, PK <- ifelse(InputsModel$LengthHydro == 0, 0,
sqrt(InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * KParam) sqrt(InputsModel$LengthHydro / max(InputsModel$LengthHydro)) * KParam)
PT <- ifelse(InputsModel$LengthHydro == 0, 0, PT <- ifelse(InputsModel$LengthHydro == 0, 0,
InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * (TParam + KParam) - PK InputsModel$LengthHydro / max(InputsModel$LengthHydro)) * (TParam + KParam) - PK
PT <- ifelse(PT < 0, 0, PT) PT <- ifelse(PT < 0, 0, PT)
PK[which(PT == 0)] <- ifelse(InputsModel$LengthHydro == 0, 0, PK[which(PT == 0)] <- ifelse(InputsModel$LengthHydro == 0, 0,
InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * (TParam + KParam) InputsModel$LengthHydro/max(InputsModel$LengthHydro)) * (TParam + KParam)
rprop <- PT - floor(PT)
## Lag set up ## Lag set up
lag <- function(t, Qroute, Qupstream){ lag <- function(t, Qroute, Qupstream){
Q <- c(Qupstream[1 : floor(t)], Qroute[1 : (length(Qroute) - floor(t))]) Q <- c(Qupstream[1 : floor(t)], Qroute[1 : (length(Qroute) - floor(t))])
......
...@@ -67,10 +67,10 @@ RunModel_LLR <- function(InputsModel, RunOptions, Param, QcontribDown) { ...@@ -67,10 +67,10 @@ RunModel_LLR <- function(InputsModel, RunOptions, Param, QcontribDown) {
## Parameters set up ## Parameters set up
TParam <- Param[1L] CParam <- Param[1L]
KParam <- Param[2L] KParam <- Param[2L]
PT <- InputsModel$LengthHydro * 1e3 / Param[1L] / RunOptions$FeatFUN_MOD$TimeStep PT <- InputsModel$LengthHydro * 1e3 / CParam / RunOptions$FeatFUN_MOD$TimeStep
PK <- KParam * PT PK <- KParam * PT
C0 <- ifelse(PK == 0, 0, exp(-1/PK)) C0 <- ifelse(PK == 0, 0, exp(-1/PK))
......
Supports Markdown
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