Merge remote-tracking branch 'firemodels/master' into master2

Author: gforney

Manuals/Bibliography/FDS_general.bib

@@ -4771,7 +4771,7 @@ @ARTICLE{Nakamura:CF2002
 
 @TECHREPORT{NASA_TM_4513,
   author       = {McBride, B. and Gordon, S. and Reno, M.},
-  title        = {{Coefficienes for Calculating Thermodynamic and Transport Properties of Individual Species}},
+  title        = {{Coefficients for Calculating Thermodynamic and Transport Properties of Individual Species}},
   institution  = {National Aeronautics and Space Administration},
   year         = {1993},
   type         = {NASA Technical Memorandum},

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -6374,7 +6374,7 @@ \subsection{Electrical Cables}
 \subsubsection{Target Particles}
 \label{info:target_particles}
 
-To declare solid particle as being just a target for heat transfer, specify \ct{TARGET_ONLY} on \ct{PART}. When this flag is set and material properties are given, then FDS will compute the particle temperature using the local gas temperature and radiative conditions; however, the particle will not exchange energy or momentum with the gas. When this flag is set, a \ct{CARTESIAN} particle's \ct{SURF} will no longer be treated as a half-thickness plate. This allows the use of backside boundary conditions on \ct{SURF} such as \ct{TMP_GAS_BACK} and \ct{EMISSIVITY_BACK}. If it is still desired to have the \ct{THICKNESS} represent half the thickness of a plate heated from both sides, then set \ct{BACKING='INSULTATED'} on the \ct{SURF} line. Target particles are treated as \ct{STATIC} but can be positioned using \ct{PATH_RAMP} and/or \ct{ORIENTATION_RAMP} on \ct{INIT}.
+To declare that a solid particle is only for recording the heat flux to a location away from a solid surface, specify \ct{TARGET_ONLY=T} on the \ct{PART} line. When this flag is set and material properties are given, FDS will compute the particle temperature using the local gas temperature and radiative conditions; however, the particle will not exchange energy or momentum with the gas. When this flag is set, a \ct{CARTESIAN} particle's \ct{SURF} will no longer be treated as a half-thickness plate. This allows the use of backside boundary conditions on \ct{SURF} such as \ct{TMP_GAS_BACK} and \ct{EMISSIVITY_BACK}. If it is still desired to have the \ct{THICKNESS} represent half the thickness of a plate heated from both sides, then set \ct{BACKING='INSULTATED'} on the \ct{SURF} line. Target particles are treated as \ct{STATIC} but can be positioned using \ct{PATH_RAMP} and/or \ct{ORIENTATION_RAMP} on \ct{INIT}.
 
 \newpage
 
@@ -10554,11 +10554,15 @@ \subsection{Heat Flux}
 \end{lstlisting}
 By default, the heat transfer coefficient, $h_{\rm c}$, in Eq.~(\ref{gauge_heat_flux}) is calculated at the solid surface to which the device is attached, based on the specified surface properties and characteristics of the surrounding flow field. However, you may specify a fixed \ct{HEAT_TRANSFER_COEFFICIENT} (W/(m$^2 \cdot$K)) for the gauge on the \ct{PROP} line.
 
-\item \ct{'GAUGE HEAT FLUX GAS'} The same as \ct{'GAUGE HEAT FLUX'}, except that it can be located anywhere within the computational domain and not just at a solid surface. It also has an arbitrary \ct{ORIENTATION} vector that points in any desired direction, much like a heat flux gauge. The \ct{ORIENTATION} vector need not be normalized, as in the following:
+\item \ct{'GAUGE HEAT FLUX GAS'} The same as \ct{'GAUGE HEAT FLUX'}, except that it can be located anywhere within the computational domain and not just at a solid surface. It also has an arbitrary \ct{ORIENTATION} vector that points in any desired direction. The \ct{ORIENTATION} vector need not be normalized, as in the following:
 \begin{lstlisting}
-&DEVC ID='hf', QUANTITY='GAUGE HEAT FLUX GAS', XYZ=..., ORIENTATION=-1,1,0/
+&DEVC ID='hf', QUANTITY='GAUGE HEAT FLUX GAS', XYZ=..., ORIENTATION=-1,1,0, 
+      PROP_ID='my gauge' /
+&PROP ID='my gauge', GAUGE_EMISSIVITY=0.85, HEAT_TRANSFER_COEFFICIENT=15. /
 \end{lstlisting}
-Note that the parameter \ct{SPATIAL_STATISTIC} is not appropriate for this quantity, meaning that you cannot integrate this quantity over a plane or volume. However, you can use the parameter \ct{POINTS} to create a one-dimensional array of these devices (see Sec.~\ref{info:line_file}). Also, the convective component of the heat flux is calculated based on the assumption that the virtual target is a flat plate normal to the \ct{ORIENTATION} vector, and that the heat transfer coefficient is a function of the local gas temperature and velocity which will not be affected by the virtual device. Alternatively, you can specify \ct{HEAT_TRANSFER_COEFFICIENT} on a \ct{PROP} line whose \ct{ID} is specified on the \ct{DEVC} line.
+The \ct{GAUGE_TEMPERATURE}, $T_{\rm gauge}$, \ct{GAUGE_EMISSIVITY}, $\epsilon_{\rm gauge}$, and \ct{HEAT_TRANSFER_COEFFICIENT}, $h_{\rm c}$ (W/(m$^2$~K)), can be specified using a \ct{PROP} line that is referenced by the \ct{DEVC} line. Their default values are \ct{TMPA}, 1, and 10, respectively.
+
+Note that the parameter \ct{SPATIAL_STATISTIC} is not appropriate for this quantity, meaning that you cannot integrate this quantity over a plane or volume. However, you can use the parameter \ct{POINTS} to create a one-dimensional array of these devices (see Sec.~\ref{info:line_file}).
 
 \item \ct{'RADIOMETER'} Similar to a water-cooled heat flux gauge, except that this quantity measures only the net radiative component:
 \be
@@ -11361,30 +11365,30 @@ \section{Gas Phase Output Quantities}
 \ct{CFL}                                       & Section~\ref{info:TIMING}                         &                   & D,I,P,S      \\ \hline
 \ct{CFL MAX}                                   & Section~\ref{info:TIMING}                         &                   & D            \\ \hline
 \ct{CHEMISTRY SUBITERATIONS}                   & Section~\ref{info:chem_integration}               &                   & D,S          \\ \hline
-\ct{CHI_R}                                    & Section~\ref{info:CHI_R}                          &                   & D,I,S        \\ \hline
+\ct{CHI_R}                                     & Section~\ref{info:CHI_R}                          &                   & D,I,S        \\ \hline
 \ct{COMBUSTION EFFICIENCY}                     & $\delta t/\tau_{\mathrm{mix}}$                    &                   & D,I,P,S      \\ \hline
 \ct{CONDUCTIVITY}                              & Section~\ref{info:CONDUCTIVITY}                   & \si{W/(m.K)}      & D,I,P,S      \\ \hline
-\ct{C_SMAG}                                   & Smagorinsky coefficient                           &                   & D,I,P,S      \\ \hline
+\ct{C_SMAG}                                    & Smagorinsky coefficient                           &                   & D,I,P,S      \\ \hline
 \ct{DENSITY}$^1$                               & Total or species density                          & kg/m$^3$          & D,I,P,S      \\ \hline
 \ct{DIFFUSIVE MASS FLUX X}$^1$                 & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline
 \ct{DIFFUSIVE MASS FLUX Y}$^1$                 & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline
 \ct{DIFFUSIVE MASS FLUX Z}$^1$                 & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline
 \ct{DIFFUSIVITY}$^1$                           & Section~\ref{info:diffusivity}                    & m$^2$/s           & D,I,P,S      \\ \hline
-\ct{DISSIPATION RATE}                          & $\mu / \rho \times $ \ct{STRAIN RATE}            & m$^2$/s$^3$       & D,I,P,S      \\ \hline
+\ct{DISSIPATION RATE}                          & $(\mu / \rho)\, S_{ij}S_{ij}$                     & m$^2$/s$^3$       & D,I,P,S      \\ \hline
 \ct{DIVERGENCE}                                & $\nabla \cdot \bu$                                & 1/s               & D,I,P,S      \\ \hline
 \ct{EFFECTIVE FLAME TEMPERATURE}               & Section~\ref{info:Flame_Temperature}              & $^\circ$C         & D,I,P,S      \\ \hline
 \ct{ENTHALPY}                                  & Section~\ref{info:Enthalpy}                       & kJ/m$^3$          & D,I,P,S      \\ \hline
-\ct{ENTHALPY FLUX X}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}          & D,I,P,S      \\ \hline
-\ct{ENTHALPY FLUX Y}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}          & D,I,P,S      \\ \hline
-\ct{ENTHALPY FLUX Z}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}          & D,I,P,S      \\ \hline
+\ct{ENTHALPY FLUX X}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}     & D,I,P,S      \\ \hline
+\ct{ENTHALPY FLUX Y}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}     & D,I,P,S      \\ \hline
+\ct{ENTHALPY FLUX Z}                           & Section~\ref{info:enthalpy_flux}                  & \unit{kW/m^2}     & D,I,P,S      \\ \hline
 \ct{EXTINCTION}                                & Section~\ref{info:extinct_out}                    &                   & D,S          \\ \hline
 \ct{EXTINCTION COEFFICIENT}                    & Section~\ref{info:visibility}                     & 1/m               & D,I,P,S      \\ \hline
