diff --git a/components/mpas-albany-landice/src/analysis_members/mpas_li_global_stats.F b/components/mpas-albany-landice/src/analysis_members/mpas_li_global_stats.F
index 961ee195ff61..0410f20d002f 100644
--- a/components/mpas-albany-landice/src/analysis_members/mpas_li_global_stats.F
+++ b/components/mpas-albany-landice/src/analysis_members/mpas_li_global_stats.F
@@ -369,7 +369,7 @@ subroutine li_compute_global_stats(domain, memberName, timeLevel, err)
             blockSumFloatingCalvingFlux = blockSumFloatingCalvingFlux + floatingCalvingThickness(iCell) * &
                areaCell(iCell) * rhoi / ( deltat / scyr )
 
-            ! mass loss due to calving (kg yr^{-1})
+            ! mass loss due to face-melting (kg yr^{-1})
             blockSumFaceMeltingFlux = blockSumFaceMeltingFlux + faceMeltingThickness(iCell) * &
                areaCell(iCell) * rhoi / ( deltat / scyr )
             blockSumGroundedFaceMeltingFlux = blockSumGroundedFaceMeltingFlux + groundedFaceMeltingThickness(iCell) * &
@@ -377,10 +377,6 @@ subroutine li_compute_global_stats(domain, memberName, timeLevel, err)
             blockSumFloatingFaceMeltingFlux = blockSumFloatingFaceMeltingFlux + floatingFaceMeltingThickness(iCell) * &
                areaCell(iCell) * rhoi / ( deltat / scyr )
 
-            ! mass loss due to face-melting (kg yr^{-1})
-            blockSumFaceMeltingFlux = blockSumFaceMeltingFlux + faceMeltingThickness(iCell) * &
-               areaCell(iCell) * rhoi / ( deltat / scyr )
-
             ! max surface speed
             if (surfaceSpeed(iCell) > blockMaxSurfaceSpeed) then
                blockMaxSurfaceSpeed = surfaceSpeed(iCell)