anchor.diff:                          do i = 1, size(all_recv_metadata,2, kind=ip) ! Every row in the receive metadata
                               ! Tell all_recv_metadata the row index in send_metadata
                               if(all_recv_metadata(IND,i,jj,k) == j .and. all_recv_metadata(IMG,i,jj,k) == ims) then
  -                                counter = counter + 1 
  +                                counter = counter + 1
                                   all_recv_metadata(7, i, jj, k) = 1.0_dp * counter
                               end if
                           end do
  
positive.code:             call communicate_p2p
              TIMER_STOP('comms1')
          end if
          
          ! Get the halo information from neighbours
          ! For coarray comms, we need to sync before to ensure the information is sent, and
          ! 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)
  
          !FIXME DEBUG
          !call md%check_for_overflow('step-after-comms')
  
          if(send_boundary_flux_data) then
  
negative.code:     ! of domains(i) on image j
      !
      ! @param domains the array of domains
      !
      subroutine create_all_bbox_and_dx(domains, periodic_xs, periodic_ys)
  
          type(domain_type), intent(inout) :: domains(:)
  