RATES Dolomite -start 2 if tot("Mg") > 1e-12 then f_dol = log(SR("Dolomite")) else f_dol = -99 10 rate = -1.2e-10 * f_dol * (m/m0)^2 # A0/V = 1/cm, decreases with m/m0 20 save rate * time -end END SOLUTION 1-10 -temp 10; pH 7 charge; Ca 1 Calcite; C(4) 1 CO2(g) -2.5 EQUILIBRIUM_PHASES 1-20 CO2(g) -2.5 ; Calcite 0 0.16 # m0 = field value / 200 END KINETICS 1-20 Dolomite; -m0 0.04 # m0 = field value / 200 END SOLUTION 0 -temp 10; pH 7 charge; C(4) 1 CO2(g) -2.5 END PRINT # Going to calculate transport now, please wait... -reset false; -status false TRANSPORT; -cells 10; -len 0.1; -time_step 3.1536e6 # 0.1 yr -diffc 0; -war false; -shift 600; -punch_frequency 300 # graph each 30 yr USER_GRAPH -init false -axis_titles "Depth / m" "Ca, Mg / (mmol/L)" "dol, cc / wt%" -head dist Ca Mg dol cc -connect_simulations false -start 10 graph_x dist 20 graph_y tot("Ca")*1e3, tot("Mg")*1e3 25 f = 100 * 0.1 / 1.8 # % * eps_w / rho_b 30 graph_sy kin("Dolomite") * 0.186 * f, equi("Calcite") * 0.1 * f -end SELECTED_OUTPUT -file 5_15.prn; -reset false; -time true; -dist true; -tot Ca Mg -kin Dolomite; -equi Calcite END