/[MITgcm]/MITgcm/pkg/compon_communic/comprecv_r8.F
ViewVC logotype

Annotation of /MITgcm/pkg/compon_communic/comprecv_r8.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Mon Dec 15 02:35:29 2003 UTC (20 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint58e_post, checkpoint57v_post, checkpoint52n_post, checkpoint52j_post, checkpoint53d_post, checkpoint58u_post, checkpoint58w_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint52e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint52l_post, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint52k_post, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint58q_post, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59h, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, checkpoint52d_post, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint53g_post, checkpoint52f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, hrcube5, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint55e_post, checkpoint58k_post, checkpoint52i_post, checkpoint52j_pre, checkpoint58v_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint58s_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, checkpoint52h_pre, checkpoint52f_pre, checkpoint57h_post, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint55d_post
coupler - multi-components interface: low-level communication S/R package.

1 jmc 1.1 !=======================================================================
2     subroutine comprecv_r8( dataname, Ni,Oi,Nj,Oj, arr )
3     implicit none
4     ! Arguments
5     character*(*) dataname
6     integer Ni,Oi,Nj,Oj,Io,Jo
7     real*8 arr(1-Oi:Ni+Oi,1-Oj:Nj+Oj)
8     ! Predefined constants/arrays
9     #include "CPLR_SIG.h"
10     ! MPI variables
11     #include "mpif.h"
12     integer count,dtype,rank,tag,comm,ierr
13     integer stat(MPI_STATUS_SIZE)
14     ! Functions
15     integer generate_tag
16     ! Local
17     integer i,j,ij,nx,ny
18     character*(MAXLEN_COMP_NAME) recvdname
19     ! ------------------------------------------------------------------
20    
21     if (HEADER_SIZE+Ni*Nj.gt.MAX_R8_BUFLEN)
22     & stop 'comprecv_r8: Nx*Ny too big'
23    
24     ! Receive message
25     count=HEADER_SIZE+MAX_R8_BUFLEN
26     dtype=MPI_DOUBLE_PRECISION
27     tag=generate_tag(121,my_rank_in_global,dataname)
28     rank=my_coupler_rank
29     comm=MPI_COMM_myglobal
30    
31     if (VERB) then
32     write(LogUnit,*) 'comprecv_r8: calling MPI_Recv rank=',rank
33     write(LogUnit,*) 'comprecv_r8: dataname=',dataname
34     call flush(LogUnit)
35     endif
36     call MPI_Recv(r8buf, count, dtype, rank, tag, comm, stat, ierr)
37     if (VERB) then
38     write(LogUnit,*) 'comprecv_r8: returned ierr=',ierr
39     call flush(LogUnit)
40     endif
41    
42     if (ierr.ne.0) then
43     write(LogUnit,*) 'comprecv_r8tiles: rank(W,G)=',
44     & my_rank_in_world,my_rank_in_global,
45     & ' ierr=',ierr
46     stop 'comprecv_r8: MPI_Recv failed'
47     endif
48    
49     ! Extract buffer
50     Io=int(0.5+r8buf(1))
51     Jo=int(0.5+r8buf(2))
52     nx=int(0.5+r8buf(3))
53     ny=int(0.5+r8buf(4))
54     call mitcplr_real2char( r8buf(9), recvdname )
55    
56     if (Io.ne.my_tile_i0(1)) stop 'comprecv_r8: bad Io'
57     if (Jo.ne.my_tile_j0(1)) stop 'comprecv_r8: bad Jo'
58     if (nx.ne.my_tile_nx(1)) stop 'comprecv_r8: bad nx'
59     if (ny.ne.my_tile_ny(1)) stop 'comprecv_r8: bad ny'
60     if (recvdname .ne. dataname) then
61     write(LogUnit,*) 'comprecv_r8: recvdname = ',recvdname
62     write(LogUnit,*) 'comprecv_r8: dataname = ',dataname
63     stop 'comprecv_r8: recvdname != dataname'
64     endif
65    
66     ! Copy buffer to interior of tile
67     do j=1,Nj
68     do i=1,Ni
69     ij=HEADER_SIZE+i+Ni*(j-1)
70     arr(i,j)=r8buf(ij)
71     enddo
72     enddo
73    
74     ! ------------------------------------------------------------------
75     return
76     end
77     !=======================================================================

  ViewVC Help
Powered by ViewVC 1.1.22