tests/testthat/test-fit_tidal_hsieh_1987.R

In jkennel/aquifer: Well solutions for aquifer tests

test_that(" works", {

  storativity <- 1e-07
  transmissivity <- 1e-05
  radius_well <- 0.1
  frequency <- 1.9322736/86400
  hsieh <- tidal_hsieh_1987(frequency, storativity, transmissivity, radius_well)

  res <- fit_tidal_hsieh_1987(frequency = 1.9322736/86400,
                              Mod(hsieh$response / storativity),       # m / nstrain
                              Arg(hsieh$response/ storativity),       # degrees
                              radius_well,
                              radius_casing = radius_well)

  expect_equivalent(res$par, log10(c(storativity, transmissivity)), tolerance = 0.05)
})