-\ct{F_X, F_Y, F_Z}                          & Momentum terms, $F_x$, $F_y$, $F_z$               & m/s$^2$           & D,I,P,S      \\ \hline
+\ct{F_X, F_Y, F_Z}                             & Momentum terms, $F_x$, $F_y$, $F_z$               & m/s$^2$           & D,I,P,S      \\ \hline
 \ct{H}                                         & $\cH=|\bu|^2/2 + \tp/\rho    $                    & (m/s)$^2$         & D,I,P,S      \\ \hline
 \ct{HRRPUV}                                    & $\dq'''$                                          & kW/m$^3$          & D,I,P,S      \\ \hline
-\ct{HRRPUV REAC}$^6$                           & $\dq'''$ for \ct{REAC_ID}                       & kW/m$^3$          & D,S          \\ \hline
+\ct{HRRPUV REAC}$^6$                           & $\dq'''$ for \ct{REAC_ID}                         & kW/m$^3$          & D,S          \\ \hline
 \ct{IDEAL GAS PRESSURE}                        & $\bar{p}=\rho R T / \overline{W}$                 & Pa                & D,I,P,S      \\ \hline
-\ct{INTEGRATED INTENSITY}                      & $U=\int I \, \d \bs$                              & \unit{kW/m^2}          & D,I,P,S      \\ \hline
+\ct{INTEGRATED INTENSITY}                      & $U=\int I \, \d \bs$                              & \unit{kW/m^2}     & D,I,P,S      \\ \hline
 \ct{INTERNAL ENERGY}                           & $\rho h - \bar{p}$                                & kJ/m$^3$          & D,I,P,S      \\ \hline
 \ct{KINETIC ENERGY}                            & Staggered $(u^2+v^2+w^2)/2$                       & (m/s)$^2$         & D,I,P,S      \\ \hline
 \ct{KOLMOGOROV LENGTH SCALE}                   & Section \ref{info:meshquality}                    & m                 & D,I,P,S      \\ \hline
@@ -11405,8 +11409,10 @@ \section{Gas Phase Output Quantities}
 \ct{PRESSURE ZONE}                             & Section~\ref{info:ZONE}                           &                   & D,S          \\ \hline
 \ct{Q CRITERION}                               & $\frac{1}{2} [\mathrm{trace}(\nabla \mathbf{u})^2 - \mathrm{trace}( (\nabla \mathbf{u})^2)]$ & 1/s$^2$ & D,I,P,S \\ \hline
 \ct{RADIAL VELOCITY}                           & $(u,v)\cdot(x,y)/\sqrt{x^2+y^2}$                  & m/s               & D,I,P,S      \\ \hline
-\ct{RADIATION LOSS}                            & $\nabla \cdot \bq_r''$                            & kW/m$^3$          & D,I,P,S      \\ \hline
-\ct{RADIATIVE HEAT FLUX GAS}                   & Section~\ref{info:heat_flux}                      & \unit{kW/m^2}          & D            \\ \hline
+\ct{RADIATION ABSORPTION}                      & $\kappa U$                                        & kW/m$^3$          & D,I,P,S      \\ \hline
+\ct{RADIATION EMISSION}                        & $4 \kappa \sigma T^4$                             & kW/m$^3$          & D,I,P,S      \\ \hline
+\ct{RADIATION LOSS}                            & $\nabla \cdot \bq_r''=\kappa (U-4 \sigma T^4)$    & kW/m$^3$          & D,I,P,S      \\ \hline
+\ct{RADIATIVE HEAT FLUX GAS}                   & Section~\ref{info:heat_flux}                      & \unit{kW/m^2}     & D            \\ \hline
 \ct{REAC SOURCE TERM}$^1$                      & $\dot{m}_\alpha^{\prime\prime\prime}$             & kg/m$^3$          & D,I,P,S      \\ \hline
 \ct{RELATIVE HUMIDITY}                         & Section~\ref{info:simple_chemistry}               & \%                & D,I,P,S      \\ \hline
 \ct{RESOLVED KINETIC ENERGY}                   & $k_{res} = (\bar{u}^2+\bar{v}^2+\bar{w}^2)/2$     & (m/s)$^2$         & D,I,P,S      \\ \hline
@@ -11423,7 +11429,7 @@ \section{Gas Phase Output Quantities}
 \ct{SUBGRID KINETIC ENERGY}                    & Section~\ref{info:meshquality}                    & \si{m^2/s^2}      & D,S          \\ \hline
 \ct{SUM LUMPED MASS FRACTIONS}                 & $\sum_i Z_i$ (should be 1)                        &                   & D,S          \\ \hline
 \ct{SUM PRIMITIVE MASS FRACTIONS}              & $\sum_\alpha Y_\alpha$ (should be 1)              &                   & D,S          \\ \hline
-\ct{TEMPERATURE}                               & Section~\ref{info:Flame_Temperature}                   & $^\circ$C         & D,I,P,S      \\ \hline
+\ct{TEMPERATURE}                               & Section~\ref{info:Flame_Temperature}              & $^\circ$C         & D,I,P,S      \\ \hline
 \ct{TOTAL MASS FLUX X}$^1$                     & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline
 \ct{TOTAL MASS FLUX Y}$^1$                     & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline
 \ct{TOTAL MASS FLUX Z}$^1$                     & Section~\ref{info:mass_flow}                      & kg/s/m$^2$        & D,I,P,S      \\ \hline

Source/cons.f90

@@ -256,6 +256,7 @@ MODULE GLOBAL_CONSTANTS
 LOGICAL :: POSITIVE_ERROR_TEST=.FALSE.
 LOGICAL :: OBST_SHAPE_AREA_ADJUST=.FALSE.
 LOGICAL :: STORE_SPECIES_FLUX=.FALSE.
+LOGICAL :: STORE_RADIATION_TERMS=.FALSE.
 LOGICAL :: STORE_PRESSURE_POISSON_RESIDUAL=.FALSE.
 LOGICAL :: OXIDATION_REACTION=.FALSE.
 LOGICAL :: PERIODIC_DOMAIN_X=.FALSE.                !< The domain is periodic \f$ x \f$

Source/data.f90

@@ -439,6 +439,16 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(77)%UNITS = ' '
 OUTPUT_QUANTITY(77)%SHORT_NAME = 'phi'
 
+! Special radiation terms
+
+OUTPUT_QUANTITY(78)%NAME = 'RADIATION EMISSION'
+OUTPUT_QUANTITY(78)%UNITS = 'kW/m3'
+OUTPUT_QUANTITY(78)%SHORT_NAME = 'emit'
+
+OUTPUT_QUANTITY(79)%NAME = 'RADIATION ABSORPTION'
+OUTPUT_QUANTITY(79)%UNITS = 'kW/m3'
+OUTPUT_QUANTITY(79)%SHORT_NAME = 'abs'
+
 ! Cell indices (for Smokeview debugging)
 
 OUTPUT_QUANTITY(80)%NAME = 'CELL INDEX I'

Source/dump.f90

@@ -7618,6 +7618,10 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
 
    CASE(77)  ! LEVEL SET VALUE
       GAS_PHASE_OUTPUT_RES = PHI_LS(II,JJ)
+   CASE(78)  ! RADIATION EMISSION
+      GAS_PHASE_OUTPUT_RES = RADIATION_EMISSION(II,JJ,KK)*0.001_EB
+   CASE(79)  ! RADIATION ABSORPTION
+      GAS_PHASE_OUTPUT_RES = RADIATION_ABSORPTION(II,JJ,KK)*0.001_EB
    CASE(80)  ! CELL INDEX I
       GAS_PHASE_OUTPUT_RES = REAL(II,EB)
    CASE(81)  ! CELL INDEX J

Source/init.f90

@@ -291,6 +291,13 @@ SUBROUTINE INITIALIZE_MESH_VARIABLES_1(DT,NM)
 ALLOCATE(M%KAPPA_GAS(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO) ; CALL ChkMemErr('INIT','KAPPA_GAS',IZERO) ; M%KAPPA_GAS = 0._EB
 ALLOCATE(M%UII(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO)       ; CALL ChkMemErr('INIT','UII',IZERO)       ; M%UII       = 0._EB
 
+IF (STORE_RADIATION_TERMS) THEN
+   ALLOCATE(M%RADIATION_EMISSION(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO) ; CALL ChkMemErr('INIT','RADIATION_EMISSION',IZERO)
+   M%RADIATION_EMISSION = 0._EB
+   ALLOCATE(M%RADIATION_ABSORPTION(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO) ; CALL ChkMemErr('INIT','RADIATION_ABSORPTION',IZERO)
+   M%RADIATION_ABSORPTION = 0._EB
+ENDIF
+
 ! Work arrays
 
 ALLOCATE(M%WORK1(0:IBP1,0:JBP1,0:KBP1),STAT=IZERO) ; CALL ChkMemErr('INIT','WORK1',IZERO)

Source/mesh.f90

@@ -46,6 +46,8 @@ MODULE MESH_VARIABLES
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: CHI_R   !< Radiative fraction, \f$ \chi_{{\rm r},ijk} \f$
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: QR      !< Radiation source term, \f$ -\nabla \cdot \dot{\mathbf{q}}_{\rm r}'' \f$
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: QR_W    !< Radiation source term, particles and droplets
+   REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: RADIATION_EMISSION   !< Radiation source term, emission
+   REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: RADIATION_ABSORPTION !< Radiation source term, absorption
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: UII     !< Integrated intensity, \f$ U_{ijk}=\sum_{l=1}^N I_{ijk}^l\delta\Omega^l\f$
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: RSUM    !< \f$ R_0 \sum_\alpha Z_{\alpha,ijk}/W_\alpha \f$
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: D_SOURCE!< Source terms in the expression for the divergence
@@ -348,7 +350,7 @@ MODULE MESH_POINTERS
 
 REAL(EB), POINTER, DIMENSION(:,:,:) :: &
    U,V,W,US,VS,WS,DDDT,D,DS,H,HS,KRES,FVX,FVY,FVZ,FVX_B,FVY_B,FVZ_B,FVX_D,FVY_D,FVZ_D,RHO,RHOS, &
-   MU,MU_DNS,TMP,Q,KAPPA_GAS,CHI_R,QR,QR_W,UII,RSUM,D_SOURCE, &
+   MU,MU_DNS,TMP,Q,KAPPA_GAS,CHI_R,QR,QR_W,RADIATION_EMISSION,RADIATION_ABSORPTION,UII,RSUM,D_SOURCE, &
    CSD2,MTR,MSR,WEM,MIX_TIME,CHEM_SUBIT,STRAIN_RATE,D_Z_MAX,PP_RESIDUAL,LES_FILTER_WIDTH
 REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZ,ZZS,REAC_SOURCE_TERM,DEL_RHO_D_DEL_Z,FX,FY,FZ, &
                                          SWORK1,SWORK2,SWORK3,SWORK4, &
@@ -525,6 +527,8 @@ SUBROUTINE POINT_TO_MESH(NM)
 CHI_R => M%CHI_R
 QR=>M%QR
 QR_W=>M%QR_W
+RADIATION_EMISSION=>M%RADIATION_EMISSION
+RADIATION_ABSORPTION=>M%RADIATION_ABSORPTION
 KAPPA_GAS=>M%KAPPA_GAS
 UII=>M%UII
 PP_RESIDUAL=>M%PP_RESIDUAL

Source/radi.f90

@@ -3546,6 +3546,10 @@ SUBROUTINE RADIATION_FVM
 
 IF (WIDE_BAND_MODEL .OR. WSGG_MODEL) THEN
    QR = 0._EB
+   IF (STORE_RADIATION_TERMS) THEN
+      RADIATION_EMISSION   = 0._EB
+      RADIATION_ABSORPTION = 0._EB
+   ENDIF
    IF (N_LP_ARRAY_INDICES>0 .AND. PARTICLES_EXISTED) QR_W = 0._EB
 ENDIF
 
@@ -4476,6 +4480,10 @@ SUBROUTINE RADIATION_FVM
 
    IF (WIDE_BAND_MODEL .OR. WSGG_MODEL) THEN
       QR = QR + KAPPA_GAS*UIID(:,:,:,IBND)-KFST4_GAS
+      IF (STORE_RADIATION_TERMS) THEN
+         RADIATION_EMISSION   = RADIATION_EMISSION + KFST4_GAS
+         RADIATION_ABSORPTION = RADIATION_ABSORPTION + KAPPA_GAS*UIID(:,:,:,IBND)
+      ENDIF
       IF (N_LP_ARRAY_INDICES>0 .AND. PARTICLES_EXISTED) QR_W = QR_W + KAPPA_PART*UIID(:,:,:,IBND) - KFST4_PART
    ENDIF
 
@@ -4506,6 +4514,10 @@ SUBROUTINE RADIATION_FVM
 IF (.NOT. (WIDE_BAND_MODEL .OR. WSGG_MODEL)) THEN
    QR = KAPPA_GAS*UII - KFST4_GAS
    IF (N_LP_ARRAY_INDICES>0) QR_W = QR_W + KAPPA_PART*UII - KFST4_PART
+   IF (STORE_RADIATION_TERMS) THEN
+      RADIATION_EMISSION   = KFST4_GAS
+      RADIATION_ABSORPTION = KAPPA_GAS*UII
+   ENDIF
 ENDIF
 
 ! Calculate the incoming radiative flux onto the solid particles

Source/read.f90

@@ -15966,6 +15966,7 @@ SUBROUTINE GET_QUANTITY_INDEX(SMOKEVIEW_LABEL,SMOKEVIEW_BAR_LABEL,OUTPUT_INDEX,O
                               Y_INDEX,Z_INDEX,PART_INDEX,DUCT_INDEX,NODE_INDEX,REAC_INDEX,MATL_INDEX,OUTTYPE, &
                               QUANTITY,QUANTITY2,SPEC_ID_IN,PART_ID,DUCT_ID,NODE_ID,REAC_ID,MATL_ID,CELL_L,&
                               DUCT_CELL_INDEX,SLICETYPE)
+
 CHARACTER(*), INTENT(INOUT) :: QUANTITY
 CHARACTER(*), INTENT(OUT) :: SMOKEVIEW_LABEL,SMOKEVIEW_BAR_LABEL
 CHARACTER(*) :: SPEC_ID_IN,PART_ID,DUCT_ID,NODE_ID
@@ -16165,6 +16166,8 @@ SUBROUTINE GET_QUANTITY_INDEX(SMOKEVIEW_LABEL,SMOKEVIEW_BAR_LABEL,OUTPUT_INDEX,O
     QUANTITY=='TOTAL MASS FLUX X'     .OR. QUANTITY=='TOTAL MASS FLUX Y'     .OR.  QUANTITY=='TOTAL MASS FLUX Z'     .OR. &
     QUANTITY=='TOTAL MASS FLUX WALL') STORE_SPECIES_FLUX = .TRUE.
 
+IF (QUANTITY=='RADIATION EMISSION' .OR. QUANTITY=='RADIATION ABSORPTION') STORE_RADIATION_TERMS = .TRUE.
+
 IF (TRIM(QUANTITY)=='HRRPUV REAC') THEN
    IF (TRIM(REAC_ID)=='null') THEN
       WRITE(MESSAGE,'(A)') 'ERROR(1011): HRRPUV REAC requires a REAC_ID.'