From 8c2bedb23943d9e5ddf469e91a1e02d4a6d01d57 Mon Sep 17 00:00:00 2001
From: mcgratta <mcgratta@gmail.com>
Date: Thu, 5 Mar 2020 14:51:34 -0500
Subject: [PATCH] FDS Source: Properly initialize CHI_R

---
 Source/fire.f90 | 2 ++
 Source/init.f90 | 2 +-
 2 files changed, 3 insertions(+), 1 deletion(-)

diff --git a/Source/fire.f90 b/Source/fire.f90
index 64b250806d7..d861fa6f7ec 100644
--- a/Source/fire.f90
+++ b/Source/fire.f90
@@ -416,6 +416,8 @@ SUBROUTINE COMBUSTION_MODEL(T,DT,ZZ_GET,Q_OUT,MIX_TIME_OUT,CHI_R_OUT,CHEM_SUBIT_
       CHI_R_SUM = CHI_R_SUM + Q_REAC_SUM(NR)*RN%CHI_R*TIME_RAMP_FACTOR
    ENDDO
    CHI_R_OUT = CHI_R_SUM/(SUM(Q_REAC_SUM))
+ELSE
+   CHI_R_OUT = REACTION(1)%CHI_R*EVALUATE_RAMP(T,0._EB,REACTION(1)%RAMP_CHI_R_INDEX)
 ENDIF
 CHI_R_OUT = MAX(CHI_R_MIN,MIN(CHI_R_MAX,CHI_R_OUT))
 
diff --git a/Source/init.f90 b/Source/init.f90
index 2ccbb1967b8..be34ff88566 100644
--- a/Source/init.f90
+++ b/Source/init.f90
@@ -291,7 +291,7 @@ SUBROUTINE INITIALIZE_MESH_VARIABLES_1(DT,NM)
 IF (.NOT.EVACUATION_ONLY(NM)) THEN
    ALLOCATE(M%CHI_R(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO)
    CALL ChkMemErr('INIT','QR',IZERO)
-   M%CHI_R = MAX(0._EB,RADIATIVE_FRACTION)
+   M%CHI_R = MAX(0._EB,RADIATIVE_FRACTION*EVALUATE_RAMP(T_BEGIN,0._EB,REACTION(1)%RAMP_CHI_R_INDEX))
    ALLOCATE(M%QR(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO)
    CALL ChkMemErr('INIT','QR',IZERO)
    M%QR = 0._EB