Merge pull request #11762 from drjfloyd/master

FDS Source: In geom init Q_RAD_OUT=0 if RADIATION=F
firemodels · Apr 17, 2023 · 9a0b483 · 9a0b483
2 parents 8483b5b + 2e79a64
commit 9a0b483
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/Source/geom.f90 b/Source/geom.f90
@@ -8660,7 +8660,11 @@ SUBROUTINE INIT_CFACE_CELL(NM,ICF,IFACE,CFACE_INDEX,SURF_INDEX,STAGE_FLG,IS_INB,
       CFA%RHO_G = CUT_CELL(ICC)%RHO(JCC)
       B1%ZZ_F(1:N_TOTAL_SCALARS)  = CUT_CELL(ICC)%ZZ(1:N_TOTAL_SCALARS,JCC)
       ! Reinitialize CFACE cell outgoing radiation for change in TMP_F
-      IF (RADIATION) B1%Q_RAD_OUT = SF%EMISSIVITY*SIGMA*B1%TMP_F**4
+      IF (RADIATION) THEN
+         B1%Q_RAD_OUT = SF%EMISSIVITY*SIGMA*B1%TMP_F**4
+      ELSE
+         B1%Q_RAD_OUT = 0._EB
+      ENDIF
       ! Assign normal velocity to CFACE from SURF input:
       B1%U_NORMAL_0 = SF%VEL
       ! Assign normal velocity from VOLUME_FLOW :