tests/testthat/test-et_predict.R

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

context("test-et_predict")

test_that("et_predict works", {
  
  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)
  et$predict(method = 'gravity')
  tide <- et$tide()
  
  expect_equal(tide$gravity, 
               c(-448.580, -564.521), 
               tolerance = .0001) 
  
  et <- Earthtide$new(utc = tms, 
                      latitude = 52.3868,
                      longitude = 9.7144,
                      elevation = 110,
                      gravity = 9.8127, 
                      cutoff = 1.0e-10,
                      wave_groups = wave_groups)
  et$predict(method = 'gravity')
  tide <- et$tide()
  
  expect_equal(tide$gravity, 
               c(-448.648, -564.549), 
               tolerance = .0001)
  
  
  
  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)
  
  et$predict(method = 'gravity', astro_update = 24)
  tide <- et$tide()
  
  expect_equal(tide$gravity[1:2], 
               c(-448.648, -564.549), 
               tolerance = 0.01)
  
  
  
})