diff --git a/applications/solvers/combustion/fireFoam/YEEqn.H b/applications/solvers/combustion/fireFoam/YEEqn.H
index 1c9c2176c5..e968600171 100644
--- a/applications/solvers/combustion/fireFoam/YEEqn.H
+++ b/applications/solvers/combustion/fireFoam/YEEqn.H
@@ -9,7 +9,6 @@ tmp<fv::convectionScheme<scalar>> mvConvection
     )
 );
 {
-    radiation->correct();
     combustion->correct();
     dQ = combustion->dQ();
     volScalarField Yt(0.0*Y[0]);
@@ -48,6 +47,8 @@ tmp<fv::convectionScheme<scalar>> mvConvection
     Y[inertIndex] = scalar(1) - Yt;
     Y[inertIndex].max(0.0);
 
+    radiation->correct();
+
     volScalarField& he = thermo.he();
 
     fvScalarMatrix EEqn
diff --git a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
index 4a976d1b58..91a6b26e2b 100644
--- a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
+++ b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
@@ -68,9 +68,14 @@ Foam::functionObjects::thermoCoupleProbes::thermoCoupleProbes
 
     // Check if the property exist (resume old calculation)
     // or of it is new.
-    if (!getProperty("Tc", Ttc_))
+    dictionary probeDict;
+    if (getDict(typeName, probeDict))
     {
-        Ttc_ = probes::sample(thermo_.T());
+        probeDict.lookup("Tc") >> Ttc_;
+    }
+    else
+    {
+       Ttc_ = probes::sample(thermo_.T());
     }
 
     // Note: can only create the solver once all samples have been found