AllGETHERV_GSD.f

      SUBROUTINE allgetherv(GRID1)
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .
! SUBPROGRAM:    AllGETHERV VERT INTRP OF MODEL LVLS TO PRESSURE
!   PRGRMMR: MING HU           ORG: GSD     DATE: 2012-01-01
!
! ABSTRACT:
!   .
!
! PROGRAM HISTORY LOG:
!
!    21-09-02  Bo Cui - Decompose UPP in X direction               
      
     use ctlblk_mod, only : im,jm,num_procs,me,jsta,jend,ista,iend,mpi_comm_comp

     implicit none

     include "mpif.h"

!
    integer i,j,ij
    integer ierr

     REAL grid1(im,jm)
     REAL ibufrecv(im*jm)
     REAL ibufsend((iend-ista+1)*(jend-jsta+1))
     integer sendcount,recvcounts(num_procs),displs(num_procs)
!
!     write(*,*) 'check mpi', im,jm,num_procs,me,jsta,jend
     sendcount=(iend-ista+1)*(jend-jsta+1)
     call mpi_allgather(sendcount, 1, mpi_integer, recvcounts,1 , &
                mpi_integer, mpi_comm_comp, ierr)
     displs(1)=0
     do i=2,num_procs
         displs(i)=displs(i-1)+recvcounts(i-1)
     enddo
!
!     write(*,*) me,'RECVCOUNTS=',RECVCOUNTS
!     write(*,*) me,'DISPLS=',DISPLS
! 
     ij=0
     ibufsend=0.0
     do j=jsta,jend
        do i=ista,iend
           ij=ij+1
           ibufsend(ij)=grid1(i,j)
        enddo
     enddo
     if(ij /= recvcounts(me+1)) then
        write(*,*) 'Error: send account is not equal to receive account',me,ij,recvcounts(me+1)
     endif
  
     call mpi_allgatherv(ibufsend, ij, mpi_real, ibufrecv, recvcounts,displs, &
                mpi_real, mpi_comm_comp, ierr)

     ij=0
     do j=1,jm
        do i=1,im
           ij=ij+1
           grid1(i,j)=ibufrecv(ij)
        enddo
     enddo
!
!     END OF ROUTINE.
!
      RETURN
      END