tests/testthat/test_datasets.R

In oce: Analysis of Oceanographic Data

## vim:textwidth=80:expandtab:shiftwidth=4:softtabstop=4
library(oce)

## FIXME: commented-out data() lines indicate things to be done.

## Part 1: various tests, some of them ad-hock

##data("adp")
##data("adv")
##data("argo")

test_that("cm", {
    data("cm")
    S <- cm[["salinity"]]
    S1 <- swSCTp(cm)
    expect_equal(S, S1, tolerance=0.001)
    S2a <- swSCTp(cm[["conductivity"]], cm[["temperature"]], cm[["pressure"]], conductivityUnit=cm[["conductivity unit"]])
    expect_equal(S1, S2a)
    S2b <- swSCTp(cm[["conductivity"]], cm[["temperature"]], cm[["pressure"]], conductivityUnit=cm[["conductivityUnit"]])
    expect_equal(S1, S2b)
    S2c <- swSCTp(cm[["conductivity"]], cm[["temperature"]], cm[["pressure"]], conductivityUnit=as.character(cm[["conductivityUnit"]]$unit))
    expect_equal(S1, S2c)
})

##data("coastlineWorld")

test_that("ctd", {
    data("ctd")
    expect_equal(ctd[["latitude"]],   44.6842666666667)
    expect_equal(ctd[["longitude"]], -63.6438833333333)
    # next two lines test issues 1460 and 1547
    expect_equal(ctd[["time"]][1], as.POSIXct("2003-10-15 15:40:47", tz="UTC"))
    expect_equal(length(ctd[["time"]]), length(ctd[["pressure"]]))
    # units are checked in test_accessors.R
    expect_equal(ctd[["pressureType"]], "sea")
})

##data("ctdRaw")
##data("echosounder")
##data("lisst")
##data("lobo")
##data("rsk")
##data("met")

test_that("sealevel", {
    data("sealevel")
    expect_equal(sealevel[["latitude"]],   44.666667)
    expect_equal(sealevel[["longitude"]], -63.583333)
    expect_equal(sealevel[["stationNumber"]], 490)
    expect_equal(sealevel[["GMTOffset"]], 0)
})

##data("sealevelTuktoyaktuk")

test_that("sealevel", {
    data("section")
    expect_equal(section[["sectionId"]], "a03")
    expect_equal(length(section@data$station), 124)
    expect_equal(section@data$station[[1]]@metadata$station, "3")
    expect_equal(section@data$station[[1]][["station"]], "3")
    expect_equal(section@data$station[[1]][["latitude"]], 36.8758)
    expect_equal(section@data$station[[1]][["longitude"]], -8.5263)
})

##data("tidedata")
test_that("topoWorld", {
    data("topoWorld")
    expect_equal(range(topoWorld[["longitude"]]), 179.75*c(-1, 1))
    expect_equal(range(topoWorld[["latitude"]]), 89.75*c(-1, 1))
})
##data("wind")


## Part 2: check units

## below constructed from the files in oce/data
test_that("units", {
    names <- c("adp", "adv", "argo", "cm", "coastlineWorld", "ocecolors", "ctd",
        "ctdRaw", "echosounder", "landsat", "lisst", "lobo", "met", "rsk",
        "sealevel", "sealevelTuktoyaktuk", "section", "tidedata",
        "topoWorld", "wind")
    for (name in names) {
        data(list=name)
        if (name == "section") {
            x <- get(name)
            x <- x[["station", 1]]
        } else {
            x <- get(name)
        }
        if ("metadata" %in% slotNames(x) && "units" %in% names(x@metadata)) {
            units <- x@metadata$units
            if (length(units) > 0) {
                unitsNames <- names(units)
                for (i in seq_along(units)) {
                    this <- units[[i]]
                    expect_equal(2, length(this))
                    expect_equal(2, sum(c("unit", "scale") %in% names(this)))
                    expect_true(is.expression(this$unit))
                    expect_true(is.character(this$scale))
                }
            }
        }
    }
})