anchor.diff:  program main
  -  use mpi
  +  use mpi_f08
     use occa
     use, intrinsic :: iso_fortran_env, only : stdout=>output_unit, &
                                               stderr=>error_unit
  
positive.code:   implicit none
  
    integer :: ierr, id
    integer :: myid, npes, gcomm, tag ! MPI variables
    integer, dimension(2) :: request
    integer, dimension(MPI_STATUS_SIZE) :: status
    integer :: otherID, offset
    integer(occaUDim_t) :: iu
    integer(occaUDim_t) :: entries = 8
  
negative.code:   call flush(stdout)
    ab_sum = myid
    ab_gather = sum(ab_ptr)
    call MPI_Gather(ab_gather, 1, MPI_REAL4, ab_sum, 1, MPI_REAL4, 0, gcomm, ierr)
    if (myid == 0) then
      if (abs(ab_sum(myid) - ab_sum(otherID)) > 1.0e-8) stop "*** Wrong result ***"
    end if
  