Do not initialize damage from stiffnessFactor on restarts

On a restart, just use the damage field in the restart file, rather than calculating it from stiffnessFactor.
MALI-Dev · Dec 15, 2022 · fe6b0e8 · fe6b0e8
1 parent 637ad9d
commit fe6b0e8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_core.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_core.F
@@ -758,7 +758,7 @@ function li_core_initial_solve(domain) result(err)
       ! initialize damage using stiffnessFactor
       call mpas_pool_get_array(velocityPool, 'stiffnessFactor', stiffnessFactor)
       call mpas_pool_get_array(geometryPool, 'damage', damage)
-      if ( config_initialize_damage_with_stiffnessFactor ) then
+      if ( config_initialize_damage_with_stiffnessFactor .and. (.not. config_do_restart) ) then
          damage(:) = 1.0_RKIND - stiffnessFactor(:)
 
          where ( damage < 0.0_RKIND )