eesupp/src/exch_z_3d_rx.template

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

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

CBOP
C     !ROUTINE: EXCH_Z_3D_RX

C     !INTERFACE:
      SUBROUTINE EXCH_Z_3D_RX(
     U                         phi,
     I                         myNz, myThid )
      IMPLICIT NONE
C     !DESCRIPTION:
C     *==========================================================*
C     | SUBROUTINE EXCH_Z_3D_RX
C     | o Handle exchanges for _RX, 3-dim, zeta coord arrays.
C     *==========================================================*
C     | Branch to appropriate exchange routine for a zeta coord
C     | array on the current grid scheme.
C     *==========================================================*

C     !USES:
C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"

C     !INPUT/OUTPUT PARAMETERS:
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     myNz   :: 3rd dimension of array to exchange
C     myThid :: My thread id.
      INTEGER myNz
      _RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,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
#endif
CEOP

#ifdef ALLOW_EXCH2
      CALL EXCH2_Z_3D_RX(
     U                    phi,
     I                    myNz, myThid )
      RETURN
#else /* ALLOW_EXCH2 */

      OLw        = OLx
      OLe        = OLx
      OLn        = OLy
      OLs        = OLy
      exchWidthX = OLx
      exchWidthY = OLy
      IF (useCubedSphereExchange) THEN
       CALL EXCH1_Z_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_z_3d_rx.template,v 1.3 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_Z_3D_RX
9
10	C !INTERFACE:
11	SUBROUTINE EXCH_Z_3D_RX(
12	U phi,
13	I myNz, myThid )
14	IMPLICIT NONE
15	C !DESCRIPTION:
16	C ==========================================================
17	C \| SUBROUTINE EXCH_Z_3D_RX
18	C \| o Handle exchanges for _RX, 3-dim, zeta coord arrays.
19	C ==========================================================
20	C \| Branch to appropriate exchange routine for a zeta coord
21	C \| array on the current grid scheme.
22	C ==========================================================
23
24	C !USES:
25	C == Global variables ==
26	#include "SIZE.h"
27	#include "EEPARAMS.h"
28
29	C !INPUT/OUTPUT PARAMETERS:
30	C == Routine arguments ==
31	C phi :: Array with overlap regions are to be exchanged
32	C Note - The interface to EXCH_RX assumes that
33	C the standard Fortran 77 sequence association rules
34	C apply.
35	C myNz :: 3rd dimension of array to exchange
36	C myThid :: My thread id.
37	INTEGER myNz
38	_RX phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy)
39	INTEGER myThid
40
41	C !LOCAL VARIABLES:
42	#ifndef ALLOW_EXCH2
43	C == Local variables ==
44	C OL[wens] :: Overlap extents in west, east, north, south.
45	C exchWidth[XY] :: Extent of regions that will be exchanged.
46	INTEGER OLw, OLe, OLn, OLs, exchWidthX, exchWidthY
47	#endif
48	CEOP
49
50	#ifdef ALLOW_EXCH2
51	CALL EXCH2_Z_3D_RX(
52	U phi,
53	I myNz, myThid )
54	RETURN
55	#else /* ALLOW_EXCH2 */
56
57	OLw = OLx
58	OLe = OLx
59	OLn = OLy
60	OLs = OLy
61	exchWidthX = OLx
62	exchWidthY = OLy
63	IF (useCubedSphereExchange) THEN
64	CALL EXCH1_Z_RX_CUBE( phi, .FALSE.,
65	I OLw, OLe, OLs, OLn, myNz,
66	I exchWidthX, exchWidthY,
67	I EXCH_UPDATE_CORNERS, myThid )
68	ELSE
69	CALL EXCH1_RX( phi,
70	I OLw, OLe, OLs, OLn, myNz,
71	I exchWidthX, exchWidthY,
72	I EXCH_UPDATE_CORNERS, myThid )
73	ENDIF
74
75	RETURN
76	#endif /* ALLOW_EXCH2 */
77	END