Function to calculate the derivative of the drawdown with respect to the derivative of log time using the central finite differences

log_derivative_central(t, s, return.pos = T, log = T)

Arguments

t

Numeric vector with the time

s

Numeric vector with the drawdown

return.pos

Logical flag to return only the positive values of the log-derivative (default = TRUE)

log

Logical flag to indicate that natural logarithm (a log to the base e) is used in the derivative calculation (default = TRUE). Logarithm to the base 10 is used if FALSE.

Value

A list with

  • x: Numeric vector with the x coordinates where the log-derivative is evaluated

  • y: Numeric vector with the value of the log-derivative

See also

Examples

data(boulton) t <- boulton$t s <- boulton$s ptest <- pumping_test('Well1', Q = 0.03, r = 20, t = t, s = s) # dptest.cn <- log_derivative_central(ptest$t, ptest$s) plot(t, s, type = "p", log = "xy", ylim = c(1e-3,2))
points(dptest.cn$x, dptest.cn$y, col = "red")