diff --git a/components/clm/src/betr/BetrBGCMod.F90 b/components/clm/src/betr/BetrBGCMod.F90
index 223f947c3c96..c3b074348c3f 100644
--- a/components/clm/src/betr/BetrBGCMod.F90
+++ b/components/clm/src/betr/BetrBGCMod.F90
@@ -717,7 +717,7 @@ subroutine do_tracer_advection(bounds, lbj, ubj, jtops, num_soilc, filter_soilc,
             qflx_adv_local(c,jtops(c)-1) = safe_div(qflx_adv(c,jtops(c)-1),aqu2bulkcef_mobile_col(c,jtops(c),j),eps=loc_eps)
             do l = jtops(c), ubj
                qflx_adv_local(c,l)     = safe_div(qflx_adv(c,l),aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
-               qflx_rootsoi_local(c,l) = safe_div(qflx_rootsoi(c,l),aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
+               qflx_rootsoi_local(c,l) = safe_div(qflx_rootsoi(c,l)*1.e-3_r8,aqu2bulkcef_mobile_col(c,l,j),eps=loc_eps)
             enddo
          enddo
 
diff --git a/components/clm/src/betr/betr_core/TracerParamsMod.F90 b/components/clm/src/betr/betr_core/TracerParamsMod.F90
index 5fb7fff80747..ed994144d3f9 100644
--- a/components/clm/src/betr/betr_core/TracerParamsMod.F90
+++ b/components/clm/src/betr/betr_core/TracerParamsMod.F90
@@ -1368,7 +1368,7 @@ subroutine diagnose_advect_water_flux(bounds, num_hydrologyc, filter_hydrologyc,
        if(j==nlevsoi)then
          qflx_adv(c,j) = qcharge(c) * 1.e-3_r8
        else
-         qflx_adv(c,j) = 1.e-3_r8 * (h2osoi_liq(c,j+1)-h2osoi_liq_copy(c,j+1))/dtime + qflx_adv(c,j+1) + qflx_rootsoi(c,j+1)
+         qflx_adv(c,j) = 1.e-3_r8 * (h2osoi_liq(c,j+1)-h2osoi_liq_copy(c,j+1))/dtime + qflx_adv(c,j+1) + qflx_rootsoi(c,j+1)*1.e-3_r8
        endif
 
      enddo
@@ -1380,7 +1380,7 @@ subroutine diagnose_advect_water_flux(bounds, num_hydrologyc, filter_hydrologyc,
      c = filter_hydrologyc(fc)
 
      !obtain the corrected infiltration
-     qflx_infl(c) = (h2osoi_liq(c,1)-h2osoi_liq_copy(c,1))/dtime + (qflx_rootsoi(c,1)+qflx_adv(c,1))*1.e3_r8
+     qflx_infl(c) = (h2osoi_liq(c,1)-h2osoi_liq_copy(c,1))/dtime + qflx_rootsoi(c,1) + qflx_adv(c,1)*1.e3_r8
      !the predicted net infiltration
      infl_tmp=qflx_gross_infl_soil(c)-qflx_gross_evap_soil(c)
      diff=qflx_infl(c)-infl_tmp
diff --git a/components/clm/src/biogeophys/SoilWaterMovementMod.F90 b/components/clm/src/biogeophys/SoilWaterMovementMod.F90
index 20758a2aefc4..edea194055dc 100644
--- a/components/clm/src/biogeophys/SoilWaterMovementMod.F90
+++ b/components/clm/src/biogeophys/SoilWaterMovementMod.F90
@@ -355,7 +355,7 @@ subroutine soilwater_zengdecker2009(bounds, num_hydrologyc, filter_hydrologyc, &
          qflx_infl         =>    waterflux_vars%qflx_infl_col       , & ! Input:  [real(r8) (:)   ]  infiltration (mm H2O /s)                          
          qflx_tran_veg_col =>    waterflux_vars%qflx_tran_veg_col   , & ! Input:  [real(r8) (:)   ]  vegetation transpiration (mm H2O/s) (+ = to atm)  
          qflx_tran_veg_pft =>    waterflux_vars%qflx_tran_veg_patch , & ! Input:  [real(r8) (:)   ]  vegetation transpiration (mm H2O/s) (+ = to atm)  
-         qflx_rootsoi      =>    waterflux_vars%qflx_rootsoi_col    , & ! Output: [real(r8) (:,:) ]  vegetation/soil water exchange (m H2O/s) (+ = to atm)
+         qflx_rootsoi      =>    waterflux_vars%qflx_rootsoi_col    , & ! Output: [real(r8) (:,:) ]  vegetation/soil water exchange (mm H2O/s) (+ = to atm)
 
          t_soisno          =>    temperature_vars%t_soisno_col        & ! Input:  [real(r8) (:,:) ]  soil temperature (Kelvin)                       
          )
@@ -786,7 +786,7 @@ subroutine soilwater_zengdecker2009(bounds, num_hydrologyc, filter_hydrologyc, &
       do j = 1, nlevsoi
          do fc = 1, num_hydrologyc
             c = filter_hydrologyc(fc)
-            qflx_rootsoi(c,j) = qflx_tran_veg_col(c) * rootr_col(c,j) * 1.e-3_r8       ![m H2O/s]
+            qflx_rootsoi(c,j) = qflx_tran_veg_col(c) * rootr_col(c,j)
          enddo
       enddo