/[MITgcm]/MITgcm/eesupp/src/exch_rx_recv_get_x.template
ViewVC logotype

Diff of /MITgcm/eesupp/src/exch_rx_recv_get_x.template

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

revision 1.1 by adcroft, Tue May 29 14:06:38 2001 UTC revision 1.14 by jmc, Mon May 17 02:28:06 2010 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
4    #undef EXCH_USE_SPINNING
5    
6    CBOP
7    C     !ROUTINE: EXCH_RX_RECV_GET_X
8    
9    C     !INTERFACE:
10        SUBROUTINE EXCH_RX_RECV_GET_X( array,        SUBROUTINE EXCH_RX_RECV_GET_X( array,
11       I            myOLw, myOLe, myOLs, myOLn, myNz,       I            myOLw, myOLe, myOLs, myOLn, myNz,
12       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
13       I            theSimulationMode, theCornerMode, myThid )       I            theSimulationMode, theCornerMode, myThid )
 C     /==========================================================\  
 C     | SUBROUTINE RECV_RX_GET_X                                 |  
 C     | o "Send" or "put" X edges for RX array.                  |  
 C     |==========================================================|  
 C     | Routine that invokes actual message passing send or      |  
 C     | direct "put" of data to update X faces of an XY[R] array.|  
 C     \==========================================================/  
14        IMPLICIT NONE        IMPLICIT NONE
15    
16    C     !DESCRIPTION:
17    C     *==========================================================*
18    C     | SUBROUTINE RECV_RX_GET_X
19    C     | o "Send" or "put" X edges for RX array.
20    C     *==========================================================*
21    C     | Routine that invokes actual message passing send or
22    C     | direct "put" of data to update X faces of an XY[R] array.
23    C     *==========================================================*
24    
25    C     !USES:
26  C     == Global variables ==  C     == Global variables ==
27  #include "SIZE.h"  #include "SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29  #include "EESUPPORT.h"  #include "EESUPPORT.h"
30  #include "EXCH.h"  #include "EXCH.h"
31    
32    C     !INPUT/OUTPUT PARAMETERS:
33  C     == Routine arguments ==  C     == Routine arguments ==
34  C     array - Array with edges to exchange.  C     array :: Array with edges to exchange.
35  C     myOLw - West, East, North and South overlap region sizes.  C     myOLw :: West, East, North and South overlap region sizes.
36  C     myOLe  C     myOLe
37  C     myOLn  C     myOLn
38  C     myOLs  C     myOLs
39  C     exchWidthX - Width of data region exchanged.  C     exchWidthX :: Width of data region exchanged.
40  C     exchWidthY  C     exchWidthY
41  C     theSimulationMode - Forward or reverse mode exchange ( provides  C     theSimulationMode :: Forward or reverse mode exchange ( provides
42  C                         support for adjoint integration of code. )  C                         support for adjoint integration of code. )
43  C     theCornerMode     - Flag indicating whether corner updates are  C     theCornerMode     :: Flag indicating whether corner updates are
44  C                         needed.  C                         needed.
45  C     myThid            - Thread number of this instance of S/R EXCH...  C     myThid            :: Thread number of this instance of S/R EXCH...
46  C     eBl               - Edge buffer level  C     eBl               :: Edge buffer level
47        INTEGER myOLw        INTEGER myOLw
48        INTEGER myOLe        INTEGER myOLe
49        INTEGER myOLs        INTEGER myOLs
# Line 48  C     eBl               - Edge buffer le Line 57  C     eBl               - Edge buffer le
57        INTEGER theSimulationMode        INTEGER theSimulationMode
58        INTEGER theCornerMode        INTEGER theCornerMode
59        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
60    
61    C     !LOCAL VARIABLES:
62  C     == Local variables ==  C     == Local variables ==
63  C     I, J, K, iMin, iMax, iB    - Loop counters and extents  C     i, j, k, iMin, iMax, iB    :: Loop counters and extents
64  C     bi, bj  C     bi, bj
65  C     biW, bjW                   - West tile indices  C     biW, bjW                   :: West tile indices
66  C     biE, bjE                   - East tile indices  C     biE, bjE                   :: East tile indices
67  C     eBl                        - Current exchange buffer level  C     eBl                        :: Current exchange buffer level
68  C     theProc, theTag, theType,  - Variables used in message building  C     theProc, theTag, theType,  :: Variables used in message building
69  C     theSize  C     theSize
70  C     westCommMode               - Working variables holding type  C     westCommMode               :: Working variables holding type
71  C     eastCommMode                 of communication a particular  C     eastCommMode                  of communication a particular
72  C                                  tile face uses.  C                                   tile face uses.
73        INTEGER I, J, K, iMin, iMax, iB, iB0        INTEGER i, j, k, iMin, iMax, iB, iB0
74        INTEGER bi, bj, biW, bjW, biE, bjE        INTEGER bi, bj, biW, bjW, biE, bjE
75        INTEGER eBl        INTEGER eBl
76        INTEGER westCommMode        INTEGER westCommMode
77        INTEGER eastCommMode        INTEGER eastCommMode
78    #ifdef EXCH_USE_SPINNING
79        INTEGER spinCount        INTEGER spinCount
80    #endif
81  #ifdef ALLOW_USE_MPI  #ifdef ALLOW_USE_MPI
82        INTEGER theProc, theTag, theType, theSize        INTEGER theProc, theTag, theType, theSize, pReqI
83        INTEGER mpiStatus(MPI_STATUS_SIZE,4), mpiRc        INTEGER mpiStatus(MPI_STATUS_SIZE,4), mpiRc
84  #endif  #endif
85    CEOP
86    
87    C--   Under a "put" scenario we
 C--   Under a "put" scenario we  
