program main
  implicit none
  integer, parameter :: n = 100
  real :: a(n)[*]
  integer :: i
  call random_number(a)
  sync all
  if (this_image() == 1) then
     do i = 2, num_images()
        a = a + a(:)[i]
     end do
     print *, "Mean :", sum(a) / (n * num_images())
  end if
end program