Numerical analysis of confined-unconfined aquifers