Skip to main content
U.S. flag

An official website of the United States government

Scripting MODFLOW model development using Python and FloPy

March 30, 2016

Graphical user interfaces (GUIs) are commonly used to construct and postprocess numerical groundwater flow and transport models. Scripting model development with the programming language Python is presented here as an alternative approach. One advantage of Python is that there are many packages available to facilitate the model development process, including packages for plotting, array manipulation, optimization, and data analysis. For MODFLOW-based models, the FloPy package was developed by the authors to construct model input files, run the model, and read and plot simulation results. Use of Python with the available scientific packages and FloPy facilitates data exploration, alternative model evaluations, and model analyses that can be difficult to perform with GUIs. Furthermore, Python scripts are a complete, transparent, and repeatable record of the modeling process. The approach is introduced with a simple FloPy example to create and postprocess a MODFLOW model. A more complicated capture-fraction analysis with a real-world model is presented to demonstrate the types of analyses that can be performed using Python and FloPy.

Publication Year 2016
Title Scripting MODFLOW model development using Python and FloPy
DOI 10.1111/gwat.12413
Authors Mark Bakker, Vincent E. A. Post, Christian D. Langevin, Joseph D. Hughes, Jeremy T. White, Jeffrey Starn, Michael N. Fienen
Publication Type Article
Publication Subtype Journal Article
Series Title Groundwater
Index ID 70170433
Record Source USGS Publications Warehouse
USGS Organization Office of Ground Water