Este ejemplo demuestra las capacidades de PHREEQC para desarrollar simulaciones geoquímicas, con las simulaciones finales dependiendo de resultados de simulaciones previas dentro de la misma corrida. El ejemplo investiga reacciones diagenéticas que pueden ocurrir en zonas donde el agua marina se mezcla con agua subterránea carbonatada. Este ejemplo se ha tomado de la documentación de PHREEQC la cual puede ser encontrada aquí:
https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/html/final-72.html
Contenido
Este ejemplo está dividido en cinco simulaciones:
a) El agua subterránea carbonatada se define equilibrando agua pura con calcita a P CO2 de 10^-2 atm
b) El agua marina se define usando la data de iones principales.
c) Las dos soluciones son mezcladas en proporciones de 70% agua subterránea y 30% agua salina.
d) La mezcla es equilibrada con calcita y dolomita.
e) La mezcla es equilibrada sólo con calcita para investigar la evolución química si se asume que la precipitación de la dolomita es despreciable.

FloPy es la librerìa de Python que construye y ejecuta modelos MODFLOW; esta librería se ha mejorado para proporcionar un soporte completo de MODFLOW 6 y la mayor parte de su desarrollo está relacionado con la funcionalidad de MODFLOW 6, herramientas para usar datos espaciales vectoriales, rasters y funcionalidades comunes de ploteo y exportación.
Python es un lenguaje de programación sencillo y potente; su simplicidad es notable en comparación con otros lenguajes de programación y su potencia se basa en la cantidad de herramientas disponibles para las diferentes áreas de estudio.
Hemos desarrollado un programa completo en modelamiento de aguas subterráneas con MODFLOW y FloPy que abarca desde los temas esenciales en Python hasta casos aplicados de modelamiento regional e intrusión marina. El programa está diseñado para brindarle al estudiante un enfoque metodológico para aprender FloPy basado en una serie de ejemplos aplicados junto con las opiniones/experiencias del profesor y la revisión del software.