eta 1-(rho_w+rho_is)/rho_w_init "Degree of pyrolysis"
k_t A_t*exp(-E_t/(R_const*T)) "Reaction rate constant for w -> t"
k_is A_is*exp(-E_is/(R_const*T)) "Reaction rate constant for w -> is"
k_g A_g*exp(-E_g/(R_const*T)) "Reaction rate constant for w -> g"
k_c A_c*exp(-E_c/(R_const*T)) "Reaction rate constant for is -> c"
k_c2 A_c2*exp(-E_c2/(R_const*T)) "Reaction rate constant for t -> c"
k_g2 A_g2*exp(-E_g2/(R_const*T)) "Reaction rate constant for t -> g"
Q "Q_w +Q_c + Q_g2 + Q_c2" "Total heat source"
Q_w "-rho_w*(k_t*DH_t + k_g*DH_g +  k_is*DH_is)" "Heat source for primary pyrolysis steps"
Q_c -k_c*rho_is*DH_c "Heat source for char production from intermediate solid"
Q_g2 -k_g2*DH_g2*tcs.rho*w_t "Heat source for gas production from tar"
Q_c2 -k_c2*DH_c2*tcs.rho*w_t "Heat source for char production from tar"