Update halos for damage before applying to stiffnessFactor

Update halos for damage before applying to stiffnessFactor. Testing
showed imprinting of the decomposition on damage and stiffnessFactor
fields before applying this halo update.

Author: trhille

components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F

@@ -3331,6 +3331,11 @@ subroutine li_calculate_damage(domain, err)
             err = ior(err, 1)
         endif
 
+        call mpas_timer_start("halo updates")
+        call mpas_dmpar_field_halo_exch(domain, 'damage')
+        call mpas_dmpar_field_halo_exch(domain, 'damageNye')
+        call mpas_timer_stop("halo updates")
+
         if (config_print_calving_info) then
            ! End of routine accounting/reporting
            localMinInfo(1) = minval(damageSource)
@@ -3538,6 +3543,13 @@ subroutine li_finalize_damage_after_advection(domain, err)
             damage = 1.0_RKIND
         end where
 
+         ! Halo update for damage before applying to stiffnessFactor
+         ! NOTE: THIS WILL NOT WORK ON MULTIPLE BLOCKS PER PROCESSOR
+         call mpas_timer_start("halo updates")
+         call mpas_dmpar_field_halo_exch(domain, 'damage')
+         call mpas_dmpar_field_halo_exch(domain, 'damageNye')
+         call mpas_timer_stop("halo updates")
+
         if (config_damage_rheology_coupling) then
             do iCell = 1, nCells
                 if (li_mask_is_floating_ice(cellMask(iCell))) then