Skip to content
GitLab
Select Git revision
  • code_refactor
  • board_v4
  • EGU23
  • revert-559ef012
  • MUX_board_v2024
  • master default
  • battery_charger
  • mqtt
  • dev2
  • Ohmpi_1.0x
  • v2023.0.0
  • V1.02
  • v1.01
  • v1
Find file BlameHistoryPermalink
Forked from reversaal / OhmPi
Source project has a limited visibility.
helper_1_RunModel.R 3.53 KiB
1 
#' Prepare useful variables for GRiwrm tests
2 
#'
3 
#' @return [environment] with the variables (See examples section)
4 
#' @noRd
5 
#'
6 
#' @examples
7 
#' # data set up
8 
#' e <- setupRunModel()
9 
# variables are copied from environment 'e' to the current environment
10 
# https://stackoverflow.com/questions/9965577/r-copy-move-one-environment-to-another
11 
#' for (x in ls(e)) assign(x, get(x, e))
12 
#'
13 
setupRunModel <-
14 
  function(runInputsModel = TRUE,
15 
           runRunOptions = TRUE,
16 
           runRunModel = TRUE,
17 
           griwrm = NULL,
18 
           Qinf = NULL,
19 
           Qrelease = NULL,
20 
           Qmin = NULL,
21 
           IsHyst = FALSE,
22 
           ParamMichel = getDefaultParamMichel()) {
23
24 
    data(Severn)
25
26 
    # Format observation
27 
    BasinsObs <- Severn$BasinsObs
28 
    DatesR <- BasinsObs[[1]]$DatesR
29 
    PrecipTot <-
30 
      cbind(sapply(BasinsObs, function(x) {
31 
        x$precipitation
32 
      }))
33 
    PotEvapTot <- cbind(sapply(BasinsObs, function(x) {
34 
      x$peti
35 
    }))
36 
    Qobs <- cbind(sapply(BasinsObs, function(x) {
37 
      x$discharge_spec
38 
    }))
39
40 
    # Set network
41 
    if (is.null(griwrm)) {
42 
      nodes <- loadSevernNodes()
43 
      griwrm <-
44 
        CreateGRiwrm(nodes)
45 
    }
46
47 
    # Convert meteo data to SD (remove upstream areas)
48 
    Precip <- ConvertMeteoSD(griwrm, PrecipTot)
49 
    PotEvap <- ConvertMeteoSD(griwrm, PotEvapTot)
50 
    if (IsHyst) {
51 
      TempMean <- PotEvap+5 # Fake temperatures
52 
    } else {
53 
      TempMean <- NULL
54 
    }
55
56 
    # set up inputs
57 
    if (!runInputsModel)
58 
      return(environment())
59 
    InputsModel <-
60 
      suppressWarnings(CreateInputsModel(griwrm, DatesR, Precip, PotEvap,
61 
                                         TempMean = TempMean,
62 
                                         Qinf = Qinf,
63 
                                         Qrelease = Qrelease,
64 
                                         Qmin = Qmin,
65 
                                         IsHyst = IsHyst))
66
67 
    # RunOptions
68 
    if (!runRunOptions)
69 
      return(environment())
70 
    e <- setupRunOptions(InputsModel)
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137
for (x in ls(e)) assign(x, get(x, e)) rm(e) # RunModel.GRiwrmInputsModel if (!runRunModel) return(environment()) OM_GriwrmInputs <- RunModel(InputsModel, RunOptions = RunOptions, Param = ParamMichel) return(environment()) } setupRunOptions <- function(InputsModel) { nTS <- 365 IndPeriod_Run <- seq(length(InputsModel[[1]]$DatesR) - nTS + 1, length(InputsModel[[1]]$DatesR)) IndPeriod_WarmUp = seq(IndPeriod_Run[1] - 365, IndPeriod_Run[1] - 1) RunOptions <- CreateRunOptions(InputsModel, IndPeriod_WarmUp = IndPeriod_WarmUp, IndPeriod_Run = IndPeriod_Run) return(environment()) } getDefaultParamMichel <- function() { list( `54057` = c( 0.781180650559296, 74.4247133147015, -1.26590474908235, 0.96012365697022, 2.51101785373787 ), `54032` = c( 0.992743594649893, 1327.19309205366, -0.505902143697436, 7.91553615210291, 2.03604525989572 ), `54001` = c( 1.03, 24.7790862245877, -1.90430150145153, 21.7584023961971, 1.37837837837838 ), `54095` = c( 256.844150254651, 0.0650458497009288, 57.523675209819, 2.71809513102128 ), `54002` = c( 419.437754485522, 0.12473266292168, 13.0379482833606, 2.12230907892238 ), `54029` = c( 219.203385553954, 0.389211590110934, 48.4242150713452, 2.00300300300301 ) ) }