The use of C eliminates nearly all limitations on array sizes, including numbers of elements, aqueous species, solutions, phases, and lengths of character strings. A new equation solver that optimizes a set of equalities subject to both equality and inequality constraints is used to determine the thermodynamically stable set of phases in equilibrium with a solution. Proper use of the program requires adequate knowledge of geochemistry and a proper formulation of the problem.
Input is arranged in keyword data blocks, which can appear in any order.
Data fields for a keyword are read in a free format, thus they are not column dependent. For 1D advective transport modeling, all the data for reaction-path modeling are needed for each cell in the modeled system. In addition, the identity and composition of a set of plausible reactants and products are needed. K are needed for each additional species and phase.