Fix accidental doubling of faceMeltingFlux

Fix accidental doubling of faceMeltingFlux when calculating global
stats.

Author: trhille

components/mpas-albany-landice/src/analysis_members/mpas_li_global_stats.F

@@ -369,18 +369,14 @@ 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) * &
                areaCell(iCell) * rhoi / ( deltat / scyr )
             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)