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()
#> Registered S3 method overwritten by 'xts':
#> method from
#> as.zoo.xts zoo
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
### Comparision of R results to all
### eight benchmark models in multiple plots (one plot for each model)
plotModelComparison(layout=c(1,1))
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rNSE' can not be computed: some elements in 'obs' are zero !
#> Warning: 'rd' can not be computed: some elements in 'obs' are zero !
res <- hantushDistancesBaseProps(baseProps = baseProps)
#> Calculating results for changed baseProperty 'time': 1.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 2.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 4.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 8.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 16.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 32.00000 ... Done!
#> Calculating results for changed baseProperty 'time': 64.00000 ... Done!