Carbonate-comp Calcite 0-comp Magnesite 0

PHASESMagnesite MgCO3 = Mg+2 + CO3-2 log_k -8.029 delta_h -6.169 kcalENDRATESCalcite-start10 k = PARM(1)20 cal = KIN("Calcite") 30 if cal <= 0 then cal = 1e-640 mag = KIN("Magnesite")50 if mag <= 0 then mag = 1e-660 x = cal/(cal + mag)70 rate = k * (1 - SR("Calcite")/x)80 SAVE rate * time-endMagnesite-start10 k = PARM(1)20 cal = KIN("Calcite") 30 if cal <= 0 then cal = 1e-640 mag = KIN("Magnesite")50 if mag <= 0 then mag = 1e-660 x = mag/(cal + mag)70 rate = k * (1 - SR("Magnesite")/x)80 SAVE rate * time-endSOLUTIONpH 7 chargeCa 10Mg 10C 40 KINETICSCalcite-m 1e-6-parm 1e-8Magnesite-m 1e-6-step 864000 in 10 -parm 1e-8USER_PRINT10 cal = KIN("Calcite")20 mag = KIN("Magnesite")30 xcal = cal/(cal + mag)40 xmag = mag/(cal + mag)50 PRINT "SRcal/xcal: ", SR("Calcite") / xcal60 PRINT "SRmag/xmag: ", SR("Magnesite") / xmagUSER_GRAPH 1 -headings days SIcal SImag Cal Mag -axis_titles "Days" "log10(SatRatio/x)" "Moles" -initial_solutions false -connect_simulations true -plot_concentration_vs x -start10 cal = KIN("Calcite")20 mag = KIN("Magnesite")30 xcal = cal/(cal + mag)40 xmag = mag/(cal + mag)50 GRAPH_X total_time / 8640060 GRAPH_Y LOG10(SR("Calcite") / xcal)70 GRAPH_Y LOG10(SR("Magnesite") / xmag)80 GRAPH_SY cal, mag -end -active true