Nothing
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")
tide <- et$tide()
expect_equal(tide$gravity[1:2],
c(-448.648, -564.549),
tolerance = 0.01
)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.