tests/testthat/test-et_analyze.R

In jkennel/earthtide: Parallel Implementation of 'ETERNA 3.40' for Prediction and Analysis of Earth Tides

context("test-et_analyze")

test_that("dimensions correct", {
  tms <- as.POSIXct('1990-01-01', tz = 'UTC') + c(0, 3600)
  
  wave_groups = data.frame(start = 0, end = 8)
  
  et <- Earthtide$new(utc = tms, 
                      latitude = 52.3868,
                      longitude = 9.7144,
                      elevation = 110,
                      gravity = 9.8127, 
                      cutoff = 1.0e-10,
                      catalog = 'hw95s',
                      wave_groups = wave_groups)
  
  out <- et$analyze(method = 'gravity')$tide()
  expect_equal(NROW(out), length(tms))
  expect_equal(NCOL(out), 3)
  
  
  tms <- as.POSIXct('1990-01-01', tz = 'UTC')
  wave_groups = data.frame(start = c(0, 1.5), end = c(1.4, 8))
  et <- Earthtide$new(utc = tms, 
                      latitude = 52.3868,
                      longitude = 9.7144,
                      elevation = 110,
                      gravity = 9.8127, 
                      cutoff = 1.0e-10,
                      wave_groups = wave_groups)
  
  out <- et$analyze(method = 'gravity')$tide()
  expect_equal(nrow(out), length(tms))
  expect_equal(ncol(out), 5)
  
  
  
  
  tms <- as.POSIXct('1990-01-01', tz = 'UTC') + 0:(24*31) * 3600
  et <- Earthtide$new(utc = tms, 
                      latitude = 52.3868,
                      longitude = 9.7144,
                      elevation = 110,
                      gravity = 9.8127, 
                      cutoff = 1.0e-10,
                      wave_groups = wave_groups)
  
  out <- et$analyze(method = 'gravity', astro_update = 24)$tide()
  expect_equal(nrow(out), length(tms))
  expect_equal(ncol(out), 5)
  
})
  
  
test_that("scaling is correct", {
  
  # scale 
  tms <- as.POSIXct('1990-01-01', tz = 'UTC') + seq(0, 86400*14, 60)
  
  wave_groups = eterna_wavegroups
  wave_groups <- na.omit(wave_groups[wave_groups$time == 'all',])
  
  et <- Earthtide$new(utc = tms, 
                      longitude = -118.67,
                      latitude = 34.23,
                      elevation = 500,
                      cutoff = 1.0e-5,
                      catalog = 'hw95s',
                      wave_groups = wave_groups)
  
  out1 <- et$analyze(method = 'gravity',  scale = FALSE, astro_update = 1)$tide()
  
  out2 <- et$analyze(method = 'gravity',  scale = FALSE, astro_update = 10)$tide()
  
  out3 <- et$analyze(method = 'gravity',  scale = TRUE, astro_update = 1)$tide()
  
  out4 <- et$analyze(method = 'gravity',  scale = TRUE, astro_update = 10)$tide()
  
  
  
  expect_equal(out1, out2, tolerance = 1e-7)
  expect_equal(out3, out4, tolerance = 1e-7)
  expect_false(isTRUE(all.equal(out1, out3, tolerance = 1e-2)))
  expect_false(isTRUE(all.equal(out2, out4, tolerance = 1e-2)))
  
  
  
})