diff --git a/physics/satmedmfvdifq.F b/physics/satmedmfvdifq.F
index 6633708a3..ae9bf32f2 100644
--- a/physics/satmedmfvdifq.F
+++ b/physics/satmedmfvdifq.F
@@ -1542,7 +1542,7 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
                  ZCAN = ZCAN-0.5  !step down in-canopy resolution of 0.5m
               END DO !end loop on canopy layers
               EDDYVEST_INT = IntegrateTrapezoid((ZCANX(COUNTCAN:1:-1)
-     &                                /FCH),EDDYVESTX(COUNTCAN:1:-1))
+     &                                ),EDDYVESTX(COUNTCAN:1:-1)) / ZFL
               dkt(i,k)= (dkt(i,k)/EDDYVEST1) * EDDYVEST_INT !Scale dkt to resolved eddy diffusivity
               dkq(i,k)= (dkq(i,k)/EDDYVEST1) * EDDYVEST_INT !Scale dkq to resolved eddy diffusivity
               dku(i,k)= (dku(i,k)/EDDYVEST1) * EDDYVEST_INT !Scale dku to resolved eddy diffusivity