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

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

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


Revision 1.1 - (show annotations) (download)
Mon Dec 15 02:35:29 2003 UTC (20 years, 5 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 !=======================================================================
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