diff --git a/biogeophys/FatesPlantHydraulicsMod.F90 b/biogeophys/FatesPlantHydraulicsMod.F90
index 2c474ccc43..b38d58a2d3 100644
--- a/biogeophys/FatesPlantHydraulicsMod.F90
+++ b/biogeophys/FatesPlantHydraulicsMod.F90
@@ -3351,7 +3351,6 @@ subroutine ImTaylorSolve1D(site_hydr,cohort,cohort_hydr,dtime,q_top, &
         ! cohort_hydr%l_aroot_layer(ilayer) is units [m/plant]
         ! site_hydr%l_aroot_layer(ilayer) is units [m/site]
  
-        ! Junyan added if statement to handle zero l_aroot_layer condition
         if (site_hydr%l_aroot_layer(ilayer)<nearzero) then
            l_aroot_layer = 2*nearzero
         else