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.
MALI-Dev · Dec 16, 2022 · aeea288 · aeea288
1 parent fe6b0e8
commit aeea288
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F
@@ -3538,6 +3538,12 @@ 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_timer_stop("halo updates")
+
         if (config_damage_rheology_coupling) then
             do iCell = 1, nCells
                 if (li_mask_is_floating_ice(cellMask(iCell))) then