updated some comments, removed unnecessary flushing

src/utils/clmfates_interfaceMod.F90

@@ -2754,7 +2754,16 @@ end subroutine wrap_WoodProducts
 
  subroutine wrap_canopy_radiation(this, bounds_clump, nc, fcansno, surfalb_inst)
 
-   ! Arguments
+
+    ! Pass boundary conditions (zenith angle, ground albedo and snow-cover)
+    ! to FATES, and have fates call normalized radiation scattering.
+    ! These methods are normalized because it is the zenith of the next
+    ! timestep, and because of atmospheric model coupling, we need to provide
+    ! an albedo. This normalized solution will be scaled by the actual
+    ! downwelling radiation during the wrap sunshade fraction call
+
+
+    ! Arguments
     class(hlm_fates_interface_type), intent(inout) :: this
     type(bounds_type),  intent(in)             :: bounds_clump
     integer            , intent(in)            :: nc ! clump index
@@ -2834,17 +2843,6 @@ subroutine wrap_canopy_radiation(this, bounds_clump, nc, fcansno, surfalb_inst)
           ftii(p,:) = this%fates(nc)%bc_out(s)%ftii_parb(ifp,:)
 
        end do
-       ! Flush other patches a default value
-       do p = col%patchi(c)+ifp+1,col%patchf(c)
-          albd(p,:) = 0._r8
-          albi(p,:) = 0._r8
-          fabd(p,:) = 0._r8
-          fabi(p,:) = 0._r8
-          ftdd(p,:) = 1._r8
-          ftid(p,:) = 1._r8
-          ftii(p,:) = 1._r8
-       end do
-
     end do
 
   end associate