anchor.diff:          ! also sync after to ensure it is not overwritten before it is used
           call md%recv_halos(sync_before=sync_before_recv, sync_after=sync_after_recv)
   
  -        !DEBUG
  -        !call check_for_overflow(md, 'step-after-comms')
  -
  +#ifdef TRACK_MULTIDOMAIN_STABILITY
  +        call check_multidomain_stability(md, 'step-after-comms')
  +#endif
           if(send_boundary_flux_data) then
               TIMER_START('nesting_flux_correction')
               ! Do flux correction
  
positive.code:             end do
              TIMER_STOP('nesting_flux_correction')
  
              !DEBUG
              !call check_for_overflow(md, 'step-after-flux_correct')
          end if
  
      end subroutine
  
negative.code:             md%domains(j)%time = tend
          end do
  
  #ifdef TRACK_MULTIDOMAIN_STABILITY
          call check_multidomain_stability(md, 'step-before-comms')
  #endif
  
          if(.not. send_halos_immediately) then
              ! Do all the coarray communication in one hit
  