Skip to content

Commit fe6b0e8

Browse files
committed
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.
1 parent 637ad9d commit fe6b0e8

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

components/mpas-albany-landice/src/mode_forward/mpas_li_core.F

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -758,7 +758,7 @@ function li_core_initial_solve(domain) result(err)
758758
! initialize damage using stiffnessFactor
759759
call mpas_pool_get_array(velocityPool, 'stiffnessFactor', stiffnessFactor)
760760
call mpas_pool_get_array(geometryPool, 'damage', damage)
761-
if ( config_initialize_damage_with_stiffnessFactor ) then
761+
if ( config_initialize_damage_with_stiffnessFactor .and. (.not. config_do_restart) ) then
762762
damage(:) = 1.0_RKIND - stiffnessFactor(:)
763763
764764
where ( damage < 0.0_RKIND )

0 commit comments

Comments
 (0)