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 index b90dcc3901b6..f63bf787b0c5 100644 --- a/components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F +++ b/components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F @@ -3535,6 +3535,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) @@ -3742,6 +3747,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