Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
CreateIniStates.R 10.94 KiB
1 
CreateIniStates <- function(FUN_MOD, InputsModel, IsHyst = FALSE, IsIntStore = FALSE,
2 
                            ProdStore = 350, RoutStore = 90, ExpStore = NULL, IntStore = NULL,
3 
                            UH1 = NULL, UH2 = NULL,
4 
                            GCemaNeigeLayers = NULL, eTGCemaNeigeLayers = NULL,
5 
                            GthrCemaNeigeLayers = NULL, GlocmaxCemaNeigeLayers = NULL,
6 
                            SD = NULL,
7 
                            verbose = TRUE) {
8
9
10 
  ObjectClass <- NULL
11
12 
  UH1n <- 20L
13 
  UH2n <- UH1n * 2L
14
15 
  nameFUN_MOD <- as.character(substitute(FUN_MOD))
16 
  FUN_MOD <- match.fun(FUN_MOD)
17
18 
  ## check FUN_MOD
19 
  BOOL <- FALSE
20 
  if (identical(FUN_MOD, RunModel_GR4H) | identical(FUN_MOD, RunModel_GR5H)) {
21 
    ObjectClass <- c(ObjectClass, "GR", "hourly")
22 
    BOOL <- TRUE
23 
  }
24 
  if (identical(FUN_MOD, RunModel_GR4J) |
25 
      identical(FUN_MOD, RunModel_GR5J) |
26 
      identical(FUN_MOD, RunModel_GR6J)) {
27 
    ObjectClass <- c(ObjectClass, "GR", "daily")
28 
    BOOL <- TRUE
29 
  }
30 
  if (identical(FUN_MOD, RunModel_GR2M)) {
31 
    ObjectClass <- c(ObjectClass, "GR", "monthly")
32 
    BOOL <- TRUE
33 
  }
34 
  if (identical(FUN_MOD, RunModel_GR1A)) {
35 
    stop("'RunModel_GR1A' does not require 'IniStates' object")
36 
  }
37 
  if (identical(FUN_MOD, RunModel_CemaNeige)) {
38 
    ObjectClass <- c(ObjectClass, "CemaNeige", "daily")
39 
    BOOL <- TRUE
40 
  }
41 
  if (identical(FUN_MOD, RunModel_CemaNeigeGR4H) | identical(FUN_MOD, RunModel_CemaNeigeGR5H)) {
42 
    ObjectClass <- c(ObjectClass, "GR", "CemaNeige", "hourly")
43 
    BOOL <- TRUE
44 
  }
45 
  if (identical(FUN_MOD, RunModel_CemaNeigeGR4J) |
46 
      identical(FUN_MOD, RunModel_CemaNeigeGR5J) |
47 
      identical(FUN_MOD, RunModel_CemaNeigeGR6J)) {
48 
    ObjectClass <- c(ObjectClass, "GR", "CemaNeige", "daily")
49 
    BOOL <- TRUE
50 
  }
51 
  if (!BOOL) {
52 
    stop("incorrect 'FUN_MOD' for use in 'CreateIniStates'")
53 
  }
54 
  if (!"CemaNeige" %in% ObjectClass & IsHyst) {
55 
    stop("'IsHyst' cannot be TRUE if CemaNeige is not used in 'FUN_MOD'")
56 
  }
57 
  if (!(identical(FUN_MOD, RunModel_GR5H) | identical(FUN_MOD, RunModel_CemaNeigeGR5H)) & IsIntStore) {
58 
    stop("'IsIntStore' cannot be TRUE if GR5H is not used in 'FUN_MOD'")
59 
  }
60
61 
  ## check InputsModel
62 
  if (!inherits(InputsModel, "InputsModel")) {
63 
    stop("'InputsModel' must be of class 'InputsModel'")
64 
  }
65 
  if ("GR" %in% ObjectClass & !inherits(InputsModel, "GR")) {
66 
    stop("'InputsModel' must be of class 'GR'")
67 
  }
68 
  if ("CemaNeige" %in% ObjectClass &
69 
      !inherits(InputsModel, "CemaNeige")) {
70 
    stop("'InputsModel' must be of class 'CemaNeige'")