pkg/compon_communic/coupsend_r8.F

C $Header: /u/gcmpack/MITgcm/pkg/compon_communic/coupsend_r8.F,v 1.2 2007/10/08 23:58:20 jmc Exp $
C $Name:  $

!=======================================================================
      subroutine coupsend_r8( component, dataname, Nx, Ny, arr )
      implicit none
! Arguments
      character*(*) component
      character*(*) dataname
      integer Nx,Ny
      real*8 arr(Nx,Ny)
! Predefined constants/arrays
#include "CPLR_SIG.h"
! MPI variables
#include "mpif.h"
      integer count,dtype,dest,tag,comm,ierr
! Functions
      integer mitcplr_match_comp
      integer generate_tag
! Local
      integer compind,numprocs
      integer i,j,ij,n
      integer Ni,Io,Nj,Jo
!     ------------------------------------------------------------------

! Establish who I am communicating with
      compind=mitcplr_match_comp( component )
      if (compind.le.0) stop 'coupsend_r8: Bad component id'
      comm=MPI_COMM_compcplr( compind )
      numprocs=num_component_procs(compind)
      if (numprocs.lt.1) then
       write(LogUnit,*) 'coupsend_r8: compind = ',compind
       stop 'coupsend_r8: numprocs < 1'
      endif
      if (VERB)
     &  write(LogUnit,*) 'coupsend_r8: ',component_Name(compind)
      if (VERB)
     &  write(LogUnit,*) 'coupsend_r8: dataname=',dataname

! Foreach component process
      do n=1,numprocs

! Create header
       Io=component_tile_i0(1,n,compind)
       Jo=component_tile_j0(1,n,compind)
       Ni=component_tile_nx(1,n,compind)
       Nj=component_tile_ny(1,n,compind)
       r8buf(1)=float( Io )
       r8buf(2)=float( Jo )
       r8buf(3)=float( Ni )
       r8buf(4)=float( Nj )
       call mitcplr_char2dbl( dataname, r8buf(9) )

! Pack data
       do j=1,Nj
        do i=1,Ni
         ij=HEADER_SIZE+i+Ni*(j-1)
         r8buf(ij)=arr(Io+i-1,Jo+j-1)
        enddo
       enddo

! Send message
       count=HEADER_SIZE+Ni*Nj
       dtype=MPI_DOUBLE_PRECISION
       tag=generate_tag(121,n,dataname)
       dest=rank_component_procs(n,compind)

       if (VERB) then
        write(LogUnit,*) 'coupsend_r8: calling MPI_Send dest=',dest,
     &    ' proc=',n,'/',numprocs
        call flush(LogUnit)
       endif
       call MPI_Send( r8buf, count, dtype, dest, tag, comm, ierr )
       if (VERB) then
        write(LogUnit,*) 'coupsend_r8: returned ierr=',ierr
        call flush(LogUnit)
       endif

       if (ierr.ne.0) then
        write(LogUnit,*) 'coupsend_r8tiles: rank(W,G)=',
     &            my_rank_in_world,my_rank_in_global,
     &            ' ierr=',ierr
        stop 'coupsend_r8: MPI_Recv failed'
       endif

      enddo ! n

!     ------------------------------------------------------------------
      return
      end
!=======================================================================
1	C $Header: /u/gcmpack/MITgcm/pkg/compon_communic/coupsend_r8.F,v 1.2 2007/10/08 23:58:20 jmc Exp $
2	C $Name: $
3
4	!=======================================================================
5	subroutine coupsend_r8( component, dataname, Nx, Ny, arr )
6	implicit none
7	! Arguments
8	character() component
9	character() dataname
10	integer Nx,Ny
11	real*8 arr(Nx,Ny)
12	! Predefined constants/arrays
13	#include "CPLR_SIG.h"
14	! MPI variables
15	#include "mpif.h"
16	integer count,dtype,dest,tag,comm,ierr
17	! Functions
18	integer mitcplr_match_comp
19	integer generate_tag
20	! Local
21	integer compind,numprocs
22	integer i,j,ij,n
23	integer Ni,Io,Nj,Jo
24	! ------------------------------------------------------------------
25
26	! Establish who I am communicating with
27	compind=mitcplr_match_comp( component )
28	if (compind.le.0) stop 'coupsend_r8: Bad component id'
29	comm=MPI_COMM_compcplr( compind )
30	numprocs=num_component_procs(compind)
31	if (numprocs.lt.1) then
32	write(LogUnit,*) 'coupsend_r8: compind = ',compind
33	stop 'coupsend_r8: numprocs < 1'
34	endif
35	if (VERB)
36	& write(LogUnit,*) 'coupsend_r8: ',component_Name(compind)
37	if (VERB)
38	& write(LogUnit,*) 'coupsend_r8: dataname=',dataname
39
40	! Foreach component process
41	do n=1,numprocs
42
43	! Create header
44	Io=component_tile_i0(1,n,compind)
45	Jo=component_tile_j0(1,n,compind)
46	Ni=component_tile_nx(1,n,compind)
47	Nj=component_tile_ny(1,n,compind)
48	r8buf(1)=float( Io )
49	r8buf(2)=float( Jo )
50	r8buf(3)=float( Ni )
51	r8buf(4)=float( Nj )
52	call mitcplr_char2dbl( dataname, r8buf(9) )
53
54	! Pack data
55	do j=1,Nj
56	do i=1,Ni
57	ij=HEADER_SIZE+i+Ni*(j-1)
58	r8buf(ij)=arr(Io+i-1,Jo+j-1)
59	enddo
60	enddo
61
62	! Send message
63	count=HEADER_SIZE+Ni*Nj
64	dtype=MPI_DOUBLE_PRECISION
65	tag=generate_tag(121,n,dataname)
66	dest=rank_component_procs(n,compind)
67
68	if (VERB) then
69	write(LogUnit,*) 'coupsend_r8: calling MPI_Send dest=',dest,
70	& ' proc=',n,'/',numprocs
71	call flush(LogUnit)
72	endif
73	call MPI_Send( r8buf, count, dtype, dest, tag, comm, ierr )
74	if (VERB) then
75	write(LogUnit,*) 'coupsend_r8: returned ierr=',ierr
76	call flush(LogUnit)
77	endif
78
79	if (ierr.ne.0) then
80	write(LogUnit,*) 'coupsend_r8tiles: rank(W,G)=',
81	& my_rank_in_world,my_rank_in_global,
82	& ' ierr=',ierr
83	stop 'coupsend_r8: MPI_Recv failed'
84	endif
85
86	enddo ! n
87
88	! ------------------------------------------------------------------
89	return
90	end
91	!=======================================================================