library(pfm)
context("Exposit runoff calculations")
test_that("Runoff PECsw are as in Exposit 3.01", {
# Expected results are from the Exposit 3.01, downloaded 2017-10-27
# 100 g/ha, Koc = 1000 L/kg, DT50 = 1000 days
PEC_spreadsheet_1 <- data.frame(dissolved = c(0.706, 0.612, 0.525, 0.367),
total = c(0.783, 0.679, 0.546, 0.377))
PEC_pfm_1 <- round(PEC_sw_exposit_runoff(100, Koc = 1000, DT50 = 1000)$PEC_sw_runoff, 3)[c("dissolved", "total")]
expect_equivalent(PEC_spreadsheet_1, PEC_pfm_1)
# 10 g/ha, Koc = 300000 L/kg, DT50 = 10 days
PEC_spreadsheet_2 <- data.frame(dissolved = c(0, 0, 0, 0),
total = c(0.141, 0.122, 0.039, 0.018))
PEC_pfm_2 <- round(PEC_sw_exposit_runoff(10, Koc = 300000, DT50 = 10)$PEC_sw_runoff, 3)[c("dissolved", "total")]
expect_equivalent(PEC_spreadsheet_2, PEC_pfm_2)
# 200 g/ha, Koc = 30 L/kg, DT50 = 100 days
PEC_spreadsheet_3 <- data.frame(dissolved = c(1.138, 0.986, 0.845, 0.592),
total = c(1.138, 0.986, 0.845, 0.592))
PEC_pfm_3 <- round(PEC_sw_exposit_runoff(200, Koc = 30, DT50 = 100)$PEC_sw_runoff, 3)[c("dissolved", "total")]
expect_equivalent(PEC_spreadsheet_3, PEC_pfm_3)
})
