anchor.diff:              ! Do the remaining sub-steps
               do i = 2, nt ! Loop never runs if nt < 2
                   call md%domains(j)%evolve_one_step(dt_local)
  -                !call check_for_overflow(md, 'innerB', j)
  +#ifdef TRACK_MULTIDOMAIN_STABILITY
  +                call check_multidomain_stability(md, 'innerB', j)
  +#endif
               end do
               md%domains(j)%max_dt = max_dt_store
   
  
positive.code:         ! there, leading to tiny numerical differences in the time.
          ! To avoid any issues, we explicitly set them to the same number
  
          !DEBUG
          !call check_for_overflow(md, 'start')
  
          ! Evolve every domain, sending the data to communicate straight after
          ! the evolve
  
negative.code: 
  
      ! Scan domains for unreasonably large numbers or NaN issues
      subroutine check_multidomain_stability(md, flag, domain_ind)
          type(multidomain_type), intent(inout) :: md
          character(len=*), intent(in) :: flag
          integer(ip), optional, intent(in):: domain_ind
  
          integer :: throw_error
          integer(ip) :: i1, i2, i3, j, d1, d2, ecw
          real(dp), parameter :: bignum = HUGE(1.0)/2.0 ! Deliberate single-precision number
  
          if(present(domain_ind)) then
  