anchor.diff:          ! A detailed domain which shares a physical boundary with the global
           ! domain. NOTE: The flow solver doesn't work with this configuration.
           !
  -        md%domains(2)%lower_left = global_ll 
  +        md%domains(2)%lower_left = global_ll
           md%domains(2)%lw = nint(inner_lw/md%domains(1)%dx)*md%domains(1)%dx ! Is a multiple of parent domain dx
           md%domains(2)%dx_refinement_factor = real(nest_ratio, dp)
           md%domains(2)%dx = md%domains(1)%dx / md%domains(2)%dx_refinement_factor
  
positive.code:     subroutine test_multidomain_global
  
  
          ! Type holding all domains 
          type(multidomain_type) :: md
  
          ! Change this to decrease the cell size by mesh_refine (i.e. for convergence testing)
  
negative.code:         !
  
          do j = 1, nd_local
  
              !
              ! Setup sends
              !
  