diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F
index 0fa0d6ed55e9..63bad581005f 100644
--- a/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F
+++ b/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F
@@ -259,19 +259,12 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
       err = ior(err, err_tmp)
       call mpas_timer_stop("advancing clock")
 
-!TODO: Determine whether grounded melting should in fact be called first
 ! === Ocean forcing extrapolation into ice-shelf cavities ===========
       call mpas_timer_start("ocean forcing extrapolation")
       call li_ocean_extrap_solve(domain, err_tmp)
       err = ior(err, err_tmp)
       call mpas_timer_stop("ocean forcing extrapolation")
 
-! === Face melting for grounded ice ===========
-      call mpas_timer_start("face melting for grounded ice")
-      call li_face_melt_grounded_ice(domain, err_tmp)
-      err = ior(err, err_tmp)
-      call mpas_timer_stop("face melting for grounded ice")
-
 ! === Basal melting for floating ice ===========
       call mpas_timer_start("basal melting for floating ice")
       call li_basal_melt_floating_ice(domain, err_tmp)
@@ -611,6 +604,12 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
       call mpas_dmpar_field_halo_exch(domain, 'vertexMask')
       call mpas_timer_stop("halo updates")
 
+! === Face melting for grounded ice ===========
+      call mpas_timer_start("face melting for grounded ice")
+      call li_face_melt_grounded_ice(domain, err_tmp)
+      err = ior(err, err_tmp)
+      call mpas_timer_stop("face melting for grounded ice")
+
 ! === Update bed topo =====================
 ! It's not clear when the best time to do this is.
 ! Seems cleaner to do it either before or after all of the time evolution of the ice