tests/testthat/test_read.r

In scheidan/SWMMR: R interface for EPA Storm Water Management Model (SWMM)

## -------------------------------------------------------
## Andreas Scheidegger
## andreas.scheidegger@eawag.ch
## -------------------------------------------------------


library(SWMMR)


context("Test reading of SWMM binary output files")

## calculate a cheap check sum
checksum <- function(x) sum(x*rep(1:7, length=length(x)))


ff <- openSWMMOutput("../testdata/output")


test_that("Check if file properties are correct", {
  expect_equal(ff$numSubc, 1)
  expect_equal(ff$numNode, 1)
  expect_equal(ff$numLink, 0)
  expect_equal(ff$numPoll, 0)
})


test_that("Check if subcatchments information is correct", {
  expect_equal(length(readSubcatchments(ff, names=c("S1"))), ff$numSubcVars)
  
  ## checksums.S1 <- c(rainfall=232.99, `snow depth`=0, `evaporation losses`=0,
  ##                   `infiltration losses`=223.4777,
  ##                   runoff=333.9936, `groundwater outflow`=0,
  ##                   `groundwater water table elevation`=0,
  ##                   `unsaturated zone moisture content`=0)
  
  ## for(var in names(checksums.S1)){
  ##   print(var)
  ##   print(checksum(readSubcatchments(ff, names=c("S1"), variables=var)[[1]]))
  ##   print(checksums.S1[var])
  ##   expect_equivalent(checksum(readSubcatchments(ff, names=c("S1"), variables=var)[[1]]),
  ##                     checksums.S1[var])
  ## }
  
})


test_that("Check if node information is correct", {
  expect_equal(length(readNodes(ff, names=c("O1"))), ff$numNodeVar) 
})


test_that("Check if link information is correct", {
  ## test file conatins no link
  ## expect_equal(length(readLinks(ff, names=c("L1"))), ff$numLinkPropSaved)
})


test_that("Check if system information is correct", {
  expect_equal(length(readSystem(ff)), ff$numSysVars)
})


closeSWMMOutput(ff)