Obtaining PHREEQC results in each cell of a 3D PHAST model

[appups]

Hello,

I am interested in reviewing the PHREEQC results in each cell of my 3-D PHAST model (especially the amount of rock dissolution). When you run PHREEQC alone, you can review this information for a 1-D model. 

Is it possible to save the PHREEQC results in each cell for the entire transient simulation? I know it would be a large file but just inquiring.

Any suggestion on how I could know the amount of rock dissolution in each cell would be very helpful.

Thanks for your help.

[dlparkhurst]

The most common way to get information out of PhreeqcRM is through SELECTED_OUTPUT and USER_PUNCH definitions. You must use RunFile or RunString to define these data blocks to the workers, and then enable retrieving the results with SetSelectedOutputOn(true). You can then retrieve and process the selected output data with a number of methods. Look at the documentation SetSelectedOutputOn in the Doxygen/html directory for related methods.

With the latest version of PhreeqcRM (3.8.5 at https://github.com/usgs-coupled/phreeqcrm) you should create a BMIPhreeqcRM object instead of a PhreeqcRM object. The BMIPhreeqcRM object is a subclass of PhreeqcRM so all of the PhreeqcRM methods are available, but additional methods are available for the BMI object. You can retrieve most data GetValues method. Look at the documentation for AddOutputVars for more information.

It is possible to get the PHREEQC output for cells, although as you say, it could be a very large file. However, you can limit the output to selected cells with the method SetPrintChemistryMask. The method SetPrintChemistryOn allows output to be written. The method SetFilePrefix defines the name for the file.

[dlparkhurst]

Sorry, I read PhreeqcRM instead of PHAST. In PHAST, there is a print option for Force chemistry print in the PRINT_FREQUENCY data block. You can restrict that print to specific times in PRINT_FREQUENCY, and you can specify a subset of cells with PRINT_LOCATIONS.