FDS Source: Add Kp to Thickened Flame Model

Author: cxp484

Source/fire.f90

@@ -886,11 +886,6 @@ SUBROUTINE COMBUSTION_MODEL(T,DT,ZZ_GET,Q_OUT,MIX_TIME_OUT,CHI_R_OUT,CHEM_SUBIT_
          Q_REAC_SUB = 0._EB
    END SELECT INTEGRATOR_SELECT
 
-   ! Consider thickened flame model consideration here
-   IF (COMBUSTION_ODE_SOLVER == CVODE_SOLVER .AND. FLAME_THICK_FACTOR > 1._EB) THEN
-      ZZ_MIXED = ZZ_0 + (ZZ_MIXED-ZZ_0)/FLAME_THICK_FACTOR
-   ENDIF
-
    CALL GET_REALIZABLE_MF(ZZ_MIXED)
    ZZ_GET = ZETA*ZZ_0 + (1._EB-ZETA)*ZZ_MIXED ! FDS Tech Guide (5.19)
    IF (NO_REACTIONS) DT_ITER = DT
@@ -928,6 +923,12 @@ SUBROUTINE COMBUSTION_MODEL(T,DT,ZZ_GET,Q_OUT,MIX_TIME_OUT,CHI_R_OUT,CHEM_SUBIT_
 
 ENDDO INTEGRATION_LOOP
 
+
+! Consider thickened flame model consideration here
+   IF (COMBUSTION_ODE_SOLVER == CVODE_SOLVER .AND. FLAME_THICK_FACTOR > 1._EB) THEN
+      ZZ_GET = ZZ_0 + (ZZ_GET-ZZ_0)/(FLAME_THICK_FACTOR)
+   ENDIF
+   
 ! Compute heat release rate
 
 Q_OUT = -RHO_IN*SUM(SPECIES_MIXTURE%H_F*(ZZ_GET-ZZ_0))/DT ! FDS Tech Guide (5.47)