diff --git a/phys/module_sf_slab.F b/phys/module_sf_slab.F
index 7304b5af44..71527aba0e 100644
--- a/phys/module_sf_slab.F
+++ b/phys/module_sf_slab.F
@@ -407,7 +407,7 @@ SUBROUTINE SLAB1D(J,T1D,QV1D,P1D,FLHC,FLQC,                  &
                 ESG=SVP1*EXP(SVP2*(TG0(I)-SVPT0)/(TG0(I)-SVP3))
                 QSG=EP2*ESG/(PSFC(I)-ESG)
                 THG=TSK(I)*(100./PSFC(I))**ROVCP
-                HFX(I)=FLHC(I)*(THG-THX(I))
+                HFX(I)=FLHC(I)*CAPG(I)*(THG-THX(I))
                 QFX(I)=FLQC(I)*(QSG-QX(I))
                 LH(I)=QFX(I)*XLV
         QS(I)=HFX(I)+QFX(I)*XLV                                
@@ -440,7 +440,7 @@ SUBROUTINE SLAB1D(J,T1D,QV1D,P1D,FLHC,FLQC,                  &
                     -SVP3))                                                      
                 QSG=EP2*ESG/(PS-ESG)                                             
 !     UPDATE FLUXES FOR NEW GROUND TEMPERATURE                                   
-                HFXT=FLHC(I)*(THG-THX(I))                                     
+                HFXT=FLHC(I)*CAPG(I)*(THG-THX(I))                                    
                 QFXT=FLQC(I)*(QSG-QX(I))
                 QS(I)=HFXT+QFXT*XLV                                
 !     SUM HFX AND QFX OVER SOIL TIMESTEPS