Equilibrium phases for more than one solution?

[Jerry]

I want to see the quantities of dissolved/precipitated mineral phases of a solution with varying content of DIC (respectively C(4)). So I can define my solutions with simple SOLUTION keyword or SOLUTION_SPREAD. However, if I calculate the equilibrium phases I always get the output of just one solution (from the last defined one).How can I calculate the equilibrium phases for more than one solution and get these results into the selected output?

[dlparkhurst]

I am afraid you have to set up the calculation for each solution:

SELECTED_OUTPUT
...
END
USE solution 1
USE equilibrium_phases 1
END
USE solution 2
USE equilibrium_phases 1
END
...

Alternatively, you can use IPhreeqc or IPhreeqcCOM to script repetitive calculations.