Calculate layerThickness after calving in RK loop

Calculate layerThickness after calving in RK loop, which is necessary when calculating weighted averages for RK integration.
MALI-Dev · Nov 2, 2023 · 5e23b98 · 5e23b98
1 parent 813ee05
commit 5e23b98
Showing 1 changed file with 2 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
@@ -403,6 +403,8 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
             call mpas_pool_get_array(thermalPool, 'temperature', temperature)
             call mpas_pool_get_array(thermalPool, 'waterFrac', waterFrac)
 
+            ! get layerThickness after calving
+            call li_calculate_layerThickness(meshPool, thickness, layerThickness)
             layerThicknessTmp(:,:) = layerThickness(:,:)
             layerThickness(:,:) = rkSSPweights(rkStage) * layerThicknessPrev(:,:) + &
                                      (1.0_RKIND - rkSSPweights(rkStage)) * layerThickness(:,:)