tests/testthat/test_lisst.R

In dankelley/oce: Analysis of Oceanographic Data

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

test_that("as.lisst()", {
    set.seed(1333334L)
    t <- seq(0, 6, 1/15) * 3600 + as.POSIXct("2012-01-01 00:00:00", tz="UTC")
    n <- length(t)
    pre <- 5 + sin(as.numeric(t - t[1]) / 12.4 / 3600 * 2 * pi) + rnorm(n, sd=0.01)
    dpdt <- c(0, diff(pre))
    tem <- 10 + 5 * sin(as.numeric(t - t[1]) / 24 / 3600 * 2 * pi) + cumsum(rnorm(n, sd=0.2))
    C <- (dpdt + rnorm(n, sd=0.1) + cumsum(rnorm(n, sd=0.5)))^2 * 2
    sd <- rep(1, length.out=32) + (1:32) / 100
    data <- matrix(nrow=n, ncol=42)
    for (i in 1:32) {
        fake <- abs(C * (1 + i / 5) + cumsum(rnorm(n, sd=sd[i]))) / 100
        data[, i] <- fake
    }
    data[, 33] <- rep(0, n) # lts
    data[, 34] <- rep(4, n) + 0.01 * cumsum(rnorm(n, 0.05)) # voltage
    data[, 35] <- rep(0.07, n) # aux
    data[, 36] <- runif(n, 3.9, 4.1) # lrs
    data[, 37] <- pre
    data[, 38] <-tem
    tt <- as.POSIXlt(t)
    data[, 39] <- 100 * tt$yday + tt$hour
    data[, 40] <- 100 * tt$min + tt$sec
    data[, 41] <- abs((pre - min(pre)) / diff(range(pre)) + cumsum(rnorm(n, sd=0.05))) # transmission
    data[, 42] <- 40 - 20*data[, 41] # beam
    lisst <- as.lisst(data, filename="(constructed)", year=2012, "UTC")
    expect_equal(lisst[["beam"]], 40-20*data[, 41])
})