I would like to use the labradorite dissolution rate provided in Carroll and Knauss (2005, Chemical Geology) in a kinetic model in GWB. The rate is given by
r = k * 10(-Ea/(2.303RT)) * [(aH+3n/aAl3+n)KT/(1+KT(aH+3n/aAl3+n))]
where the rate constant k is 10-5.69 mol cm-2 s-1, the activation energy Ea is 10.06 kcal mol-1, n is 0.31, and KT is a complex formation constant that varies with temperature. How do I incorporate this rate equation in GWB?
Thanks!