tests/testthat/test_evaporation.R

In KNMI/scenarioIndices: Calculate indices for KNMI climate scenarios

context("Evaporation sums calculations")

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

context("evmk sum calc - Entire station set")

inputTemp <- KnmiRefFile("KNMI14____ref_tg___19810101-20101231_v3.2.txt")
inputRad  <- KnmiRefFile("KNMI14____ref_rsds___19810101-20101231_v3.2.txt")
ofile     <- "tmp.txt" # output file - used only temporary
var <- "tg"
regions <- MatchRegionsOnStationId(ReadInput(var, inputTemp)$header[1, -1])


test_that("evaporation sums reference", {

  tmp <- evmk_sums_reference(ofile = ofile)

  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_ref_evmk_sums.rds")
})

test_that("2030 decadal prediction", {
  scenario = "GL"

  horizon = 2030

  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario = scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14___2030_evmk_sums.rds")
})

test_that("scenarioenario WL", {
  scenario = "WL"

  horizon = 2050
  tmp <- evmkSumsRelchange(inputTemp = inputTemp, inputRad = inputRad,
                             scenario = scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_WL_2050_evmk_sums.rds")

  horizon = 2085
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_WL_2085_evmk_sums.rds")
})

test_that("scenarioenario WH", {
  scenario = "WH"

  horizon = 2050
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_WH_2050_evmk_sums.rds")

  horizon = 2085
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_WH_2085_evmk_sums.rds")
})

test_that("scenarioenario GH", {
  scenario = "GH"

  horizon = 2050
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_GH_2050_evmk_sums.rds")

  horizon = 2085
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_GH_2085_evmk_sums.rds")
})

test_that("scenarioenario GL", {
  scenario = "GL"

  horizon = 2050
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_GL_2050_evmk_sums.rds")

  horizon = 2085
  tmp <- evmkSumsRelchange(inputTemp=inputTemp, inputRad=inputRad,
                             scenario=scenario,
                             horizon = horizon,
                             regions = regions)



  expect_equal_to_reference(tmp, "regressionOutput/evaporation/KNMI14_GL_2085_evmk_sums.rds")
})