tests/testthat/test-runLWFB90_errors.R

In pschmidtwalter/LWFBrook90R: Simulate Evapotranspiration and Soil Moisture with the SVAT Model LWF-Brook90

# Set up the input data
data("slb1_soil")
data("slb1_meteo")
opts <- set_optionsLWFB90(startdate = as.Date("2002-06-01"), enddate = as.Date("2002-06-30"))
soil <- cbind(slb1_soil, hydpar_wessolek_tab(texture = slb1_soil$texture))

parms <- set_paramLWFB90()

test_that("timelimit works",{
expect_error(
  run_LWFB90(options_b90 = set_optionsLWFB90(startdate = as.Date("1990-01-01")),
          param_b90 = parms,
          climate = slb1_meteo,
          soil = soil,
          rtrn_output = FALSE,
          rtrn_input = FALSE,
          timelimit = 0.1)
)
  })

soil_bad <- soil
soil_bad$ths[10] <- NA
parms_badsoil <- parms
parms_badsoil$soil_nodes <- soil_to_param(soil_bad)[[1]]
parms_badsoil$soil_materials <- soil_to_param(soil_bad)[[2]]

test_that("NAs in soil are refused",{
  expect_error(
    run_LWFB90(options_b90 = set_optionsLWFB90(startdate = as.Date("1990-01-01")),
               param_b90 = parms,
               climate = slb1_meteo,
               soil = soil_bad,
               rtrn_output = FALSE,
               rtrn_input = FALSE)
  )
  expect_error(
    run_LWFB90(options_b90 = set_optionsLWFB90(startdate = as.Date("1990-01-01")),
               param_b90 = parms_badsoil,
               climate = slb1_meteo,
               rtrn_output = FALSE,
               rtrn_input = FALSE)
  )


  })


# manipulate climate
clim_bad <- slb1_meteo[data.table::year(slb1_meteo$dates) == 2002,]
clim_bad$dates[170] <- as.Date("2002-06-30")
clim_na <- slb1_meteo[data.table::year(slb1_meteo$dates) == 2002,]
clim_na[160:165,c(5,8)] <- NA

opts <- set_optionsLWFB90(startdate = as.Date("2002-06-01"), enddate = as.Date("2002-06-30"))
parms <- set_paramLWFB90()

test_that("errorhandling works",{
  expect_error(
    run_LWFB90(options_b90 = opts,
              param_b90 = set_paramLWFB90(),
              climate = clim_bad,
              soil = soil,
              rtrn_input = FALSE, rtrn_output = FALSE)
  )
  expect_error(
    run_LWFB90(options_b90 = opts,
               param_b90 = set_paramLWFB90(),
               climate = clim_na,
               soil = soil,
               rtrn_input = FALSE)
  )

  expect_error(
    run_LWFB90(options_b90 = opts,
              param_b90 = set_paramLWFB90(psiini = c(rep(-6.3, nrow(soil)-1), 999)),
              climate = slb1_meteo,
              soil = soil,
              rtrn_input = FALSE, rtrn_output = FALSE)
  )
})