The R package kwb.wtaq provides a programmatic interface to the well drawdown model WTAQ version 2.1 (developed by the U.S. Geological Survey). The WTAQ 2.1 model engine is included in the package, but also available for download as self-extracting win-zip file from the USGS .
Launch tutorial in cloud RStudio session (and install kwb.wtaq from master
branch):
For installing the latest release of this R package run the following code below:
# Enable repository from kwb-r
options(repos = c(
kwbr = 'https://kwb-r.r-universe.dev',
CRAN = 'https://cloud.r-project.org'))
# Download and install kwb.wtaq in R
install.packages('kwb.wtaq')
# Browse the kwb.wtaq manual pages
help(package = 'kwb.wtaq')
Compilation of the WTAQ-2.1 source code (/src) throws the following messages:
-
"Error: Arithmetic overflow converting INTEGER(16) to INTEGER(4) at (1)" (fixed in /src/Makevars with PKG_FFLAGS = -fno-range-check)
-
"Warning: Named COMMON block 'par6' at (1) shall be of the same size" (unfixed help in solving this is highly appreciated!)
However, at least these issues do not cause calculation problems for the test carried out under windows:
### Load/install required packages
if(!require("testthat")) { install.packages("testthat") }
library(kwb.wtaq)
library(testthat)
### Download test file from Github
test_url <- "https://raw.githubusercontent.com/KWB-R/kwb.wtaq/master/tests/testthat/test_windowsCompilation.R"
test_path <- file.path(tempdir(),"test_windowsCompilation.R")
download.file(test_url,destfile = test_path)
### Run test (only working on windows!)
testthat::test_file(test_path)