88  C--     i. set completetion signal for buffer we put into.  C--     i. set completetion signal for buffer we put into.
89  C--    ii. wait for completetion signal indicating data has been put in  C--    ii. wait for completetion signal indicating data has been put in
90  C--        our buffer.  C--        our buffer.
91  C--   Under a messaging mode we "receive" the message.  C--   Under a messaging mode we "receive" the message.
92  C--   Under a "get" scenario we  C--   Under a "get" scenario we
# Line 83  C--     i. Check that the data is ready. Line 94  C--     i. Check that the data is ready.
94  C--    ii. Read the data.  C--    ii. Read the data.
95  C--   iii. Set data read flag + memory sync.  C--   iii. Set data read flag + memory sync.
96    
   
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         ebL = exchangeBufLevel(1,bi,bj)  
         westCommMode  = _tileCommModeW(bi,bj)  
         eastCommMode  = _tileCommModeE(bi,bj)  
         biE =  _tileBiE(bi,bj)  
         bjE =  _tileBjE(bi,bj)  
         biW =  _tileBiW(bi,bj)  
         bjW =  _tileBjW(bi,bj)  
         IF ( westCommMode .EQ. COMM_MSG ) THEN  
97  #ifdef ALLOW_USE_MPI  #ifdef ALLOW_USE_MPI
98  #ifndef ALWAYS_USE_MPI  #ifndef ALWAYS_USE_MPI
99           IF ( usingMPI ) THEN        IF ( usingMPI ) THEN
100  #endif  #endif
101    C--   Receive buffer data: Only Master Thread do proc communication
102          _BEGIN_MASTER(myThid)
103    
104          DO bj=1,nSy
105           DO bi=1,nSx
106            eBl = exchangeBufLevel(1,bi,bj)
107            westCommMode = _tileCommModeW(bi,bj)
108            eastCommMode = _tileCommModeE(bi,bj)
109            biE = _tileBiE(bi,bj)
110            bjE = _tileBjE(bi,bj)
111            biW = _tileBiW(bi,bj)
112            bjW = _tileBjW(bi,bj)
113            theType = _MPI_TYPE_RX
114            theSize = sNy*exchWidthX*myNz
115    
116            IF ( westCommMode .EQ. COMM_MSG ) THEN
117           theProc = tilePidW(bi,bj)           theProc = tilePidW(bi,bj)
118           theTag  = _tileTagRecvW(bi,bj)           theTag  = _tileTagRecvW(bi,bj)
119           theType = MPI_DOUBLE_PRECISION  # ifndef ALLOW_AUTODIFF_OPENAD_AMPI
120           theSize = sNy*exchWidthX*myNz           CALL MPI_Recv( westRecvBuf_RX(1,eBl,bi,bj), theSize,
121           CALL MPI_Recv( westRecvBuf_RX(1,eBl,bi,bj), theSize, theType,       &                  theType, theProc, theTag, MPI_COMM_MODEL,
      &                  theProc, theTag, MPI_COMM_MODEL,  
122       &                  mpiStatus, mpiRc )       &                  mpiStatus, mpiRc )
123  #ifndef ALWAYS_USE_MPI  # else
124          ENDIF                           pReqI=exchNReqsX(1,bi,bj)+1
125  #endif           CALL ampi_recv_RX(
126  #endif /* ALLOW_USE_MPI */       &        westRecvBuf_RX(1,eBl,bi,bj) ,
127         &        theSize ,
128         &        theType ,
129         &        theProc ,
130         &        theTag ,
131         &        MPI_COMM_MODEL ,
132         &        exchReqIdX(pReqI,1,bi,bj),
133         &        exchNReqsX(1,bi,bj),
134         &        mpiStatus ,
135         &        mpiRc )
136    # endif /* ALLOW_AUTODIFF_OPENAD_AMPI */
137             westRecvAck(eBl,bi,bj) = 1
138          ENDIF          ENDIF
139    
140          IF ( eastCommMode .EQ. COMM_MSG ) THEN          IF ( eastCommMode .EQ. COMM_MSG ) THEN
 #ifdef ALLOW_USE_MPI  
 #ifndef ALWAYS_USE_MPI  
          IF ( usingMPI ) THEN  
 #endif  
