GwOptim: The identification of aquifer parameters using inverse...

Description Usage Arguments Value Examples

View source: R/FittingR.R

Description

The identification of aquifer parameters using differential evolution

Usage

1
GwOptim(obs = NA, DimlesTime = NA, OPtFunType = "MSE", ...)

Arguments

obs

observed dimensionless drwandown

DimlesTime

dimensionless time

OPtFunType

type of inverse problem minimization based on MSE - MSE, MAE - MAE, MAx error MAXER

...

additional parameters of RcppDE

Value

list, well pumping test parameters Cbar and Skin, simulated dimensionles drawdown

Examples

1
2
3
4
5
DimlesTime =c (1.88, 3.74, 11.83,47.1,118.31,470.98,1489.37,3741.12, 14893.66)
obsSnormDra = c( 0.09899987,0.19699975,0.60699923,2.18499723,4.58399418,9.32698816,11.05998596,11.64,12.38)
dd=GwOptim(obsSnormDra,DimlesTime,"MSE",lower=c(0,0),upper=c(200000,200),control = DEoptim.control(itermax=100))
plot(log10(DimlesTime), obsSnormDra,pch=19,ylab="Dimensionless Drawdown [-]",xlab = "log10(Dimensionless Time) [-]")
lines(log10(DimlesTime), dd$dimDrawdown,col="red")

petrmaca/GTest documentation built on May 29, 2019, 11:44 a.m.