Skip to content
GitLab
Select Git revision
  • develop
  • main default
Find file BlameHistoryPermalink
CreateInputsModel.R 10.47 KiB
1 
CreateInputsModel <- function(FUN_MOD,
2 
                              DatesR,
3 
                              Precip, PrecipScale = TRUE,
4 
                              PotEvap = NULL,
5 
                              TempMean = NULL, TempMin = NULL, TempMax = NULL,
6 
                              ZInputs = NULL, HypsoData = NULL, NLayers = 5,
7 
                              verbose = TRUE) {
8
9
10 
    ObjectClass <- NULL
11
12 
    FUN_MOD <- match.fun(FUN_MOD)
13
14 
    ##check_FUN_MOD
15 
    BOOL <- FALSE
16 
    if (identical(FUN_MOD, RunModel_GR4H)) {
17 
      ObjectClass <- c(ObjectClass, "hourly", "GR")
18
19 
      TimeStep <- as.integer(60 * 60)
20
21 
      BOOL <- TRUE
22 
    }
23 
    if (identical(FUN_MOD, RunModel_GR4J) |
24 
        identical(FUN_MOD, RunModel_GR5J) |
25 
        identical(FUN_MOD, RunModel_GR6J)) {
26 
      ObjectClass <- c(ObjectClass, "daily", "GR")
27
28 
      TimeStep <- as.integer(24 * 60 * 60)
29
30 
      BOOL <- TRUE
31 
    }
32 
    if (identical(FUN_MOD, RunModel_GR2M)) {
33 
      ObjectClass <- c(ObjectClass, "GR", "monthly")
34
35 
      TimeStep <- as.integer(c(28, 29, 30, 31) * 24 * 60 * 60)
36
37 
      BOOL <- TRUE
38 
    }
39 
    if (identical(FUN_MOD, RunModel_GR1A)) {
40 
      ObjectClass <- c(ObjectClass, "GR", "yearly")
41
42 
      TimeStep <- as.integer(c(365, 366) * 24 * 60 * 60)
43
44 
      BOOL <- TRUE
45 
    }
46 
    if (identical(FUN_MOD, RunModel_CemaNeige)) {
47 
      ObjectClass <- c(ObjectClass, "daily", "CemaNeige")
48
49 
      TimeStep <- as.integer(24 * 60 * 60)
50
51 
      BOOL <- TRUE
52 
    }
53 
    if (identical(FUN_MOD, RunModel_CemaNeigeGR4J) |
54 
        identical(FUN_MOD, RunModel_CemaNeigeGR5J) |
55 
        identical(FUN_MOD, RunModel_CemaNeigeGR6J)) {
56 
      ObjectClass <- c(ObjectClass, "daily", "GR", "CemaNeige")
57
58 
      TimeStep <- as.integer(24 * 60 * 60)
59
60 
      BOOL <- TRUE
61 
    }
62 
    if (identical(FUN_MOD, RunModel_CemaNeigeGR4H)) {
63 
      ObjectClass <- c(ObjectClass, "hourly", "GR", "CemaNeige")
64
65 
      TimeStep <- as.integer(60 * 60)
66
67 
      BOOL <- TRUE
68 
    }
69 
    if (!BOOL) {
70 
      stop("incorrect 'FUN_MOD' for use in 'CreateInputsModel'")
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
} ##check_arguments if ("GR" %in% ObjectClass | "CemaNeige" %in% ObjectClass) { if (is.null(DatesR)) { stop("'DatesR' is missing") } if (!"POSIXlt" %in% class(DatesR) & !"POSIXct" %in% class(DatesR)) { stop("'DatesR' must be defined as 'POSIXlt' or 'POSIXct'") } if (!"POSIXlt" %in% class(DatesR)) { DatesR <- as.POSIXlt(DatesR) } if (!difftime(tail(DatesR, 1), tail(DatesR, 2), units = "secs")[[1]] %in% TimeStep) { TimeStepName <- grep("hourly|daily|monthly|yearly", ObjectClass, value = TRUE) stop(paste0("the time step of the model inputs must be ", TimeStepName, "\n")) } if (any(duplicated(DatesR))) { stop("'DatesR' must not include duplicated values") } LLL <- length(DatesR) } if ("GR" %in% ObjectClass) { if (is.null(Precip)) { stop("Precip is missing") } if (is.null(PotEvap)) { stop("'PotEvap' is missing") } if (!is.vector(Precip) | !is.vector(PotEvap)) { stop("'Precip' and 'PotEvap' must be vectors of numeric values") } if (!is.numeric(Precip) | !is.numeric(PotEvap)) { stop("'Precip' and 'PotEvap' must be vectors of numeric values") } if (length(Precip) != LLL | length(PotEvap) != LLL) { stop("'Precip', 'PotEvap' and 'DatesR' must have the same length") } } if ("CemaNeige" %in% ObjectClass) { if (is.null(Precip)) { stop("'Precip' is missing") } if (is.null(TempMean)) { stop("'TempMean' is missing") } if (!is.vector(Precip) | !is.vector(TempMean)) { stop("'Precip' and 'TempMean' must be vectors of numeric values") } if (!is.numeric(Precip) | !is.numeric(TempMean)) { stop("'Precip' and 'TempMean' must be vectors of numeric values") } if (length(Precip) != LLL | length(TempMean) != LLL) { stop("'Precip', 'TempMean' and 'DatesR' must have the same length") } if (is.null(TempMin) != is.null(TempMax)) { stop("'TempMin' and 'TempMax' must be both defined if not null") } if (!is.null(TempMin) & !is.null(TempMax)) { if (!is.vector(TempMin) | !is.vector(TempMax)) { stop("'TempMin' and 'TempMax' must be vectors of numeric values") } if (!is.numeric(TempMin) | !is.numeric(TempMax)) { stop("'TempMin' and 'TempMax' must be vectors of numeric values") } if (length(TempMin) != LLL | length(TempMax) != LLL) { stop("'TempMin', 'TempMax' and 'DatesR' must have the same length") } } if (!is.null(HypsoData)) {