tests/testthat/test-precipitation_centr.R

In MartinRoth/knmitransformer: Transformation program

context("Transformation of precipitation")

library(futile.logger)
flog.threshold(DEBUG)
flog.appender(appender.file("knmitransformer_precipitation.log"))
library(data.table)

context("rr transformation (centr) - Entire station set")

subscenario <- "centr"
input       <- ReadInput("rr",
    KnmiRefFile("KNMI14____ref_rrcentr___19810101-20101231_v3.2.txt"))
ofile       <- NA
rounding    <- FALSE

test_that("2030 decadal prediction", {
  scenario <- "GL"
  horizon  <- 2030
  tmp <- TransformPrecip(input = KnmiRefFile("KNMI14____ref_rrcentr___19810101-20101231_v3.2.txt"), #nolint
                         ofile = "tmp.txt",
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14___2030_rr_centr.rds")

  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = TRUE)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14___2030_rr_centr_rounded.rds")
})

test_that("Scenario GL", {
  scenario <- "GL"
  horizon  <- 2050
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_GL_2050_rr_centr.rds")

  horizon <- 2085
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_GL_2085_rr_centr.rds")
})

test_that("Scenario GH", {
  scenario <- "GH"
  horizon  <- 2050
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_GH_2050_rr_centr.rds")

  horizon <- 2085
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_GH_2085_rr_centr.rds")
})

test_that("Scenario WH", {
  scenario <- "WH"
  horizon <- 2050
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_WH_2050_rr_centr.rds")

  horizon <- 2085
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_WH_2085_rr_centr.rds")
})

test_that("Scenario WL", {
  scenario <- "WL"
  horizon  <- 2050
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_WL_2050_rr_centr.rds")

  horizon <- 2085
  tmp <- TransformPrecip(input = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14_WL_2085_rr_centr.rds")
})