tests/testthat/test-pole_lod.R

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

context("test-pole_lod")

test_that("pole and LOD tides work", {
  tms <- as.POSIXct("2011-01-01", tz = "UTC")
  wave_groups <- data.frame(start = 0, end = 8)

  et <- Earthtide$new(
    utc = tms,
    latitude = 14.8861,
    longitude = 103.5203,
    elevation = 20,
    gravity = 0,
    cutoff = 1.0e-10,
    catalog = "hw95s",
    wave_groups = wave_groups
  )
  et$pole_tide()
  et$lod_tide()

  tide <- et$tide()

  # This needs to be checked
  expect_equal(tide$pole_tide, -21.651, tolerance = 0.05)
  expect_equal(tide$lod_tide, 0.248, tolerance = 0.05)

})