k1 A1*exp(-E1/Rg/mod1.T) Rate constant k2 A2*exp(-E2/Rg/mod1.T) Rate constant r1 k1*mod1.c_A Reaction rate r2 k2*mod1.c_B Reaction rate Q_rxn -r1*H1-r2*H2 heat source of reaction