anchor.diff:  
       ! Scan domains for unreasonably large numbers or NaN issues
       subroutine check_for_overflow(md, flag, domain_ind)
  -        class(multidomain_type), intent(inout) :: md
  +        type(multidomain_type), intent(inout) :: md
           character(len=*), intent(in) :: flag
           integer(ip), optional, intent(in):: domain_ind
   
  -        logical :: throw_error
  +        integer :: throw_error
           integer(ip) :: i1, i2, i3, j, d1, d2
           real(dp), parameter :: bignum = HUGE(1.0)/2.0 ! Deliberate single-precision number
   
  
positive.code:         end do
  
          !DEBUG
          !call md%check_for_overflow('step-before-comms')
  
          if(.not. send_halos_immediately) then
              ! Do all the coarray communication in one hit
  
negative.code: 
              ! Step once
              call md%domains(j)%evolve_one_step(dt_local)
              !call check_for_overflow(md, 'inner', j)
  
              ! Report max_dt as the peak dt during the first timestep
              ! In practice max_dt may cycle between time-steps
  