reactingEulerFoam: Moved fluid.correctThermo() into energy corrector loop

to update the interfacial temperature.
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=2029

applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H

@@ -33,9 +33,10 @@ for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
             EEqn->solve();
         }
     }
+
+    fluid.correctThermo();
 }
 
-fluid.correctThermo();
 
 forAll(phases, phasei)
 {

applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H

@@ -47,9 +47,9 @@ for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
             E2eqn->solve();
         }
     }
-}
 
-fluid.correctThermo();
+    fluid.correctThermo();
+}
 
 Info<< phase1.name() << " min/max T "
     << min(phase1.thermo().T()).value()