Update thickness and tracer halos before velocity solve

Update thickness and tracer halos before velocity solve. This allows eismint2 decomposition test to pass validation.
MALI-Dev · Oct 6, 2023 · 35977fc · 35977fc
1 parent 72f3b83
commit 35977fc
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_time_integration_fe_rk.F
@@ -480,6 +480,19 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
          floatingBasalMassBalAccum = floatingBasalMassBalAccum + rkTendWeights(rkStage) * floatingBasalMassBalApplied
          fluxAcrossGroundingLineAccum = fluxAcrossGroundingLineAccum + rkTendWeights(rkStage) * fluxAcrossGroundingLine
 
+         ! Halo updates
+         call mpas_timer_start("halo updates")
+
+         call mpas_dmpar_field_halo_exch(domain, 'thickness')
+         call mpas_dmpar_field_halo_exch(domain, 'temperature')
+         call mpas_dmpar_field_halo_exch(domain, 'waterFrac')
+         call mpas_dmpar_field_halo_exch(domain, 'enthalpy')
+         call mpas_dmpar_field_halo_exch(domain, 'damage')
+         call mpas_dmpar_field_halo_exch(domain, 'passiveTracer2d')
+
+         call mpas_timer_stop("halo updates")
+
+
          ! Update velocity for each RK step
          ! === Solve Velocity =====================
          call li_velocity_solve(domain, solveVelo=.true., err=err_tmp)