Stress and temperature relaxation towards equilibrium, hotspot burn with separate rate: "separate_surface_reactions"

Reaction model

As for `surface_reactions', except that the surface reactions have their own rate law rather than using the bulk rate law.

Equilibration method

As for `surface_reactions'.

Input

As for `surface_reactions', with in addition:
For each reaction:
   surface_reaction_rate (see rate laws