Commit 9a398d72 authored by David's avatar David
Browse files

fix: User defined donor for Delavenne

- node with user defined donor for a priori parameters works
- rename InputsModel item isUngauged to inUngaugedCluster

Refs #157
2 merge requests!93Draft: Version 0.7.0,!89Resolve "CreateInputsCrit: allow apriori node not only at upstream"
Pipeline #55593 passed with stage
in 8 minutes and 53 seconds
Showing with 10 additions and 6 deletions
+10 -6
......@@ -31,7 +31,7 @@ Calibration.GRiwrmInputsModel <- function(InputsModel,
OutputsModel <- list()
class(OutputsModel) <- append("GRiwrmOutputsModel", class(OutputsModel))
b <- sapply(InputsModel, function(IM) !IM$isUngauged)
b <- sapply(InputsModel, function(IM) !IM$inUngaugedCluster)
gaugedIds <- names(b[b])
for (id in gaugedIds) {
......
......@@ -59,7 +59,7 @@ CreateInputsCrit.GRiwrmInputsModel <- function(InputsModel,
"\nIf possible, set this apriori id as the donor of the node \"",
id,"\" to force the calibration sequence order")
}
if (InputsModel[[AprioriIds[id]]]$isUngauged &
if (InputsModel[[AprioriIds[id]]]$inUngaugedCluster &
InputsModel[[AprioriIds[id]]]$gaugedId == id) {
stop("'AprioriIds': the node \"", AprioriIds[id],
"\" is ungauged, use a gauged node instead")
......
......@@ -331,10 +331,14 @@ CreateOneGRiwrmInputsModel <- function(id, griwrm, DatesR, ..., Qobs, Qmin, Qrel
# Add the model function
InputsModel$FUN_MOD <- FUN_MOD
featModel <- .GetFeatModel(InputsModel, IsHyst)
InputsModel$isUngauged <- node$id != node$donor &&
isNodeDownstream(griwrm, id, node$donor)
InputsModel$isReceiver <- node$id != node$donor &&
!isNodeDownstream(griwrm, id, node$donor)
# inUngaugedCluster: Ungauged node with downstream donor
# including reservoirs between ungauged nodes and donor
InputsModel$inUngaugedCluster <- (node$model == "Ungauged" || np$Reservoir) &&
node$id != node$donor &&
isNodeDownstream(griwrm, id, node$donor)
# isReceiver: Ungauged node with not downstream donor
InputsModel$isReceiver <- node$model == "Ungauged" &&
!isNodeDownstream(griwrm, id, node$donor)
InputsModel$gaugedId <- node$donor
InputsModel$hasUngaugedNodes <- hasUngaugedNodes(id, griwrm)
InputsModel$model <-
......
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