Fix: correct area of intermediate basin passed to lag model

Refs #2

Fix: correct area of intermediate basin passed to lag model
Refs #2
3a1003f6 · Dorchies David · f495d16f · 3a1003f6
Commit 3a1003f6 authored 4 years ago by Dorchies David
Hide whitespace changes
Inline Side-by-side

Showing

with 7 additions and 1 deletion
+7 -1
--- a/R/RunModelGriwrm.R
+++ b/R/RunModelGriwrm.R
@@ -79,7 +79,13 @@ SetAirGrInputsAndOptions <- function(id, ginet, girop, gits, OutputsModels, IndP
      }
    }
    LengthHydro <- matrix(ginet$length[girop$id %in% UpstrNodes] , nrow = 1)
-    BasinAreas <- matrix(c(girop$area[girop$id %in% UpstrNodes], girop$area[girop$id == id]) , nrow = 1)
+    BasinAreas <- matrix(
+      c(
+        girop$area[girop$id %in% UpstrNodes],
+        girop$area[girop$id == id] - sum(girop$area[girop$id %in% UpstrNodes])
+      ),
+      nrow = 1
+    )
  }
  # Set model inputs
  InputsModel <- CreateInputsModel(