From 8fe44f0f5e72bba16fc0a4e9af95fec25833e261 Mon Sep 17 00:00:00 2001 From: Trevor Ray Hillebrand Date: Thu, 19 Oct 2023 20:34:19 -0600 Subject: [PATCH] Update layerThickness halos after advection Update layerThickness halos after advection. This reduces the error in a test with 500 m/yr face-melt, so it seems to be necessary, but it does not solve the issue entirely. --- .../mpas-albany-landice/src/mode_forward/mpas_li_advection.F | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_advection.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_advection.F index 9b9f416164ff..bab4ec37f6c9 100644 --- a/components/mpas-albany-landice/src/mode_forward/mpas_li_advection.F +++ b/components/mpas-albany-landice/src/mode_forward/mpas_li_advection.F @@ -560,6 +560,10 @@ subroutine li_advection_thickness_tracers(& enddo endif + ! Update halos after advection. + call mpas_timer_start("halo updates") + call mpas_dmpar_field_halo_exch(domain, 'layerThickness') + call mpas_timer_stop("halo updates") ! Calculate dynamicThickening (layerThickness is updated by advection at this point, while thickness is still old) dynamicThickening = (sum(layerThickness, 1) - thickness) / dt * scyr ! units of m/yr