inst/data_generation/EFSA_GW_interception.R

library(here)

bbch <- paste0(0:9, "x")
crops <- c(
  "Beans (field + vegetable)",
  "Peas",
  "Summer oilseed rape", "Winter oilseed rape",
  "Tomatoes",
  "Spring cereals", "Winter cereals")
EFSA_GW_interception_2014 <- matrix(NA, length(crops), length(bbch),
  dimnames = list(Crop = crops, BBCH = bbch))
EFSA_GW_interception_2014["Beans (field + vegetable)", ] <-
  c(0, 0.25, rep(0.4, 2), rep(0.7, 5), 0.8)
EFSA_GW_interception_2014["Peas", ] <-
  c(0, 0.35, rep(0.55, 2), rep(0.85, 5), 0.85)
EFSA_GW_interception_2014["Summer oilseed rape", ] <-
  c(0, 0.4, rep(0.8, 2), rep(0.8, 5), 0.9)
EFSA_GW_interception_2014["Winter oilseed rape", ] <-
  c(0, 0.4, rep(0.8, 2), rep(0.8, 5), 0.9)
EFSA_GW_interception_2014["Tomatoes", ] <-
  c(0, 0.5, rep(0.7, 2), rep(0.8, 5), 0.5)
EFSA_GW_interception_2014["Spring cereals", ] <-
  c(0, 0, 0.2, 0.8, rep(0.9, 3), rep(0.8, 2), 0.8)
EFSA_GW_interception_2014["Winter cereals", ] <-
  c(0, 0, 0.2, 0.8, rep(0.9, 3), rep(0.8, 2), 0.8)

save(EFSA_GW_interception_2014,
  file = here("data/EFSA_GW_interception_2014.RData"))
jranke/pfm documentation built on March 14, 2024, 3:18 a.m.