FDS Source: Issue firemodels#13440. Remove unneeded loop

Source/radi.f90

@@ -2789,7 +2789,7 @@ SUBROUTINE INIT_RADIATION
 USE WSGG_ARRAYS
 REAL(EB) :: THETAUP,THETALOW,PHIUP,PHILOW,F_THETA,PLANCK_C2,KSI,LT,RCRHO,YY,YY2,BBF,AP0,AMEAN,RADIANCE,TRANSMISSIVITY,X_N2,&
             THETA,PHI
-INTEGER  :: N,I,J,K,IPC,IZERO,NN,NI,II,JJ,IIM,JJM,IBND,NS,NS2,NRA,NSB,RADCAL_TEMP(16)=0,RCT_SKIP=-1,OR_IN,I1,I2,IO
+INTEGER  :: N,I,J,K,IPC,IZERO,NN,NI,II,JJ,IIM,JJM,IBND,NS,NS2,NRA,NSB,RADCAL_TEMP(16)=0,RCT_SKIP=-1,IO
 TYPE (LAGRANGIAN_PARTICLE_CLASS_TYPE), POINTER :: LPC
 REAL(EB), ALLOCATABLE, DIMENSION(:) :: COSINE_ARRAY
 TYPE (RAD_FILE_TYPE), POINTER :: RF
@@ -3391,22 +3391,6 @@ SUBROUTINE INIT_RADIATION
 ! Determine angle factors for Lagrangian particles with ORIENTATION
 
 IF (SOLID_PARTICLES) THEN
-   PARTICLE_CLASS_LOOP: DO IPC=1,N_LAGRANGIAN_CLASSES
-      LPC => LAGRANGIAN_PARTICLE_CLASS(IPC)
-      IF (LPC%N_ORIENTATION==0) CYCLE PARTICLE_CLASS_LOOP
-      I1 = LPC%ORIENTATION_INDEX
-      I2 = LPC%ORIENTATION_INDEX+LPC%N_ORIENTATION-1
-      ALLOCATE(COSINE_ARRAY(I1:I2))
-      ANGLE_LOOP: DO N=1,NRA
-         ORIENTATION_LOOP: DO OR_IN=I1,I2
-            COSINE_ARRAY(OR_IN) = ORIENTATION_VECTOR(1,OR_IN)*DLX(N) + &
-                                  ORIENTATION_VECTOR(2,OR_IN)*DLY(N) + &
-                                  ORIENTATION_VECTOR(3,OR_IN)*DLZ(N)
-         ENDDO ORIENTATION_LOOP
-      ENDDO ANGLE_LOOP
-      DEALLOCATE(COSINE_ARRAY)
-   ENDDO PARTICLE_CLASS_LOOP
-
    ALLOCATE(COSINE_ARRAY(1:NRA))
    ALLOCATE(NEAREST_RADIATION_ANGLE(N_ORIENTATION_VECTOR))
    ALLOCATE(VIEW_ANGLE_AREA(N_ORIENTATION_VECTOR))
@@ -3424,9 +3408,7 @@ SUBROUTINE INIT_RADIATION
       NEAREST_RADIATION_ANGLE(IO) = MINLOC(COSINE_ARRAY,DIM=1)
       VIEW_ANGLE_AREA(IO) = PI/VIEW_ANGLE_AREA(IO)
    ENDDO
-
    DEALLOCATE(COSINE_ARRAY)
-
 ENDIF
 
 ! Allocate array needed by angle-specific RADF output files