Nothing
knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
if(!require("remotes")) { install.packages("remotes") } remotes::install_github(repo = "KWB-R/kwb.hantush", dependencies = TRUE)
library(kwb.hantush)
For comparing the implementation of the Hantush equation, 1967 the results are compared to the USGS benchmark example (page 25). The results can be visualized with the following code:
### Comparision of R results to all ### eight benchmark models in one plot plotModelComparison() ### Comparision of R results to all ### eight benchmark models in multiple plots (one plot for each model) plotModelComparison(layout=c(1,1))
baseProps <- baseProperties( time = 2^(0:6), ## day, for 6 different times ! infiltrationRate = 1, ## meter / day basinWidth = 10, ## meter basinLength = 50, ## meter horizConductivity = 10, ## meter / day iniHead = 10, ## meter specificYield = 0.2)
res <- hantushDistancesBaseProps(baseProps = baseProps)
cols <- length(unique(res$dat[[res$changedBaseProp.Name]])) mainTxt <- sprintf("Changed baseProperty: %s", res$changedBaseProp.Name) xyplot(WLincrease ~ x, groups=res$dat[[res$changedBaseProp.Name]], data=res$dat, type="b", auto.key=list(columns=cols), main=mainTxt)
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.