anchor.diff:      include 'mpif.h'
       integer :: ierr
   
  +!    call MPI_BARRIER(MPI_COMM_WORLD)
       call MPI_BCAST (N_states, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
       if (ierr /= MPI_SUCCESS) then
         print *,  'Unable to broadcast N_states'
         stop -1
       endif
   
  +!    call MPI_BARRIER(MPI_COMM_WORLD)
       call MPI_BCAST (N_det, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
       if (ierr /= MPI_SUCCESS) then
         print *,  'Unable to broadcast N_det'
         stop -1
       endif
   
  +!    call MPI_BARRIER(MPI_COMM_WORLD)
       call MPI_BCAST (psi_det_size, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
       if (ierr /= MPI_SUCCESS) then
         print *,  'Unable to broadcast psi_det_size'
  
positive.code:       TOUCH N_det_selectors
      endif
  
      call MPI_BCAST (energy, size(energy), MPI_DOUBLE_PRECISION, 0, MPI_COMM_WORLD, ierr)
      if (ierr /= MPI_SUCCESS) then
        print *,  'Unable to broadcast energy'
        stop -1
      endif
  
    IRP_ENDIF
  end
  
negative.code:       TOUCH N_det_generators
      endif
  
  !    call MPI_BARRIER(MPI_COMM_WORLD)
      call MPI_BCAST (N_det_selectors, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
      if (ierr /= MPI_SUCCESS) then
        print *,  'Unable to broadcast N_det_selectors'
  