diff --git a/ufs/glc_elevclass_mod.F90 b/ufs/glc_elevclass_mod.F90
index 6524f064f..626bb3ee0 100644
--- a/ufs/glc_elevclass_mod.F90
+++ b/ufs/glc_elevclass_mod.F90
@@ -37,6 +37,7 @@ subroutine glc_get_elevation_classes_without_bareland(glc_topo, glc_elevclass, l
     real(r8), intent(in)  :: glc_topo(:)      ! topographic height
     integer , intent(out) :: glc_elevclass(:) ! elevation class
     integer , intent(in)  :: logunit
+    glc_elevclass = 0
   end subroutine glc_get_elevation_classes_without_bareland
 
   !-----------------------------------------------------------------------
@@ -45,6 +46,7 @@ subroutine glc_get_elevation_classes_with_bareland(glc_ice_covered, glc_topo, gl
     real(r8), intent(in)  :: glc_topo(:)        ! ice topographic height
     integer , intent(out) :: glc_elevclass(:)   ! elevation class
     integer , intent(in)  :: logunit
+    glc_elevclass = 0
   end subroutine glc_get_elevation_classes_with_bareland
 
   !-----------------------------------------------------------------------
@@ -57,11 +59,12 @@ end function glc_mean_elevation_virtual
 
   !-----------------------------------------------------------------------
   subroutine glc_get_fractional_icecov(nec, glc_topo, glc_icefrac, glc_icefrac_ec, logunit)
-    integer , intent(in)  :: nec              ! number of elevation classes 
+    integer , intent(in)  :: nec              ! number of elevation classes
     real(r8), intent(in)  :: glc_topo(:)      ! topographic height
     real(r8), intent(in)  :: glc_icefrac(:)
     real(r8), intent(out) :: glc_icefrac_ec(:,:)
     integer , intent(in)  :: logunit
+    glc_icefrac_ec = 0.0_r8
   end subroutine glc_get_fractional_icecov
 
 end module glc_elevclass_mod