Click here to donate to keep PhreeqcUsers open

100 REM equilibrium_phases ====================== 110 e = SYS("equi", count, equi$, type$, moles) 120 FOR i = 1 TO count 130 min_vol = EQUI(equi$(i)) * PHASE_VM(equi$(i)) 140 min_vol_tot = min_vol_tot + min_vol 150 NEXT i

PHASESFix_H+H+ = H+log_k 0.0SOLUTION 1pH 7.5Alkalinity 3Ca 1 Calcite 0.0 10Na 1 chargeENDUSE solution 1EQUILIBRIUM_PHASES 1Fix_H+ -10 Ca(OH)2 10#CO2(g) -3.4END

PHASESFix_H+H+ = H+log_k 0.0SOLUTION 1pH 7.5Alkalinity 3Ca 1 Calcite 0.0 10Na 1 chargeENDUSE solution 1REACTION 1Ca(OH)2 10.0015 moleEQUILIBRIUM_PHASES 1Fix_H+ -10 CO2 10END

PHASESFix_H+H+ = H+log_k 0.0SOLUTION 1pH 7.5Alkalinity 3Ca 1 Calcite 0.0 10Na 1 chargeEND#USE solution 1#EQUILIBRIUM_PHASES 1#CO2(g) -3.4 10#SAVE solution 1#ENDUSE solution 1EQUILIBRIUM_PHASES 1Fix_H+ -10 Ca(OH)2 10SAVE solution 2ENDUSE solution 2EQUILIBRIUM_PHASES 2Calcite 0 0#CO2(g) -3.4 10END