From aeea2882e32da9d2b4001f119ac6cff988a7832e Mon Sep 17 00:00:00 2001 From: Trevor Hillebrand Date: Fri, 16 Dec 2022 10:20:45 -0800 Subject: [PATCH] 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. --- .../mpas-albany-landice/src/mode_forward/mpas_li_calving.F | 6 ++++++ 1 file changed, 6 insertions(+) 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 62bf8b048873..412ceca7c4e9 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 @@ -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