Debug waterPressureSmooth

fixes minor bug in waterPRessureSmooth calculation
glemieux · May 22, 2024 · 1f94eb5 · 1f94eb5
1 parent e76d9a9
commit 1f94eb5
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F
@@ -2674,7 +2674,7 @@ subroutine calc_waterPressureSmooth(block,err)
          do iCell = 1, nCells
             if (li_mask_is_grounded_ice(cellMask(iCell))) then !smooth over all grounded cells
 
-               totalPressure = pressureField(iCell)
+               totalPressure = pressureField(iCell)*areaCell(iCell) 
                totalArea = areaCell(iCell)
           
                do jCell = 1, nEdgesOnCell(iCell)
@@ -2688,7 +2688,7 @@ subroutine calc_waterPressureSmooth(block,err)
                   endif
                end do
           
-               waterPressureSmooth(iCell) = totalPressure / totalArea !area-weighted average height
+               waterPressureSmooth(iCell) = totalPressure / totalArea !area-weighted average pressure
             endif
       end do