Merge pull request #13916 from drjfloyd/master

FDS Source: Fix array assignment for new SWORK allocation Issue #13913

Author: drjfloyd

Source/fire.f90

@@ -1617,7 +1617,7 @@ SUBROUTINE CONDENSATION_EVAPORATION(DT,NM)
 CALL POINT_TO_MESH(NM)
 
 ZZ_INTERIM=> SWORK1
-ZZ_INTERIM = ZZ
+ZZ_INTERIM(:,:,:,1:) = ZZ(:,:,:,1:) ! Lower bound may be zero for MW flux correction
 RHO_INTERIM => WORK1
 RHO_INTERIM = RHO
 TMP_INTERIM => WORK2

Source/part.f90

@@ -3634,7 +3634,7 @@ SUBROUTINE PARTICLE_MASS_ENERGY_TRANSFER(T,DT,NM)
 
    RHO_INTERIM => WORK1 ; RHO_INTERIM = RHO
    TMP_INTERIM => WORK2 ; TMP_INTERIM = TMP
-   ZZ_INTERIM  => SWORK1 ; ZZ_INTERIM = ZZ
+   ZZ_INTERIM  => SWORK1 ; ZZ_INTERIM(:,:,:,1:) = ZZ(:,:,:,1:)
 
 ENDIF