Merge pull request #14613 from mcgratta/master

FDS Source: Update all rad angles if OBST created or removed

Author: mcgratta

Source/cons.f90

@@ -177,6 +177,7 @@ MODULE GLOBAL_CONSTANTS
 
 LOGICAL :: HVAC_DEBUG=.FALSE.               !< Output known hvac values to smokeview
 LOGICAL :: RADIATION=.TRUE.                 !< Perform radiation transport
+LOGICAL :: UPDATE_ALL_ANGLES=.FALSE.        !< Update all radiation angles the next time the solver is called
 LOGICAL :: INCLUDE_PYROLYSIS=.FALSE.        !< Solid phase pyrolysis is included in the simulation
 LOGICAL :: EXCHANGE_RADIATION=.FALSE.       !< Do an MPI radiation exchange at this time step
 LOGICAL :: EXCHANGE_OBST_MASS=.FALSE.       !< Exchange mass loss information for obstructions bordering interpolated meshes

Source/main.f90

@@ -1573,7 +1573,7 @@ SUBROUTINE CREATE_OR_REMOVE_OBSTRUCTIONS
    DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
       CALL REASSIGN_WALL_CELLS(T,NM)
    ENDDO
-
+   UPDATE_ALL_ANGLES = .TRUE.  ! Update all radiation angles the next time radiation is updated
 ENDIF
 
 END SUBROUTINE CREATE_OR_REMOVE_OBSTRUCTIONS

Source/radi.f90

@@ -3525,7 +3525,7 @@ SUBROUTINE RADIATION_FVM
 
 ! Check if it time to update radiation intensity field
 
-IF ( MOD(RAD_CALL_COUNTER,TIME_STEP_INCREMENT)==0 .OR. INITIALIZATION_PHASE .OR. ICYC==1) THEN
+IF ( MOD(RAD_CALL_COUNTER,TIME_STEP_INCREMENT)==0 .OR. INITIALIZATION_PHASE .OR. ICYC==1 .OR. UPDATE_ALL_ANGLES) THEN
    UPDATE_INTENSITY   = .TRUE.
    EXCHANGE_RADIATION = .TRUE.
 ELSE
@@ -3558,7 +3558,7 @@ SUBROUTINE RADIATION_FVM
 IF (UPDATE_INTENSITY) THEN
    DO IW=1,N_INTERNAL_WALL_CELLS+N_EXTERNAL_WALL_CELLS
       WC => WALL(IW)
-      IF (WC%B1_INDEX==0) CYCLE
+      IF (WC%B1_INDEX==0 .OR. WC%BOUNDARY_TYPE==NULL_BOUNDARY) CYCLE
       B1 => BOUNDARY_PROP1(WC%B1_INDEX)
       B1%Q_RAD_IN = 0._EB
       SF  => SURFACE(WALL(IW)%SURF_INDEX)
@@ -3949,7 +3949,7 @@ SUBROUTINE RADIATION_FVM
       ! If updating intensities first time, sweep ALL angles
 
       N_UPDATES = 1
-      IF (INITIALIZATION_PHASE .OR. ICYC==1) N_UPDATES = ANGLE_INCREMENT
+      IF (INITIALIZATION_PHASE .OR. ICYC==1 .OR. UPDATE_ALL_ANGLES) N_UPDATES = ANGLE_INCREMENT
 
       UPDATE_LOOP: DO I_UIID = 1,N_UPDATES
 
@@ -4580,6 +4580,10 @@ SUBROUTINE RADIATION_FVM
 
 ENDIF
 
+! UPDATE_ALL_ANGLES is a one-time only logical. If needed again, it will be reset again.
+
+UPDATE_ALL_ANGLES = .FALSE.
+
 END SUBROUTINE RADIATION_FVM