Merge pull request #13829 from mcgratta/master

FDS Source: Lower C_MIN for RTE correction

Author: mcgratta

Source/cons.f90

@@ -467,9 +467,9 @@ MODULE GLOBAL_CONSTANTS
 REAL(EB) :: RTE_SOURCE_CORRECTION_FACTOR=1._EB   !< Multiplicative factor used in correcting RTE source term
 REAL(EB) :: RAD_Q_SUM=0._EB   !< \f$ \sum_{ijk} \left( \chi_{\rm r} \dot{q}_{ijk}''' + \kappa_{ijk} U_{ijk} \right) V_{ijk} \f$
 REAL(EB) :: KFST4_SUM=0._EB   !< \f$ \sum_{ijk} 4 \kappa_{ijk} \sigma T_{ijk}^4 V_{ijk} \f$
-REAL(EB) :: QR_CLIP           !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
+REAL(EB) :: QR_CLIP=10._EB    !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
 REAL(EB) :: C_MAX=100._EB     !< Maximum value of RAD_Q_SUM/KFST4_SUM
-REAL(EB) :: C_MIN=1._EB       !< Minimum value of RAD_Q_SUM/KFST4_SUM
+REAL(EB) :: C_MIN=0.1_EB      !< Minimum value of RAD_Q_SUM/KFST4_SUM
 
 ! Ramping parameters
 

Source/read.f90

@@ -9851,7 +9851,6 @@ SUBROUTINE READ_RADI
 MIE_MAXIMUM_DIAMETER = 0._EB
 MIE_MINIMUM_DIAMETER = 0._EB
 MIE_NDG              = 50
-QR_CLIP              = 10._EB  ! kW/m3, lower bound for radiation source correction
 OPTICALLY_THIN       = .FALSE.
 
 ! Read radiation parameters
@@ -12519,7 +12518,7 @@ SUBROUTINE READ_INIT
 
    ELSE
 
-      ! Use information from DEVC line to create an INIT line for 'RADIATIVE HEAT FLUX GAS', 'ADIABATIC SURFACE TEMPERATURE GAS', etc.
+      ! Use information from DEVC line to create an INIT for 'RADIATIVE HEAT FLUX GAS', 'ADIABATIC SURFACE TEMPERATURE GAS', etc.
 
       CALL SET_INIT_DEFAULTS
       DV => DEVICE(INIT_RESERVED(N-N_INIT_READ)%DEVC_INDEX)  ! First device in the line of POINTS