Fix accidental doubling of faceMeltingFlux

Fix accidental doubling of faceMeltingFlux when calculating global stats.
MALI-Dev · Oct 18, 2023 · fb51715 · fb51715
1 parent 479ba60
commit fb51715
Showing 1 changed file with 1 addition and 5 deletions.
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
@@ -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)