pkg/compon_communic/compsend_r4.F

C $Header:  $
C $Name:  $

!=======================================================================
      subroutine compsend_r4( dataname, Ni,Oi,Nj,Oj, arr )
      implicit none
! Arguments
      character*(*) dataname
      integer Ni,Oi,Nj,Oj
      real*4 arr(1-Oi:Ni+Oi,1-Oj:Nj+Oj)
! Predefined constants/arrays
#include "CPLR_SIG.h"
! MPI variables
#include "mpif.h"
      integer count,datatype,dest,tag,comm,ierr
! Functions
      integer generate_tag
! Local
      integer i,j,ij,Io,Jo
!     ------------------------------------------------------------------

      if (HEADER_SIZE+Ni*Nj.gt.MAX_R4_BUFLEN)
     &    stop 'compsend_r4: Nx*Ny too big'

! Set up buffer
      Io=my_tile_i0(1)
      Jo=my_tile_j0(1)
      r4buf(1)=float(Io)
      r4buf(2)=float(Ni)
      r4buf(3)=float(Jo)
      r4buf(4)=float(Nj)
      call mitcplr_char2real( dataname, r4buf(9) )

! Copy interior of tile to buffer
      do j=1,Nj
       do i=1,Ni
        ij=HEADER_SIZE+i+Ni*(j-1)
        r4buf(ij)=arr(i,j)
       enddo
      enddo

! Send message
      count=HEADER_SIZE+Ni*Nj
      datatype=MPI_REAL
      dest=my_coupler_rank
      tag=generate_tag(102,my_rank_in_global,dataname)
      comm=MPI_COMM_myglobal

      if (VERB) then
       write(LogUnit,*) 'compsend_r4: calling MPI_Send dest=',dest
       write(LogUnit,*) 'compsend_r4: dataname=',dataname
       call flush(LogUnit)
      endif
      call MPI_Send( r4buf, count, datatype, dest, tag, comm, ierr )
      if (VERB) then
       write(LogUnit,*) 'compsend_r4: returned ierr=',ierr
       call flush(LogUnit)
      endif

      if (ierr.ne.0) then
       write(LogUnit,*) 'compsend_r4: rank(W,G,L)=',
     &            my_rank_in_world,my_rank_in_global,my_rank_in_local,
     &            ' ierr=',ierr
       stop 'compsend_r4: MPI_Send failed'
      endif

!     ------------------------------------------------------------------
      return
      end
!=======================================================================
1	C $Header: $
2	C $Name: $
3
4	!=======================================================================
5	subroutine compsend_r4( dataname, Ni,Oi,Nj,Oj, arr )
6	implicit none
7	! Arguments
8	character() dataname
9	integer Ni,Oi,Nj,Oj
10	real*4 arr(1-Oi:Ni+Oi,1-Oj:Nj+Oj)
11	! Predefined constants/arrays
12	#include "CPLR_SIG.h"
13	! MPI variables
14	#include "mpif.h"
15	integer count,datatype,dest,tag,comm,ierr
16	! Functions
17	integer generate_tag
18	! Local
19	integer i,j,ij,Io,Jo
20	! ------------------------------------------------------------------
21
22	if (HEADER_SIZE+Ni*Nj.gt.MAX_R4_BUFLEN)
23	& stop 'compsend_r4: Nx*Ny too big'
24
25	! Set up buffer
26	Io=my_tile_i0(1)
27	Jo=my_tile_j0(1)
28	r4buf(1)=float(Io)
29	r4buf(2)=float(Ni)
30	r4buf(3)=float(Jo)
31	r4buf(4)=float(Nj)
32	call mitcplr_char2real( dataname, r4buf(9) )
33
34	! Copy interior of tile to buffer
35	do j=1,Nj
36	do i=1,Ni
37	ij=HEADER_SIZE+i+Ni*(j-1)
38	r4buf(ij)=arr(i,j)
39	enddo
40	enddo
41
42	! Send message
43	count=HEADER_SIZE+Ni*Nj
44	datatype=MPI_REAL
45	dest=my_coupler_rank
46	tag=generate_tag(102,my_rank_in_global,dataname)
47	comm=MPI_COMM_myglobal
48
49	if (VERB) then
50	write(LogUnit,*) 'compsend_r4: calling MPI_Send dest=',dest
51	write(LogUnit,*) 'compsend_r4: dataname=',dataname
52	call flush(LogUnit)
53	endif
54	call MPI_Send( r4buf, count, datatype, dest, tag, comm, ierr )
55	if (VERB) then
56	write(LogUnit,*) 'compsend_r4: returned ierr=',ierr
57	call flush(LogUnit)
58	endif
59
60	if (ierr.ne.0) then
61	write(LogUnit,*) 'compsend_r4: rank(W,G,L)=',
62	& my_rank_in_world,my_rank_in_global,my_rank_in_local,
63	& ' ierr=',ierr
64	stop 'compsend_r4: MPI_Send failed'
65	endif
66
67	! ------------------------------------------------------------------
68	return
69	end
70	!=======================================================================