tests/testthat/test-precipitation_lower.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 (lower) - Entire station set")

subscenario <- "lower"
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 = input,
                         ofile = ofile,
                         scenario = scenario,
                         horizon = horizon,
                         subscenario = subscenario,
                         rounding = rounding)
  expect_equal_to_reference(tmp,
      "regressionOutput/precipitation/KNMI14___2030_rr_lower.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_lower.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_lower.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_lower.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_lower.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_lower.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_lower.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_lower.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_lower.rds")
})