these Fortran90 source files comprise the unconfined (moving water
table - aka Neuman 1972) slug test solution given in 

Malama, Kuhlman, Barrash, Cardiff & Thoma, 2011. "Modeling slug tests
in unconfined aquifers taking into account water table kinematics,
wellbore skin and inertial effects", Journal of Hydrology

http://dx.doi.org/10.1016/j.hydrol.2011.7.028

The code is available under the MIT open-source license, which allows
free use and redistribution, as long as you acknowledge the original
authors (by citing our original paper, please!)

The code requires a very recent version of the gfortran (>= 4.6)
compiler to properly compile (one that includes libquadmath), but I
have gotten it to work with some minor tweaks using the latest version
of the Intel Fortran compiler (which supports quad precision, but not
for the complex hyperbolic trig functions -- so they must be
implemented in Fortran, rather than in the compiler library).

The code can be compiled, given you have the correct compiler set up 
already, by calling the compiler-specific makefile:

make -f Makefile.gfortran debug_driver
or
make -f Makefile.gfortran driver

these two makefile targets build the debug target (slower, checks for 
bounds overflows, etc) and the standard target (builds using OpenMP 
for a parallel build which is optimized and much faster. 

The input.dat file has all the input parameters read by the program.
Parameters include physical properties of the formation, well geometry
parameters, numerical convergence parameters and the locations/times
to compute the solution at.

Anything to the right is a comment and is ignored by the program when
reading in the input.  Units are arbitrary but consistent.

=== example of input.dat ===

F  F                 :: suppress output to screen?, output dimensionless t?
20.00                :: initial saturated thickness
1.93D+1  1.90D+1     :: l,d;  depth to bottom, top of test interval
2.50D-2  1.30D-2     :: rw,rc;  wellbore, tubing radius
5.924E-4  1.0D0      :: Kr,kappa
6.000E-4  .25000     :: Ss,Sy
9.81D+0  1.20D-6     :: g,nu;  gravitational acceleration, kinematic viscosity at 15C
40  1.0D-8  1.0D-9   :: deHoog,Knight&Stokes inverse Laplace transform algorithm;  M,alpha,tol
7  4                 :: tanh-sinh quadrature;   2^k-1 order, # extrapollation steps
1  1  12  100        :: Gauss-Lobatto quadrature;  min/max zero split, # zeros to integrate, # abcissa/zero
-2  4  400           :: logspace times;  lo, hi, # times to compute
T  26  data.dat      :: compute times?, # times to read, filename for times (1 time/line)
data.out                             :: filename for output

=== end of example ===

If you have any questions or require assistance using the code, please
contact me (Kris Kuhlman) at klkuhlm <at> sandia <dot> gov.  I will
try to help you.  I mostly work in Linux and on a Mac, but I have
gotten codes to work in Windows before.

July, 2011