141           theProc = tilePidE(bi,bj)           theProc = tilePidE(bi,bj)
142           theTag  = _tileTagRecvE(bi,bj)           theTag  = _tileTagRecvE(bi,bj)
143           theType = MPI_DOUBLE_PRECISION  # ifndef ALLOW_AUTODIFF_OPENAD_AMPI
144           theSize = sNy*exchWidthX*myNz           CALL MPI_Recv( eastRecvBuf_RX(1,eBl,bi,bj), theSize,
145           CALL MPI_Recv( eastRecvBuf_RX(1,eBl,bi,bj), theSize, theType,       &                  theType, theProc, theTag, MPI_COMM_MODEL,
      &                  theProc, theTag, MPI_COMM_MODEL,  
146       &                  mpiStatus, mpiRc )       &                  mpiStatus, mpiRc )
147  #ifndef ALWAYS_USE_MPI  # else
148          ENDIF                           pReqI=exchNReqsX(1,bi,bj)+1
149  #endif           CALL ampi_recv_RX(
150  #endif /* ALLOW_USE_MPI */       &        eastRecvBuf_RX(1,eBl,bi,bj) ,
151         &        theSize ,
152         &        theType ,
153         &        theProc ,
154         &        theTag ,
155         &        MPI_COMM_MODEL ,
156         &        exchReqIdX(pReqI,1,bi,bj),
157         &        exchNReqsX(1,bi,bj),
158         &        mpiStatus ,
159         &        mpiRc )
160    # endif /* ALLOW_AUTODIFF_OPENAD_AMPI */
161             eastRecvAck(eBl,bi,bj) = 1
162          ENDIF          ENDIF
163         ENDDO         ENDDO
164        ENDDO        ENDDO
165    
166  C--   Wait for buffers I am going read to be ready.  C--   Processes wait for buffers I am going to read to be ready.
167          IF ( .NOT.exchUsesBarrier  ) THEN
168           DO bj=1,nSy
169            DO bi=1,nSx
170             IF ( exchNReqsX(1,bi,bj) .GT. 0 ) THEN
171    # ifndef ALLOW_AUTODIFF_OPENAD_AMPI
172              CALL MPI_Waitall( exchNReqsX(1,bi,bj), exchReqIdX(1,1,bi,bj),
173         &                      mpiStatus, mpiRC )
174    # else
175              CALL ampi_waitall(
176         &         exchNReqsX(1,bi,bj),
177         &         exchReqIdX(1,1,bi,bj),
178         &         mpiStatus,
179         &         mpiRC )
180    # endif /* ALLOW_AUTODIFF_OPENAD_AMPI */
181             ENDIF
182    C        Clear outstanding requests counter
183             exchNReqsX(1,bi,bj) = 0
184            ENDDO
185           ENDDO
186          ENDIF
187    
188          _END_MASTER(myThid)
189    C--   need to sync threads after master has received data ;
190    C     (done after mpi waitall in case waitall is really needed)
191          _BARRIER
192    
193    #ifndef ALWAYS_USE_MPI
194          ENDIF
195    #endif
196    #endif /* ALLOW_USE_MPI */
197    #
198    C--   Threads wait for buffers I am going to read to be ready.
199    C     note: added BARRIER in exch_send_put S/R and here above (message mode)
200    C           so that we no longer needs this (undef EXCH_USE_SPINNING)
201    #ifdef EXCH_USE_SPINNING
202        IF ( exchUsesBarrier  ) THEN        IF ( exchUsesBarrier  ) THEN
203  C      o On some machines ( T90 ) use system barrier rather than spinning.  C      o On some machines ( T90 ) use system barrier rather than spinning.
204         CALL BARRIER( myThid )         CALL BARRIER( myThid )
# Line 139  C      o Spin waiting for completetion f Line 207  C      o Spin waiting for completetion f
207  C        i.e. we only lock waiting for data that we need.  C        i.e. we only lock waiting for data that we need.
208         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
209          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
210    
211           spinCount = 0           spinCount = 0
212           ebL = exchangeBufLevel(1,bi,bj)           eBl = exchangeBufLevel(1,bi,bj)
213           westCommMode = _tileCommModeW(bi,bj)           westCommMode = _tileCommModeW(bi,bj)
214           eastCommMode = _tileCommModeE(bi,bj)           eastCommMode = _tileCommModeE(bi,bj)
215    # ifndef ALLOW_AUTODIFF_OPENAD_AMPI
216     10    CONTINUE     10    CONTINUE
217            CALL FOOL_THE_COMPILER            CALL FOOL_THE_COMPILER( spinCount )
218            spinCount = spinCount+1            spinCount = spinCount+1
219  C         IF ( myThid .EQ. 1 .AND. spinCount .GT. _EXCH_SPIN_LIMIT ) THEN  C         IF ( myThid .EQ. 1 .AND. spinCount .GT. _EXCH_SPIN_LIMIT ) THEN
220  C          WRITE(*,*) ' eBl = ', ebl  C          WRITE(*,*) ' eBl = ', ebl
221  C          STOP ' S/R EXCH_RECV_GET_X: spinCount .GT. _EXCH_SPIN_LIMIT'  C          STOP ' S/R EXCH_RECV_GET_X: spinCount .GT. _EXCH_SPIN_LIMIT'
222  C         ENDIF  C         ENDIF
223            IF ( westRecvAck(eBl,bi,bj) .EQ. 0. ) GOTO 10            IF ( westRecvAck(eBl,bi,bj) .EQ. 0 ) GOTO 10
224            IF ( eastRecvAck(eBl,bi,bj) .EQ. 0. ) GOTO 10            IF ( eastRecvAck(eBl,bi,bj) .EQ. 0 ) GOTO 10
225    # else
226             DO WHILE ((westRecvAck(eBl,bi,bj) .EQ. 0
227         &             .OR.
228         &              eastRecvAck(eBl,bi,bj) .EQ. 0 ))
229              CALL FOOL_THE_COMPILER( spinCount )
230              spinCount = spinCount+1
231             ENDDO
232    # endif /* ALLOW_AUTODIFF_OPENAD_AMPI */
233  C        Clear outstanding requests  C        Clear outstanding requests
234           westRecvAck(eBl,bi,bj) = 0.           westRecvAck(eBl,bi,bj) = 0
235           eastRecvAck(eBl,bi,bj) = 0.           eastRecvAck(eBl,bi,bj) = 0
   
          IF ( exchNReqsX(1,bi,bj) .GT. 0 ) THEN  
 #ifdef ALLOW_USE_MPI  
 #ifndef ALWAYS_USE_MPI  
          IF ( usingMPI ) THEN  
 #endif  
           CALL MPI_Waitall( exchNReqsX(1,bi,bj), exchReqIdX(1,1,bi,bj),  
      &                      mpiStatus, mpiRC )  
 #ifndef ALWAYS_USE_MPI  
         ENDIF                  
 #endif  
 #endif /* ALLOW_USE_MPI */  
          ENDIF  
 C        Clear outstanding requests counter  
          exchNReqsX(1,bi,bj) = 0  
