Premium
A Python Script to Compute Isochrones for MODFLOW
Author(s) -
Feo Alessandra,
Zanini Andrea,
Petrella Emma,
Celico Fulvio
Publication year - 2017
Publication title -
groundwater
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.84
H-Index - 94
eISSN - 1745-6584
pISSN - 0017-467X
DOI - 10.1111/gwat.12588
Subject(s) - modflow , python (programming language) , aquifer , scripting language , computer science , visualization , graphical user interface , groundwater flow , geology , groundwater model , hydrology (agriculture) , computational science , groundwater , programming language , data mining , geotechnical engineering
MODFLOW constitutes today the most popular modeling tool in the study of water flow in aquifers and in modeling aquifers. To simplify the interface to MODFLOW various GUI have been developed for the creation of model definition files and for the visualization and interpretation of results. Recently Bakker et al. (2016) developed the FloPy interface to MODFLOW that allows to import and use the produced simulation data using Python. This allows to construct model input files, run the models, read and plot simulations results through Python scripts. In this note, we present a Python program (that uses FloPy) interface that allows us to generate time‐related capture zones (isochrones) for confined 2D steady‐state groundwater flow in unbounded domains, with one or more wells. As an application, we show a validation of the approach and the results of four basic test cases: a homogenous aquifer with one well, a heterogeneous aquifer with one well, an aquifer with four wells located both longitudinal and perpendicular to the flow direction.