diff --git a/Source/fire.f90 b/Source/fire.f90
index 36052b18e5..05538f99ae 100644
--- a/Source/fire.f90
+++ b/Source/fire.f90
@@ -1422,10 +1422,11 @@ SUBROUTINE EXTINCT_2(EXTINCT,ZZ_0,ZZ_IN,TMP_IN,TAU_MIX)
 USE PHYSICAL_FUNCTIONS, ONLY: GET_ENTHALPY
 REAL(EB),INTENT(IN) :: TMP_IN,ZZ_IN(1:N_TRACKED_SPECIES),ZZ_0(1:N_TRACKED_SPECIES),TAU_MIX
 LOGICAL, INTENT(INOUT) :: EXTINCT
-REAL(EB) :: ZZ_HAT_0(1:N_TRACKED_SPECIES),ZZ_HAT(1:N_TRACKED_SPECIES),H_0,H_CRIT,PHI_TILDE,CFT
+REAL(EB) :: ZZ_HAT_0(1:N_TRACKED_SPECIES),ZZ_HAT(1:N_TRACKED_SPECIES),H_0,H_CRIT,PHI_TILDE,CFT,TMP_INF
 INTEGER :: NS,NR
 REAL(EB) :: SUM_ZZ,SUM_CFT
 TYPE(REACTION_TYPE), POINTER :: RN,R1
+REAL(EB), PARAMETER :: DELTA_TMP_INF=-425._EB, TMP_EXT_0 = 1700._EB, TAU_CHEM_LOC = 1.45E-3_EB
 
 ! Get the weighted average of the critical flame temperature (CFT) based on the relative amounts of fuels of the primary reactions
 
@@ -1444,7 +1445,8 @@ SUBROUTINE EXTINCT_2(EXTINCT,ZZ_0,ZZ_IN,TMP_IN,TAU_MIX)
 ENDIF
 
 IF (VARIABLE_CFT) THEN
-   CFT = 1522.1_EB*(1000._EB*TAU_MIX)**(-0.038_EB)  ! Experimental mechanism, U of Ghent
+   TMP_INF = SUM_CFT/SUM_ZZ + DELTA_TMP_INF
+   CFT = TMP_INF + (TMP_EXT_0-TMP_INF)*EXP(-SQRT(TAU_MIX/TAU_CHEM_LOC)) ! Experimental mechanism, UGent
 ELSE
    CFT = SUM_CFT/SUM_ZZ
 ENDIF