eesupp/src/exch_z_xyz_rx.template

C $Header: /u/gcmpack/models/MITgcmUV/eesupp/src/Attic/exch_z_xyz_rx.template,v 1.1.2.2 2001/04/12 10:52:48 cnh Exp $
C $Name: pre38-close $

#include "CPP_EEOPTIONS.h"

      SUBROUTINE EXCH_Z_XYZ_RX(
     U                       phi,
     I                       myThid )
C     /==========================================================\
C     | SUBROUTINE EXCH_Z_XYZ_RX                                 |
C     | o Handle exchanges for _RX,      3-dimensional arrays.   |
C     |==========================================================|
C     | Do true shared-memory data transfers and "messaging"     |
C     | tranfers for blocking case of data transfers.            |
C     | Applications call this routine using                     |
C     |  CALL EXCH..( x, myThid )                                |
C     | where x is a two-dimensional array with overlaps.        |
C     | This routine does true-shared-memory copies for blocks   |
C     | within a thread. It will also do MPI meesaging between   |
C     | different processes.                                     |
C     | Note:                                                    |
C     | =====                                                    |
C     | This routine is a default mode interface to the more     |
C     | general exchange routine. Under default mode the         |
C     | exchange is for an edge width of OLx, and OLy and        |
C     | updates corner points as well as edges ( slightly more   |
C     | synchronisation required).                               |
C     \==========================================================/
      IMPLICIT NONE

C     === Global data ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "EESUPPORT.h"

C     === Routine arguments ===
C     phi    - Array with overlap regions are to be exchanged
C              Note - The interface to EXCH_RX assumes that
C              the standard Fortran 77 sequence association rules
C              apply.
C     myThid - My thread id.
      _RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
C     OL[wens]       - Overlap extents in west, east, north, south.
C     exchWidth[XY]  - Extent of regions that will be exchanged.
      INTEGER OLw, OLe, OLn, OLs, exchWidthX, exchWidthY, myNz

      OLw        = OLx
      OLe        = OLx
      OLn        = OLy
      OLs        = OLy
      exchWidthX = OLx
      exchWidthY = OLy
      myNz       = Nr
C     ** NOTE ** The exchange routine we use here does not
C                require the preceeding and following barriers.
C                However, the slow, simple exchange interface
C                that is calling it here is meant to ensure
C                that threads are synchronised before exchanges
C                begine.
      IF (useCubedSphereExchange) THEN
       CALL EXCH_Z_RX_CUBE( phi,
     I            OLw, OLe, OLs, OLn, myNz,
     I            exchWidthX, exchWidthY,
     I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
      ELSE
       CALL EXCH_RX( phi,
     I            OLw, OLe, OLs, OLn, myNz,
     I            exchWidthX, exchWidthY,
     I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
      ENDIF

      RETURN
      END
1	adcroft	1.2	C $Header: /u/gcmpack/models/MITgcmUV/eesupp/src/Attic/exch_z_xyz_rx.template,v 1.1.2.2 2001/04/12 10:52:48 cnh Exp $
2			C $Name: pre38-close $
3
4			#include "CPP_EEOPTIONS.h"
5
6			SUBROUTINE EXCH_Z_XYZ_RX(
7			U phi,
8			I myThid )
9			C /==========================================================\
10			C \| SUBROUTINE EXCH_Z_XYZ_RX \|
11			C \| o Handle exchanges for _RX, 3-dimensional arrays. \|
12			C \|==========================================================\|
13			C \| Do true shared-memory data transfers and "messaging" \|
14			C \| tranfers for blocking case of data transfers. \|
15			C \| Applications call this routine using \|
16			C \| CALL EXCH..( x, myThid ) \|
17			C \| where x is a two-dimensional array with overlaps. \|
18			C \| This routine does true-shared-memory copies for blocks \|
19			C \| within a thread. It will also do MPI meesaging between \|
20			C \| different processes. \|
21			C \| Note: \|
22			C \| ===== \|
23			C \| This routine is a default mode interface to the more \|
24			C \| general exchange routine. Under default mode the \|
25			C \| exchange is for an edge width of OLx, and OLy and \|
26			C \| updates corner points as well as edges ( slightly more \|
27			C \| synchronisation required). \|
28			C \==========================================================/
29			IMPLICIT NONE
30
31			C === Global data ===
32			#include "SIZE.h"
33			#include "EEPARAMS.h"
34			#include "EESUPPORT.h"
35
36			C === Routine arguments ===
37			C phi - Array with overlap regions are to be exchanged
38			C Note - The interface to EXCH_RX assumes that
39			C the standard Fortran 77 sequence association rules
40			C apply.
41			C myThid - My thread id.
42			_RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43			INTEGER myThid
44			CEndOfInterface
45
46			C == Local variables ==
47			C OL[wens] - Overlap extents in west, east, north, south.
48			C exchWidth[XY] - Extent of regions that will be exchanged.
49			INTEGER OLw, OLe, OLn, OLs, exchWidthX, exchWidthY, myNz
50
51			OLw = OLx
52			OLe = OLx
53			OLn = OLy
54			OLs = OLy
55			exchWidthX = OLx
56			exchWidthY = OLy
57			myNz = Nr
58			C NOTE The exchange routine we use here does not
59			C require the preceeding and following barriers.
60			C However, the slow, simple exchange interface
61			C that is calling it here is meant to ensure
62			C that threads are synchronised before exchanges
63			C begine.
64			IF (useCubedSphereExchange) THEN
65			CALL EXCH_Z_RX_CUBE( phi,
66			I OLw, OLe, OLs, OLn, myNz,
67			I exchWidthX, exchWidthY,
68			I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
69			ELSE
70			CALL EXCH_RX( phi,
71			I OLw, OLe, OLs, OLn, myNz,
72			I exchWidthX, exchWidthY,
73			I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
74			ENDIF
75
76			RETURN
77			END