diff --git a/src/SurfaceEnergyFluxVegetatedMod.F90 b/src/SurfaceEnergyFluxVegetatedMod.F90
index 5bf43dae..1283553a 100644
--- a/src/SurfaceEnergyFluxVegetatedMod.F90
+++ b/src/SurfaceEnergyFluxVegetatedMod.F90
@@ -305,14 +305,14 @@ subroutine SurfaceEnergyFluxVegetated(noahmp)
        EnergyResTmp         = RadSwAbsVeg - RadLwNetCanopy - HeatSensibleCanopy - &
                               HeatLatentCanEvap - HeatLatentCanTransp + HeatPrecipAdvCanopy
        FluxTotCoeff         = VegFrac * (4.0*LwCoeffCan*TemperatureCanopy**3 + ShCoeff + &
-                                        (LhCoeff+TranspHeatCoeff)*VapPresSatCanTempD) + HeatCapacCan/MainTimeStep     ! volumetric heat capacity
+                                        (LhCoeff+TranspHeatCoeff)*VapPresSatCanTempD + HeatCapacCan/MainTimeStep)     ! volumetric heat capacity
        TemperatureCanChg    = EnergyResTmp / FluxTotCoeff
        ! update fluxes with temperature change
-       RadLwNetCanopy       = RadLwNetCanopy + VegFrac * 4.0 * LwCoeffCan * TemperatureCanopy**3 * TemperatureCanChg
-       HeatSensibleCanopy   = HeatSensibleCanopy + VegFrac * ShCoeff * TemperatureCanChg
-       HeatLatentCanEvap    = HeatLatentCanEvap + VegFrac * LhCoeff * VapPresSatCanTempD * TemperatureCanChg
-       HeatLatentCanTransp  = HeatLatentCanTransp  + VegFrac * TranspHeatCoeff * VapPresSatCanTempD * TemperatureCanChg
-       HeatCanStorageChg    = TemperatureCanChg * HeatCapacCan / MainTimeStep   ! canopy heat storage change [W/m2]
+       RadLwNetCanopy       = RadLwNetCanopy      + VegFrac * 4.0 * LwCoeffCan * TemperatureCanopy**3 * TemperatureCanChg
+       HeatSensibleCanopy   = HeatSensibleCanopy  + VegFrac * ShCoeff * TemperatureCanChg
+       HeatLatentCanEvap    = HeatLatentCanEvap   + VegFrac * LhCoeff * VapPresSatCanTempD * TemperatureCanChg
+       HeatLatentCanTransp  = HeatLatentCanTransp + VegFrac * TranspHeatCoeff * VapPresSatCanTempD * TemperatureCanChg
+       HeatCanStorageChg    = VegFrac * HeatCapacCan / MainTimeStep * TemperatureCanChg  ! canopy heat storage change [W/m2]
        ! update vegetation temperature
        TemperatureCanopy    = TemperatureCanopy + TemperatureCanChg
       !TemperatureCanopyAir = TempShGhTmp + ExchCoeffShFrac * TemperatureCanopy                        ! canopy air T; update here for consistency