Commit 116ba07e authored by David's avatar David
Browse files

docs: improve Calibration messages

Showing with 11 additions and 7 deletions
+11 -7
...@@ -37,11 +37,6 @@ Calibration.GRiwrmInputsModel <- function(InputsModel, ...@@ -37,11 +37,6 @@ Calibration.GRiwrmInputsModel <- function(InputsModel,
for (id in gaugedIds) { for (id in gaugedIds) {
IM <- InputsModel[[id]] IM <- InputsModel[[id]]
if (inherits(InputsCrit[[id]], "InputsCritLavenneFunction")) {
IC <- getInputsCrit_Lavenne(id, OutputsModel, InputsCrit)
} else {
IC <- InputsCrit[[id]]
}
hasUngauged <- IM$hasUngaugedNodes hasUngauged <- IM$hasUngaugedNodes
if (hasUngauged) { if (hasUngauged) {
l <- updateParameters4Ungauged(id, l <- updateParameters4Ungauged(id,
...@@ -62,6 +57,12 @@ Calibration.GRiwrmInputsModel <- function(InputsModel, ...@@ -62,6 +57,12 @@ Calibration.GRiwrmInputsModel <- function(InputsModel,
} }
} }
if (inherits(InputsCrit[[id]], "InputsCritLavenneFunction")) {
IC <- getInputsCrit_Lavenne(id, OutputsModel, InputsCrit)
} else {
IC <- InputsCrit[[id]]
}
if (!is.null(IM$isReservoir) && IM$isReservoir & any(is.na(CalibOptions[[id]]$FixedParam))) { if (!is.null(IM$isReservoir) && IM$isReservoir & any(is.na(CalibOptions[[id]]$FixedParam))) {
stop("Parameters of node '", id, "' using `RunModel_Reservoir` can't be calibrated", stop("Parameters of node '", id, "' using `RunModel_Reservoir` can't be calibrated",
"Fix its parameters by using the command:\n", "Fix its parameters by using the command:\n",
......
...@@ -28,7 +28,7 @@ getInputsCrit_Lavenne <- function(id, OutputsModel, InputsCrit) { ...@@ -28,7 +28,7 @@ getInputsCrit_Lavenne <- function(id, OutputsModel, InputsCrit) {
AprParamR[featMod$iX4] <- AprParamR[featMod$iX4] * featMod$X4Ratio AprParamR[featMod$iX4] <- AprParamR[featMod$iX4] * featMod$X4Ratio
} }
AprParamR <- AprParamR[featMod$indexParamUngauged] AprParamR <- AprParamR[featMod$indexParamUngauged]
message("A priori parameters from node ", AprioriId, ": ", paste(round(AprParamR, 3), collapse = ", ")) message("Parameter regularization: get a priori parameters from node ", AprioriId, ": ", paste(round(AprParamR, 3), collapse = ", "))
AprCrit <- ErrorCrit(InputsCrit[[AprioriId]], OutputsModel[[AprioriId]])$CritValue AprCrit <- ErrorCrit(InputsCrit[[AprioriId]], OutputsModel[[AprioriId]])$CritValue
return(Lavenne_FUN(AprParamR, AprCrit)) return(Lavenne_FUN(AprParamR, AprCrit))
} }
...@@ -216,7 +216,7 @@ transferGRparams <- function(InputsModel, Param, donor, receiver, default_param ...@@ -216,7 +216,7 @@ transferGRparams <- function(InputsModel, Param, donor, receiver, default_param
missing_params <- setdiff(InputsModel[[receiver]]$model$indexParamUngauged, missing_params <- setdiff(InputsModel[[receiver]]$model$indexParamUngauged,
InputsModel[[donor]]$model$indexParamUngauged) InputsModel[[donor]]$model$indexParamUngauged)
if (verbose) { if (verbose) {
message("Tranfering parameters from node '", donor, "' to node '", receiver, "'") message("Tranferring parameters from node '", donor, "' to node '", receiver, "'")
} }
if (length(missing_params) > 0) { if (length(missing_params) > 0) {
if (is.null(default_param)) { if (is.null(default_param)) {
...@@ -253,6 +253,9 @@ transferGRparams <- function(InputsModel, Param, donor, receiver, default_param ...@@ -253,6 +253,9 @@ transferGRparams <- function(InputsModel, Param, donor, receiver, default_param
0.5 0.5
) )
} }
if (verbose) {
message(message("\t Param = ", paste(sprintf("%8.3f", p), collapse = ", ")))
}
return(p) return(p)
} }
......
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