Skip to content

Commit

Permalink
Restore calving front within RK loop
Browse files Browse the repository at this point in the history
Restore calving front within RK loop to prevent velocity solutions
on extended ice geometry.
  • Loading branch information
trhille committed Oct 30, 2023
1 parent 8d5f28e commit 769dc33
Showing 1 changed file with 5 additions and 6 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -475,6 +475,11 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
call mpas_timer_stop("halo updates")
if (config_restore_calving_front) then
! restore the calving front to its initial position before velocity solve.
call li_restore_calving_front(domain, err_tmp)
err = ior(err, err_tmp)
endif
! Update velocity for each RK step
! === Solve Velocity =====================
Expand Down Expand Up @@ -531,12 +536,6 @@ subroutine li_time_integrator_forwardeuler_rungekutta(domain, err)
call li_calve_ice(domain, err_tmp)
err = ior(err, err_tmp)
if (config_restore_calving_front) then
! restore the calving front to its initial position; calving options are ignored
call li_restore_calving_front(domain, err_tmp)
err = ior(err, err_tmp)
endif
call mpas_timer_stop("calve_ice")
call mpas_timer_start("halo updates")
Expand Down

0 comments on commit 769dc33

Please sign in to comment.