FDS Source : Remove computation of CP when CONSTANT_SPECIFIC_HEAT_RAT…

…IO=T, use CPOPR.

Source/divg.f90

@@ -427,21 +427,15 @@ SUBROUTINE DIVERGENCE_PART_1(T,DT,NM)
          ENDDO
       ENDDO
    ENDDO
-   IF (CONSTANT_SPECIFIC_HEAT_RATIO) THEN ! We have to compute CP, not done previously for constant CP ratio.
-      CP => WORK5
-      DO K=1,KBAR
-         DO J=1,JBAR
-            DO I=1,IBAR
-               IF (SOLID(CELL_INDEX(I,J,K))) CYCLE
-               ZZ_GET(1:N_TRACKED_SPECIES) = ZZP(I,J,K,1:N_TRACKED_SPECIES)
-               CALL GET_SPECIFIC_HEAT(ZZ_GET,CP(I,J,K),TMP(I,J,K))
-            ENDDO
-         ENDDO
-      ENDDO
-   ENDIF
    DEALLOCATE(ZZ_GET)
 
-   IF (SIM_MODE==LES_MODE) KP = KP + MAX(0._EB,(MU-MU_DNS))*CP*RPR
+   IF (SIM_MODE==LES_MODE) THEN
+      IF(.NOT.CONSTANT_SPECIFIC_HEAT_RATIO) THEN
+         KP = KP + MAX(0._EB,(MU-MU_DNS))*CP*RPR
+      ELSE
+         KP = KP + MAX(0._EB,(MU-MU_DNS))*CPOPR
+      ENDIF
+   ENDIF
 
    BOUNDARY_LOOP: DO IW=1,N_EXTERNAL_WALL_CELLS
       WC=>WALL(IW)