RATES Qtz-start10 if TOTAL_TIME > 60*60*4 then goto 10020 rate = 2e-5/3600 # mol/s30 goto 200100 REM after 4 hours110 rate = 1e-5/3600 # mol/s200 moles = rate * TIME210 SAVE moles-endENDSOLUTIONENDINCREMENTAL_REACTIONSKINETICSQtz-formula SiO2 1-step 86400 in 48 stepsENDUSER_GRAPH 1 -headings Hours Si -axis_titles "Time, in hours" "Si, mmol/kgw" "" -initial_solutions false -connect_simulations true -plot_concentration_vs x -start10 graph_x TOTAL_TIME/360020 graph_y TOT("Si")*1000 -end -active trueUSE solution 1USE kinetics 1END