El futuro más prometedor de los recursos de aguas subterráneas podría depender en gran medida de las redes sociales y plataformas de difusión masiva, con estrategias de interacción rápida y sensibilización. Sin embargo, para participar en esta nueva tendencia nuestras evaluaciones deben estar en un nivel superior de texto preciso, gráficos agradables y figuras 3D impresionantes.
La visualización avanzada del flujo de agua subterránea y el transporte se considera un punto clave en nuestras evaluaciones. No sabemos si la visualización de las aguas subterráneas puede considerarse una ciencia, aunque dudamos si un artículo científico de gráficos de aguas subterráneas puede ser elegible para una conferencia de hidrogeología o una revista indexada en hidrogeología. Sin embargo, a menos que pueda motivar a las partes interesadas con diferencias finitas o análisis de incertidumbre, se necesitará gráficos potentes (y, por supuesto, precisos) que garanticen la apropiación de conceptos de las características principales del flujo y la calidad del agua subterránea o cualquier tema relacionado que aborde la investigación .
Una nueva versión no oficial de Flopy está disponible como Hataripy. Este "fork" tiene herramientas para la representación de muchas características del modelo, condiciones de borde y datos de salida con opciones para suavizar la cuadrícula. La documentación de esta biblioteca está en progreso, ya que esperamos introducir otras características e identificar algunos errores.
Actualmente la biblioteca puede:
Crear la malla del modelo en la zona activa con valores de punto y celda.
Crear un malla de la capa freática.
Crear la geometría de cualquier condición de contorno.
Instalación
Puede instalar la biblioteca escribiendo esta línea en Anaconda Prompt si está en Windows o en el Shell de Linux:
pip install -i https://test.pypi.org/simple/ hataripy
También puede ver el link del código en Github:
https://github.com/hatarilabs/hataripy
Galeria de fotos
Ejemplo
Para hacer este código, descargue los archivos de este enlace, descomprimirlos en la carpeta "Documentos" y cree un nuevo cuaderno de la carpeta "Scripts". El script ubicará los archivos del modelo y generará los archivos VTU en la dirección especificada. A description of the code can be seen on this video:
Import packages define paths and load model
import os, re, sys, hataripy
import numpy as nphataripy is installed in E:\Software\Anaconda3\lib\site-packages\hataripymodPath = '../Model/'
modName = 'Model1'
exeName = '../Exe/MODFLOW-NWT_64.exe'
mfModel = hataripy.modflow.Modflow.load(modName+'.nam', model_ws=modPath,
exe_name=exeName)# get a list of the model packages
mfModel.get_package_list()['DIS', 'NWT', 'BAS6', 'UPW', 'RCH', 'EVT', 'DRN', 'OC']Define objects that will be represented on the VTKs and add them to geometry object
# read heads from the model output
headArray = hataripy.utils.binaryfile.HeadFile(modPath+modName+'.hds').get_data()
# get information about the drain cells
drnCells = mfModel.drn.stress_period_data[0]# add the arrays to the vtkObject
vtkObject = hataripy.export.vtk.Vtk3D(mfModel,'../vtuFiles/',verbose=True)
vtkObject.add_array('head',headArray)
vtkObject.add_array('drn',drnCells)Create the VTKs for model output, boundary conditions and water table
vtkObject.modelMesh('modelMesh.vtu',smooth=True,cellvalues=['head'])
vtkObject.modelMesh('modelDrn.vtu',smooth=True,cellvalues=['drn'],boundary='drn',avoidpoint=True)
vtkObject.waterTable('waterTable.vtu',smooth=True)Removing existing vtk file: modelMesh.vtu
Writing vtk file: modelMesh.vtu
Number of point is 255920, Number of cells is 31990
Removing existing vtk file: modelDrn.vtu
Writing vtk file: modelDrn.vtu
Number of point is 7432, Number of cells is 929
Removing existing vtk file: waterTable.vtu
Writing vtk file: waterTable.vtu
Number of point is 25592, Number of cells is 6398El procesamiento de los VTK en Paraview se puede ver en este video:
