First pass at sorting reactions.

Uses discrepancy function stored in the thermo
object.