A Python class to export Modflow 2005 model 3D geometry and parameters to Vtk - Tutorial
/Generating 3D visualizations of groundwater models is essential to analyze the flow regime, perform quality checks and see the interaction of the groundwater body with external factors / boundary conditions. The Flopy library has tools to export the parameters, boundary conditions and results that we have modified and compiled within a Python class. The use of this class allows the generation of Vtk files on a friendly way and in few steps. The tutorial also includes a representation of the parameters generated in ParaView.
We have another Python class for Modflow 6 model available on this link.
Tutorial
Code
from workingTools import Mf2005VtkGenerator
# load model for examples
namFile = "ModelWebinar1d.nam"
modelWs = "../Model/"
exeName = "../Bin/"
vtkDir = "../Vtk"
vtkGen = Mf2005VtkGenerator(namFile=namFile,
modelWs=modelWs,
exeName=exeName,
vtkDir=vtkDir)
#load model
vtkGen.loadModel()
Check for files in '../Vtk'.
Deleted file: ../Vtk\modelGeometry.vtk
Deleted file: ../Vtk\modelParameters.vtk
Deleted file: ../Vtk\param_hk.npy
Deleted file: ../Vtk\param_ibound.npy
Deleted file: ../Vtk\param_ss.npy
Deleted file: ../Vtk\param_strt.npy
Deleted file: ../Vtk\param_sy.npy
Deleted file: ../Vtk\param_vka.npy
All files in '../Vtk' have been deleted.
C:\Users\saulm\anaconda3\Lib\site-packages\flopy\mbase.py:97: UserWarning: The program ../Bin/ does not exist or is not executable.
warn(
Models in simulation: ModelWebinar1d
#get package list
vtkGen.getPackageList()
Package list: ['DIS', 'BAS6', 'OC', 'NWT', 'UPW', 'DRN', 'RCH', 'EVT']
#generate model geometry as vtk and parameter array
vtkGen.generateGeometryArrays()
#generate parameter vtk
vtkGen.generateParamVtk()
Parameter Vtk Generated
#not implemented yet!
#generate bc and obs vtk
#vtkGen.generateBcObsVtk()
Input data
You can download the input data on this link.