diff --git a/Source/dump.f90 b/Source/dump.f90
index 57196a47b7c..e8feaab6462 100644
--- a/Source/dump.f90
+++ b/Source/dump.f90
@@ -7854,7 +7854,7 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
          ENDDO
       ENDIF
       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)
+      DV%TMP_L = MAX(MIN(TMP_G,PY%INITIAL_TEMPERATURE) , DV%TMP_L + DT*RHS)
       GAS_PHASE_OUTPUT_RES = DV%TMP_L - TMPM
 
    CASE(157) ! LINK TEMPERATURE