diff --git a/Source/dump.f90 b/Source/dump.f90
index 5bf9a6fcd40..57196a47b7c 100644
--- a/Source/dump.f90
+++ b/Source/dump.f90
@@ -7181,7 +7181,7 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
 REAL(EB), INTENT(IN) :: T,DT
 INTEGER, INTENT(IN) :: II,JJ,KK,IND,IND2,NM,Y_INDEX,Z_INDEX,PART_INDEX,VELO_INDEX,PIPE_INDEX,PROP_INDEX,REAC_INDEX,MATL_INDEX
 INTEGER, INTENT(IN), OPTIONAL :: ICC_IN,JCC_IN
-REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,DROPLET_COOLING,&
+REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,&
             Q_SUM,TMP_G,UU,VV,WW,VEL2,Y_MF_INT,PATHLENGTH,EXT_COEF,MASS_EXT_COEF,ZZ_FUEL,ZZ_OX,&
             VELSR,WATER_VOL_FRAC,RHS,DT_C,DT_E,T_RATIO,Y_E_LAG, H_G,H_G_SUM,CPBAR,CP,ZZ_GET(1:N_TRACKED_SPECIES),RCON,&
             EXPON,Y_SPECIES,MEC,Y_SPECIES2,Y_H2O,R_Y_H2O,R_DN,SGN,Y_ALL(N_SPECIES),H_S,D_Z_N(0:I_MAX_TEMP),&
@@ -7853,13 +7853,8 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
                AVG_DROP_DEN(I,J,K,LAGRANGIAN_PARTICLE_CLASS(NN)%ARRAY_INDEX)/LAGRANGIAN_PARTICLE_CLASS(NN)%DENSITY
          ENDDO
       ENDIF
-      IF (DV%TMP_L>PY%INITIAL_TEMPERATURE) THEN
-         DROPLET_COOLING = C_DIMARZO*VEL*WATER_VOL_FRAC
-      ELSE
-         DROPLET_COOLING = 0._EB
-      ENDIF
-      RHS      = ( VELSR*(TMP_G-DV%TMP_L) - PY%C_FACTOR*(DV%TMP_L-PY%INITIAL_TEMPERATURE) - DROPLET_COOLING )/PY%RTI
-      DV%TMP_L = DV%TMP_L + DT*RHS
+      RHS      = (VELSR*(TMP_G-DV%TMP_L) - PY%C_FACTOR*(DV%TMP_L-PY%INITIAL_TEMPERATURE) - C_DIMARZO*VEL*WATER_VOL_FRAC)/PY%RTI
+      DV%TMP_L = MAX(PY%INITIAL_TEMPERATURE , DV%TMP_L + DT*RHS)
       GAS_PHASE_OUTPUT_RES = DV%TMP_L - TMPM
 
    CASE(157) ! LINK TEMPERATURE