tests/testthat/test_misc.R

In bgctw/ncdfTools: Handling NetCDF Files

#require(testthat)
context('misc')

# TRUE during testthat::check()
isNotOnCRAN <- identical(Sys.getenv("NOT_CRAN"), "true")

test_that('classR2Ncdf',{
  expect_equal( classR2Ncdf(1L), "NC_INT")
  expect_equal( classR2Ncdf(1.01), "NC_DOUBLE")
  expect_equal( classR2Ncdf(33000), "NC_DOUBLE")
  expect_equal( classR2Ncdf('site'), "NC_CHAR")
})


test_that('POSIX2ncdf and back',{
  #if (require(udunits2)) {
  timesOrig <- seq(
    ISOdatetime(2010,1,1,0,0,0,tz = 'UTC'), by = "30 min", length.out = 10)
  unit <- "days since 1582-10-15 00:00"
  timesDays <- POSIX2ncdf(timesOrig, unit)
  timesPosix <- ncdf2POSIX(timesDays, unit)
  #expect_equal( timesPosix, timesOrig) # fails only on Windows i386
  # cannot trace the error by RUnit.out on winbuild.
  expect_true( all(abs(timesPosix - timesOrig) < 1)) # less than one second diff
  #}
})

test_that('.POSIX2DaysSinceOrigin and back',{
  timesOrig <- seq(
    ISOdatetime(2010,1,1,0,0,0,tz = 'UTC'), by = "30 min", length.out = 10)
  timesDays <- ncdfTools:::.POSIX2DaysSinceOrigin(timesOrig)
  unit <- "days since 1582-10-15 00:00"
  timesDaysExp <- POSIX2ncdf(timesOrig, unit)
  expect_equal( timesDays, timesDaysExp)
  timesPosix <- ncdfTools:::.daysSinceOrigin2POSIX(timesDays)
  expect_equal( timesPosix, timesOrig) 
})

test_that('read and write ncdf time',{
  nRec <- 30L
  times <- seq(ISOdatetime(2010,1,1,0,0,0, tz = "UTC")
               , by = "30 min", length.out = nRec)
  dsTest <- data.frame(
    time = times
    , soilResp = rnorm(nRec)
    , sdSoilResp = 0.1
  )
  tmpDir <- tempdir()
  #tmpDir <- "tmp"
  unitsOrig <- c('g/m2','g2/m4')
  ncFileDs <- createStdNcdfFile(
    data = dsTest
    , file.name = file.path(tmpDir,"dsTestTime.nc")
    , lat.values = 53.0
    , long.values = 13
    , time.values = times
    , units = unitsOrig
  )
  ans <- readNcdfTime(ncFileDs)  
  expect_equal(ans, times)
  infoVars <- infoNcdfVars(ncFileDs)
  units <- infoVars$unit[match(names(dsTest)[-1], infoVars$name)]
  expect_equal(units, unitsOrig)
  #
  # modify file to use different origin
  origin2 <-  ISOdatetime(1970,1,1,0,0,0, tz = "UTC")
  ans <- modifyNcdfWriteTime(ncFileDs, origin = origin2)
  expect_equal( ans, ncdfTools:::.POSIX2DaysSinceOrigin(times, origin = origin2))
  ansInfo <- infoNcdfAtts(ncFileDs, var.id = "time")
  expect_equal(ansInfo["units","value"], "days since 1970-01-01")
  ans <- readNcdfTime(ncFileDs)  
  expect_equal(ans, times)
  #
  # convert to standard time units again
  ans <- modifyNcdfWriteTime(ncFileDs)
  expect_equal( ans, ncdfTools:::.POSIX2DaysSinceOrigin(times))
  ansInfo <- infoNcdfAtts(ncFileDs, var.id = "time")
  expect_equal(ansInfo["units","value"], "days since 1582-10-15")
  ans <- readNcdfTime(ncFileDs)  
  expect_equal(ans, times)
})


test_that('defining a factor dimension',{
  #require(testthat)
  nRec <- 30L
  times <- seq(ISOdatetime(2010,1,1,0,0,0, tz = "UTC")
               , by = "30 min", length.out = nRec)
  dsTest <- data.frame(
    time = times
    , soilResp = rnorm(nRec)
  )
  tmpDir <- tempdir()
  #tmpDir <- "tmp"
  ncFile <- createStdNcdfFile(
    "soilResp"
    , file.name = file.path(tmpDir,"soilResp.1.1.30.nc")
    , lat.values = 53.0
    , long.values = 13
    , time.values = dsTest$time
    , units = 'g/m2'
    , data = dsTest$soilResp
  )
  createFactorDim(ncFile,"treatment",c("ctrl","Nadd","NPadd"))
  infoVars <- infoNcdfVars(ncFile, dimvars = TRUE)
  expect_true("treatment" %in% infoVars$name)
  infoAtts <- infoNcdfAtts(ncFile, "treatment")
})

test_that('.parseLevelAttr',{
  levelsOrig <- c("ctrl","Nadd","NPadd")
  nLevel <- length(levelsOrig)
  levelsAttr <- paste(paste(1:nLevel, levelsOrig, sep = "="), collapse = ",")
  levelsParsed <- ncdfTools:::.parseLevelAttr(levelsAttr)
  expect_equal(levelsParsed, levelsOrig)
  expect_equal(ncdfTools:::.parseLevelAttr("1=bla"), "bla")
})