gatherNrfai.f

Go to the documentation of this file.

!            we note the last index may be actual max used
!
      do k = 1, ansites
         do j = 1, 4
            do l = 1,nfai
               do i = 1, nrbin
                  nrfaireca(i, l, j, k)=0.
               enddo
            enddo
         enddo
      enddo
      do k = 1, maxnoofsites
         do j = 1, 4
            do l = 1,nfai
               do i = 1, nrbin
                  nrfaialla(i, l, j, k)=0.
               enddo
            enddo
         enddo
      enddo
      do k = 1,  maxnoofsites
         do l = 1,nfai
            do i = 1, nrbin
               derfaia(i, l, k) = 0.
            enddo
         enddo
      enddo
      call mpi_reduce(nrfairec, nrfaireca, ansites*4*nfai*nrbin, 
     *   mpi_real, mpi_sum, 0, mpi_comm_world, icon)

      call mpi_reduce(nrfaiall, nrfaialla, maxnoofsites*4*nfai*nrbin, 
     *   mpi_real, mpi_sum, 0, mpi_comm_world, icon)

      call mpi_reduce(derfai, derfaia, maxnoofsites*nfai*nrbin, 
     *   mpi_real, mpi_sum, 0, mpi_comm_world, icon)
!          gatherred rec data is broadcasted; not needed; we work at rank 0 
!      call mpi_bcast(nrfairecA, ansites*4*nfai*nrbin, MPI_REAL,
!     *      0, MPI_COMM_WORLD, icon)