PROGRAM Anisotropic_Pumping_and_Slug_Test
IMPLICIT NONE
REAL, PARAMETER :: pi = 3.14159265358979323846
REAL :: S, Tx, Ty, phi_rad
INTEGER :: TestType_command
INTEGER :: n
REAL, DIMENSION(:), ALLOCATABLE :: xx, yy, rc, rw, Q, Ho
REAL :: s
REAL :: Tx_est, Ty_est, Phi_est, S_est
REAL, DIMENSION(:), ALLOCATABLE :: x, y
REAL, DIMENSION(:,:), ALLOCATABLE :: time, PAR
REAL, DIMENSION(:,:), ALLOCATABLE :: g, K, gwell, h, hh, H_calc
REAL, DIMENSION(:,:), ALLOCATABLE :: M, Det, tempM, Num, q
REAL, DIMENSION(:), ALLOCATABLE :: A, C, d
REAL, POINTER :: temp(:)
INTEGER :: i, j, kt
REAL :: omegaa
CHARACTER(LEN=100) :: FF
REAL, DIMENSION(:), ALLOCATABLE :: alfa, beta, bdif
! Initial estimates for aquifer parameters
S = 1.0e-7
Tx = 1.0e-6
Ty = 1.0e-3
phi_rad = 0.0
! Test type: 1 = slug test, 2 = pumping test
TestType_command = 1
! Well parameters
n = 2
ALLOCATE(xx(n), yy(n), rc(n), rw(n), Q(n), Ho(n))
xx = [0.0, 2.0]
yy = [0.0, 2.0]
rc = [0.0762, 0.0762]
rw = [0.0762, 0.0762]
Q = [0.03648, 0.0]
Ho = [2.0, 0.0]
IF (TestType_command == 1) THEN
Q = 0.0
ELSE IF (TestType_command == 2) THEN
Ho = 0.0
ENDIF
! Declare symbolic variables
REAL :: Tx_est, Ty_est, Phi_est, S_est