TITLE Hydrothermal dolomiteSOLUTION 1#seawater brine-units mol/kgw-temp 40pH 6.77#anionsAlkalinity 0.000260977Cl 6.15E+00S(6) 2.20E-01#cationsCa 1.10E-02Mg 5.92E-01K 1.15E-01Na 5.28E+00SELECTED_OUTPUT-file Pressure.txt-reset false-alk-solution-temp-pH-alkalinity-totals S(6) Ca K Na Mg Cl #concentrations (moles/kgw)-saturation_indices Calcite Dolomite Anhydrite Halite CO2(g)ENDUSE solution 1EQUILIBRIUM_PHASES#Calcite#DolomiteREACTION_TEMPERATURE250REACTION_PRESSURE1 250 500 750 1000END
TITLE Hydrothermal dolomiteSOLUTION 1#seawater brine-units mol/kgw-temp 40pH 6.77#anionsAlkalinity 0.000260977Cl 6.15E+00S(6) 2.20E-01#cationsCa 1.10E-02Mg 5.92E-01K 1.15E-01Na 5.28E+00ENDUSE solution 1REACTION_TEMPERATURE250SAVE solution 1USER_PRINT10 PUT(MOL("Ca+2"), 1)20 PUT(ACT("Ca+2"), 2)30 PUT(GAMMA("Ca+2"), 3)110 PUT(MOL("CO3-2"), 11)120 PUT(ACT("CO3-2"), 12)130 PUT(GAMMA("CO3-2"), 13)210 PUT(SR("Calcite"), 21)220 PUT(ACT("Ca+2")*ACT("CO3-2"), 22)230 PUT(10^LK_PHASE("Calcite"), 23)ENDPRINT-user_print falseENDUSE solution 1REACTION_PRESSURE1 250 500 750 1000USER_GRAPH 1 -headings pressure SR(Calcite) IAP(Calcite) K(Calcite) -axis_titles "Pressure" "Relative change" -start10 GRAPH_X PRESSURE20 GRAPH_Y SR("Calcite")/GET(21) 30 GRAPH_Y ACT("Ca+2")*ACT("CO3-2") / GET(22)40 GRAPH_Y 10^LK_PHASE("Calcite") / GET(23) -endUSER_GRAPH 2 -headings pressure MOL(Ca+2) ACT(Ca+2) GAMMA(Ca+2) -axis_titles "Pressure" "Relative change" "" -start10 GRAPH_X PRESSURE30 GRAPH_Y MOL("Ca+2")/GET(1), ACT("Ca+2")/GET(2), GAMMA("Ca+2")/GET(3)USER_GRAPH 3 -headings pressure MOL(CO3-2) ACT(CO3-2) GAMMA(CO3-2) -axis_titles "Pressure" "Relative change" "" -start10 GRAPH_X PRESSURE40 GRAPH_Y MOL("CO3-2")/GET(11), ACT("CO3-2")/GET(12), GAMMA("CO3-2")/GET(13) -endEND