eesupp/src/exch_xyz_rx.template

C $Header: /u/gcmpack/MITgcm/eesupp/src/exch_xyz_rx.template,v 1.9 2007/08/17 18:34:17 jmc Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_EEOPTIONS.h"

CBOP
C     !ROUTINE: EXCH_XYZ_RX

C     !INTERFACE:
      SUBROUTINE EXCH_XYZ_RX(
     U                       phi,
     I                       myThid )

C     !DESCRIPTION:
C     *==========================================================*
C     | SUBROUTINE EXCH_XYZ_RX
C     | o Handle exchanges for _RX, three-dim scalar arrays.
C     *==========================================================*
C     | Invoke appropriate exchange routine depending on type
C     | of grid (cube or globally indexed) to be operated on.
C     *==========================================================*

C     !USES:
      IMPLICIT NONE

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

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     phi    :: Array with overlap regions are to be exchanged
C     myThid :: My thread id.
      _RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
      INTEGER myThid

C     !LOCAL VARIABLES:
#ifndef ALLOW_EXCH2
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
#endif
CEOP

#ifdef ALLOW_EXCH2
      CALL EXCH2_3D_RX( phi, Nr, myThid )
      RETURN
#else /* ALLOW_EXCH2 */

      OLw        = OLx
      OLe        = OLx
      OLn        = OLy
      OLs        = OLy
      exchWidthX = OLx
      exchWidthY = OLy
      myNz       = Nr
      IF (useCubedSphereExchange) THEN
       CALL EXCH1_RX_CUBE( phi, .FALSE.,
     I            OLw, OLe, OLs, OLn, myNz,
     I            exchWidthX, exchWidthY,
     I            EXCH_UPDATE_CORNERS, myThid )
      ELSE
       CALL EXCH1_RX( phi,
     I            OLw, OLe, OLs, OLn, myNz,
     I            exchWidthX, exchWidthY,
     I            EXCH_UPDATE_CORNERS, myThid )
      ENDIF

      RETURN
#endif /* ALLOW_EXCH2 */
      END
1	C $Header: /u/gcmpack/MITgcm/eesupp/src/exch_xyz_rx.template,v 1.9 2007/08/17 18:34:17 jmc Exp $
2	C $Name: $
3
4	#include "PACKAGES_CONFIG.h"
5	#include "CPP_EEOPTIONS.h"
6
7	CBOP
8	C !ROUTINE: EXCH_XYZ_RX
9
10	C !INTERFACE:
11	SUBROUTINE EXCH_XYZ_RX(
12	U phi,
13	I myThid )
14
15	C !DESCRIPTION:
16	C ==========================================================
17	C \| SUBROUTINE EXCH_XYZ_RX
18	C \| o Handle exchanges for _RX, three-dim scalar arrays.
19	C ==========================================================
20	C \| Invoke appropriate exchange routine depending on type
21	C \| of grid (cube or globally indexed) to be operated on.
22	C ==========================================================
23
24	C !USES:
25	IMPLICIT NONE
26
27	C === Global data ===
28	#include "SIZE.h"
29	#include "EEPARAMS.h"
30
31	C !INPUT/OUTPUT PARAMETERS:
32	C === Routine arguments ===
33	C phi :: Array with overlap regions are to be exchanged
34	C myThid :: My thread id.
35	_RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
36	INTEGER myThid
37
38	C !LOCAL VARIABLES:
39	#ifndef ALLOW_EXCH2
40	C == Local variables ==
41	C OL[wens] :: Overlap extents in west, east, north, south.
42	C exchWidth[XY] :: Extent of regions that will be exchanged.
43	INTEGER OLw, OLe, OLn, OLs, exchWidthX, exchWidthY, myNz
44	#endif
45	CEOP
46
47	#ifdef ALLOW_EXCH2
48	CALL EXCH2_3D_RX( phi, Nr, myThid )
49	RETURN
50	#else /* ALLOW_EXCH2 */
51
52	OLw = OLx
53	OLe = OLx
54	OLn = OLy
55	OLs = OLy
56	exchWidthX = OLx
57	exchWidthY = OLy
58	myNz = Nr
59	IF (useCubedSphereExchange) THEN
60	CALL EXCH1_RX_CUBE( phi, .FALSE.,
61	I OLw, OLe, OLs, OLn, myNz,
62	I exchWidthX, exchWidthY,
63	I EXCH_UPDATE_CORNERS, myThid )
64	ELSE
65	CALL EXCH1_RX( phi,
66	I OLw, OLe, OLs, OLn, myNz,
67	I exchWidthX, exchWidthY,
68	I EXCH_UPDATE_CORNERS, myThid )
69	ENDIF
70
71	RETURN
72	#endif /* ALLOW_EXCH2 */
73	END