diff --git a/Source/fire.f90 b/Source/fire.f90
index 0729337e8b5..953bb8ba95e 100644
--- a/Source/fire.f90
+++ b/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
diff --git a/Source/part.f90 b/Source/part.f90
index 527451a667f..5b1dbd156c4 100644
--- a/Source/part.f90
+++ b/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