Transient Stream Drawdown and Depletion

A collection of 5 scripts for computing transient aquifer and stream drawdown induced by groundwater pumping from a well near a stream.

There are 5 scripts, four of which are functions.

1. Prerun_sp.m: Frontend script in which most model parameters are specified as well as the optimization option (using lsqcurvefit) for parameter estimation. This script call the function streamdrawdown(parameters, times);
2. streamdrawdown.m: The script is of the function streamdrawdown(params,times) where params are aquifer parameters, namely Kx=Hydraulic conductivity in x-direction, Ss= Aquifer specific storage, Cr = Stream storage coefficient, \beta = streambed conductance, and \kappa = anisotropy ratio, i.e., \kappa = Ky/Kx. In this script, the fixed variables are specified, namely the pumping rate Q, aquifer thickness b, stream channel width w, (x,y) positions of observation location, and R the perpendicular (shortest) distance from the pumping well to the stream bank. The script calls the function dehoog2sDr(t,x,y,tmax,N,params) to perform the Laplace transform inversion. Note: for stream ; pumping well side of stream is , and other side of stream is .
3. dehoog2sDr.m Inverts Laplace transform to get the time solution. This can be parallelized.
4. inverse_cosFourier.m: performs inversion of the Fourier cosine transform by simple integration relative to the transform variable \xi from 0 to inf to get the solution in y.
5. fourierlaplace_sD.m: This is the solution in the Laplace-Fourier cosine space.