From fb51715ee3855fe03edb46f7103a5a17410e84fe Mon Sep 17 00:00:00 2001 From: Trevor Ray Hillebrand Date: Wed, 18 Oct 2023 17:05:25 -0600 Subject: [PATCH] Fix accidental doubling of faceMeltingFlux Fix accidental doubling of faceMeltingFlux when calculating global stats. --- .../src/analysis_members/mpas_li_global_stats.F | 6 +----- 1 file changed, 1 insertion(+), 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 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)