Evaporation and Homogeneous Redox Reactions Simulation with Phreeqc
/Evaporation is accomplished by removing water from the chemical system. Water can be removed as an irreversible reactant with a negative reaction coefficient in the REACTION keyword input. This example uses this method; the REACTION data block is used to simulate concentration of rain water by approximately 20 fold by removing 95 percent of the water. The resulting solution contains only about 0.05 kg of water. In a subsequent simulation, the MIX keyword is used to generate a solution that has the same concentrations as the evaporated solution, but has a total of mass of water of approximately 1 kg.

FloPy is the Python library that builds and executes MODFLOW models; this library has been enhanced to provide full support of MODFLOW 6 with most of its recent development is related to functionality for MODFLOW 6, tools to use vector and raster spatial data and common plotting and export functionality.