diff --git a/physics/module_sf_noahmplsm.F90 b/physics/module_sf_noahmplsm.F90
index 1888a26e8..7813130b6 100644
--- a/physics/module_sf_noahmplsm.F90
+++ b/physics/module_sf_noahmplsm.F90
@@ -4246,7 +4246,7 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
 	end if
 
 ! canopy heat capacity
-        hcv = parameters%cbiom*vaie*cwat + canliq*cwat/denh2o + canice*cice/denice !j/m2/k
+        hcv = fveg*(parameters%cbiom*vaie*cwat + canliq*cwat/denh2o + canice*cice/denice) !j/m2/k
 
         b   = sav-irc-shc-evc-tr+pahv                          !additional w/m2
 !       a   = fveg*(4.*cir*tv**3 + csh + (cev+ctr)*destv) !volumetric heat capacity