236  C        Update statistics  C        Update statistics
237           IF ( exchCollectStatistics ) THEN           IF ( exchCollectStatistics ) THEN
238            exchRecvXExchCount(1,bi,bj) = exchRecvXExchCount(1,bi,bj)+1            exchRecvXExchCount(1,bi,bj) = exchRecvXExchCount(1,bi,bj)+1
239            exchRecvXSpinCount(1,bi,bj) =            exchRecvXSpinCount(1,bi,bj) =
240       &    exchRecvXSpinCount(1,bi,bj)+spinCount       &    exchRecvXSpinCount(1,bi,bj)+spinCount
241            exchRecvXSpinMax(1,bi,bj) =            exchRecvXSpinMax(1,bi,bj) =
242       &    MAX(exchRecvXSpinMax(1,bi,bj),spinCount)       &    MAX(exchRecvXSpinMax(1,bi,bj),spinCount)
243            exchRecvXSpinMin(1,bi,bj) =            exchRecvXSpinMin(1,bi,bj) =
244       &    MIN(exchRecvXSpinMin(1,bi,bj),spinCount)       &    MIN(exchRecvXSpinMin(1,bi,bj),spinCount)
245           ENDIF           ENDIF
246    
247          ENDDO          ENDDO
248         ENDDO         ENDDO
249        ENDIF        ENDIF
250    #endif /* EXCH_USE_SPINNING */
251    
252    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
253    
254  C--   Read from the buffers  C--   Read from the buffers
255        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
256         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
257    
258          ebL = exchangeBufLevel(1,bi,bj)          eBl = exchangeBufLevel(1,bi,bj)
259          biE =  _tileBiE(bi,bj)          biE = _tileBiE(bi,bj)
260          bjE =  _tileBjE(bi,bj)          bjE = _tileBjE(bi,bj)
261          biW =  _tileBiW(bi,bj)          biW = _tileBiW(bi,bj)
262          bjW =  _tileBjW(bi,bj)          bjW = _tileBjW(bi,bj)
263          westCommMode = _tileCommModeW(bi,bj)          westCommMode = _tileCommModeW(bi,bj)
264          eastCommMode = _tileCommModeE(bi,bj)          eastCommMode = _tileCommModeE(bi,bj)
265    
266          IF     ( _theSimulationMode .EQ. FORWARD_SIMULATION ) THEN          IF     ( _theSimulationMode .EQ. FORWARD_SIMULATION ) THEN
267           iMin = sNx+1           iMin = sNx+1
268           iMax = sNx+exchWidthX           iMax = sNx+exchWidthX
# Line 202  C--   Read from the buffers Line 270  C--   Read from the buffers
270           IF (     eastCommMode .EQ. COMM_PUT           IF (     eastCommMode .EQ. COMM_PUT
271       &       .OR. eastCommMode .EQ. COMM_MSG ) THEN       &       .OR. eastCommMode .EQ. COMM_MSG ) THEN
272            iB  = 0            iB  = 0
273            DO K=1,myNz            DO k=1,myNz
274             DO J=1,sNy             DO j=1,sNy
275              DO I=iMin,iMax              DO i=iMin,iMax
276               iB = iB + 1               iB = iB + 1
277               array(I,J,K,bi,bj) = eastRecvBuf_RX(iB,eBl,bi,bj)               array(i,j,k,bi,bj) = eastRecvBuf_RX(iB,eBl,bi,bj)
278              ENDDO              ENDDO
279             ENDDO             ENDDO
280            ENDDO            ENDDO
281           ELSEIF ( eastCommMode .EQ. COMM_GET ) THEN           ELSEIF ( eastCommMode .EQ. COMM_GET ) THEN
282            DO K=1,myNz            DO k=1,myNz
283             DO J=1,sNy             DO j=1,sNy
284              iB = iB0              iB = iB0
285              DO I=iMin,iMax              DO i=iMin,iMax
286               iB = iB+1               iB = iB+1
287               array(I,J,K,bi,bj) = array(iB,J,K,biE,bjE)               array(i,j,k,bi,bj) = array(iB,j,k,biE,bjE)
288              ENDDO              ENDDO
289             ENDDO             ENDDO
290            ENDDO            ENDDO
# Line 225  C--   Read from the buffers Line 293  C--   Read from the buffers
293           iMin = sNx-exchWidthX+1           iMin = sNx-exchWidthX+1
294           iMax = sNx           iMax = sNx
295           iB0  = 1-exchWidthX-1           iB0  = 1-exchWidthX-1
296           IF (     eastCommMode .EQ. COMM_PUT           IF (     eastCommMode .EQ. COMM_PUT
297       &       .OR. eastCommMode .EQ. COMM_MSG ) THEN       &       .OR. eastCommMode .EQ. COMM_MSG ) THEN
298            iB  = 0            iB  = 0
299            DO K=1,myNz            DO k=1,myNz
300             DO J=1,sNy             DO j=1,sNy
301              DO I=iMin,iMax              DO i=iMin,iMax
302               iB = iB + 1               iB = iB + 1
303               array(I,J,K,bi,bj) =               array(i,j,k,bi,bj) =
304       &       array(I,J,K,bi,bj)+eastRecvBuf_RX(iB,eBl,bi,bj)       &       array(i,j,k,bi,bj) + eastRecvBuf_RX(iB,eBl,bi,bj)
305              ENDDO              ENDDO
306             ENDDO             ENDDO
307            ENDDO            ENDDO
308           ELSEIF ( eastCommMode .EQ. COMM_GET ) THEN           ELSEIF ( eastCommMode .EQ. COMM_GET ) THEN
309            DO K=1,myNz            DO k=1,myNz
310             DO J=1,sNy             DO j=1,sNy
311              iB = iB0              iB = iB0
312              DO I=iMin,iMax              DO i=iMin,iMax
313               iB = iB+1               iB = iB+1
314               array(I,J,K,bi,bj) =               array(i,j,k,bi,bj) =
315       &       array(I,J,K,bi,bj)+array(iB,J,K,biE,bjE)       &       array(i,j,k,bi,bj) + array(iB,j,k,biE,bjE)
316                 array(iB,j,k,biE,bjE) = 0.0
317              ENDDO              ENDDO
318             ENDDO             ENDDO
319            ENDDO            ENDDO
320           ENDIF           ENDIF
321          ENDIF          ENDIF
322    
323          IF     ( _theSimulationMode .EQ. FORWARD_SIMULATION ) THEN          IF     ( _theSimulationMode .EQ. FORWARD_SIMULATION ) THEN
324           iMin = 1-exchWidthX           iMin = 1-exchWidthX
325           iMax = 0           iMax = 0
# Line 257  C--   Read from the buffers Line 327  C--   Read from the buffers
327           IF (      westCommMode .EQ. COMM_PUT           IF (      westCommMode .EQ. COMM_PUT
328       &        .OR. westCommMode .EQ. COMM_MSG ) THEN       &        .OR. westCommMode .EQ. COMM_MSG ) THEN
329            iB  = 0            iB  = 0
330            DO K=1,myNz            DO k=1,myNz
331             DO J=1,sNy             DO j=1,sNy
332              DO I=iMin,iMax              DO i=iMin,iMax
333               iB = iB + 1               iB = iB + 1
334               array(I,J,K,bi,bj) = westRecvBuf_RX(iB,eBl,bi,bj)               array(i,j,k,bi,bj) = westRecvBuf_RX(iB,eBl,bi,bj)
335              ENDDO              ENDDO
336             ENDDO             ENDDO
337            ENDDO            ENDDO
338           ELSEIF ( westCommMode .EQ. COMM_GET ) THEN           ELSEIF ( westCommMode .EQ. COMM_GET ) THEN
339            DO K=1,myNz            DO k=1,myNz
340             DO J=1,sNy             DO j=1,sNy
341              iB = iB0              iB = iB0
342              DO I=iMin,iMax              DO i=iMin,iMax
343               iB = iB+1               iB = iB+1
344               array(I,J,K,bi,bj) = array(iB,J,K,biW,bjW)               array(i,j,k,bi,bj) = array(iB,j,k,biW,bjW)
345              ENDDO              ENDDO
346             ENDDO             ENDDO
347            ENDDO            ENDDO
# Line 280  C--   Read from the buffers Line 350  C--   Read from the buffers
350           iMin = 1           iMin = 1
351           iMax = 1+exchWidthX-1           iMax = 1+exchWidthX-1
352           iB0  = sNx           iB0  = sNx
353           IF (      westCommMode .EQ. COMM_PUT           IF (     westCommMode .EQ. COMM_PUT
354       &        .OR. westCommMode .EQ. COMM_MSG ) THEN       &       .OR. westCommMode .EQ. COMM_MSG ) THEN
355            iB  = 0            iB  = 0
356            DO K=1,myNz            DO k=1,myNz
357             DO J=1,sNy             DO j=1,sNy
358              DO I=iMin,iMax              DO i=iMin,iMax
359               iB = iB + 1               iB = iB + 1
360               array(I,J,K,bi,bj) =               array(i,j,k,bi,bj) =
361       &       array(I,J,K,bi,bj)+westRecvBuf_RX(iB,eBl,bi,bj)       &       array(i,j,k,bi,bj) + westRecvBuf_RX(iB,eBl,bi,bj)
362              ENDDO              ENDDO
363             ENDDO             ENDDO
364            ENDDO            ENDDO
365           ELSEIF ( westCommMode .EQ. COMM_GET ) THEN           ELSEIF ( westCommMode .EQ. COMM_GET ) THEN
366            DO K=1,myNz            DO k=1,myNz
367             DO J=1,sNy             DO j=1,sNy
368              iB = iB0              iB = iB0
369              DO I=iMin,iMax              DO i=iMin,iMax
370               iB = iB+1               iB = iB+1
371               array(I,J,K,bi,bj) =               array(i,j,k,bi,bj) =
372       &       array(I,J,K,bi,bj)+array(iB,J,K,biW,bjW)       &       array(i,j,k,bi,bj) + array(iB,j,k,biW,bjW)
373                 array(iB,j,k,biW,bjW) = 0.0
374              ENDDO              ENDDO
375             ENDDO             ENDDO
376            ENDDO            ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22