FDS Source: add SKIP_INRAD to SURF

Author: rmcdermo

Source/radi.f90

@@ -4443,6 +4443,7 @@ SUBROUTINE RADIATION_FVM
          IF (WC%BOUNDARY_TYPE/=SOLID_BOUNDARY) CYCLE
          B1 => BOUNDARY_PROP1(WC%B1_INDEX)
          SF => SURFACE(WC%SURF_INDEX)
+         IF (SF%SKIP_INRAD) INRAD_W(IW) = 0._EB
          IF (SF%EXTERNAL_FLUX > TWO_EPSILON_EB) THEN
             IF (ABS(T_BEGIN) <= SPACING(B1%T_IGN)) THEN
                TSI = T
@@ -4461,6 +4462,7 @@ SUBROUTINE RADIATION_FVM
          IF (CFACE(ICF)%BOUNDARY_TYPE==NULL_BOUNDARY) CYCLE
          B1 => BOUNDARY_PROP1(CFA%B1_INDEX)
          SF => SURFACE(CFA%SURF_INDEX)
+         IF (SF%SKIP_INRAD) INRAD_W(ICF) = 0._EB
          IF (SF%EXTERNAL_FLUX > TWO_EPSILON_EB) THEN
             IF (ABS(T_BEGIN) <= SPACING(B1%T_IGN)) THEN
                TSI = T

Source/read.f90

@@ -7907,7 +7907,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
             VEG_LSET_WIND_EXP,VEG_LSET_BETA,VEG_LSET_HT,VEG_LSET_SIGMA,VEG_LSET_ROS_00, &
             VEG_LSET_M1,VEG_LSET_M10,VEG_LSET_M100,VEG_LSET_MLW,VEG_LSET_MLH,VEG_LSET_SURF_LOAD,VEG_LSET_FIREBASE_TIME,&
             VEG_LSET_CHAR_FRACTION,VEG_LSET_WIND_HEIGHT,VEL_PART,INIT_PER_AREA,TIME_STEP_FACTOR
-LOGICAL :: DEFAULT,VEG_LSET_SPREAD,VEG_LSET_TAN2,TGA_ANALYSIS,VARIABLE_THICKNESS,HT3D,THERM_THICK,VEG_LSET_ROS_FIXED
+LOGICAL :: DEFAULT,VEG_LSET_SPREAD,VEG_LSET_TAN2,TGA_ANALYSIS,VARIABLE_THICKNESS,HT3D,THERM_THICK,VEG_LSET_ROS_FIXED,SKIP_INRAD
 LOGICAL, ALLOCATABLE, DIMENSION(:) :: DUPLICATE
 ! Ember generating variables
 REAL(EB) :: EMBER_GENERATION_HEIGHT(2),EMBER_IGNITION_POWER_MEAN,EMBER_IGNITION_POWER_SIGMA,EMBER_TRACKING_RATIO,EMBER_YIELD
@@ -7936,7 +7936,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
                 RAMP_V,RAMP_V_X,RAMP_V_Y,RAMP_V_Z,&
                 REFERENCE_HEAT_FLUX,REFERENCE_HEAT_FLUX_TIME_INTERVAL,MINIMUM_SCALING_HEAT_FLUX,MAXIMUM_SCALING_HEAT_FLUX,&
                 REFERENCE_THICKNESS,REMESH_RATIO,&
-                RGB,ROUGHNESS,SHAPE_FACTOR,SPEC_ID,&
+                RGB,ROUGHNESS,SHAPE_FACTOR,SKIP_INRAD,SPEC_ID,&
                 SPREAD_RATE,STRETCH_FACTOR,SURFACE_VOLUME_RATIO,&
                 TAU_EF,TAU_MF,TAU_PART,TAU_Q,TAU_T,TAU_V,TEXTURE_HEIGHT,TEXTURE_MAP,TEXTURE_WIDTH,&
                 TGA_ANALYSIS,TGA_CONVERSION_FACTOR,TGA_DT,TGA_DUMP,TGA_FINAL_TEMPERATURE,TGA_HEATING_RATE,THICKNESS,&
@@ -8463,6 +8463,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
    SF%NO_SLIP              = NO_SLIP
    SF%FYI                  = FYI
    SF%EXTERNAL_FLUX        = 1000._EB*EXTERNAL_FLUX
+   SF%SKIP_INRAD           = SKIP_INRAD
    SF%HORIZONTAL           = HORIZONTAL
    SF%INIT_IDS             = INIT_IDS
    IF (INIT_IDS(1)/='null') INIT_INVOKED_BY_SURF = .TRUE.
@@ -9152,6 +9153,7 @@ SUBROUTINE SET_SURF_DEFAULTS
    ROUGHNESS = -1._EB !4.5E-5_EB  ! meters, commercial steel
 ENDIF
 SHAPE_FACTOR            = 0.25_EB
+SKIP_INRAD              = .FALSE.
 SPEC_ID                 = 'null'
 SPREAD_RATE             = -1._EB
 STRETCH_FACTOR          = 2._EB

Source/type.f90

@@ -998,6 +998,7 @@ MODULE TYPES
    LOGICAL :: INERT_Q_REF                            !< Treat REFERENCE_HEAT_FLUX as an inert atmosphere test
    LOGICAL :: ALLOW_UNDERSIDE_PARTICLES=.FALSE.      !< Allow droplets to move along downward facing surfaces
    LOGICAL :: ALLOW_SURFACE_PARTICLES=.TRUE.         !< Allow particles to live on a solid surface
+   LOGICAL :: SKIP_INRAD = .FALSE.                   !< Only apply external flux to a surface
    INTEGER :: GEOMETRY,BACKING,PROFILE,HEAT_TRANSFER_MODEL=0,NEAR_WALL_TURB_MODEL=5
    CHARACTER(LABEL_LENGTH) :: PART_ID
    CHARACTER(LABEL_LENGTH) :: ID,TEXTURE_MAP,LEAK_PATH_ID(2)