eesupp/src/exch_uv_xy_rx.template

C $Header: /u/gcmpack/MITgcm/eesupp/src/exch_uv_xy_rx.template,v 1.8 2007/07/28 16:07:35 heimbach Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_EEOPTIONS.h"
#ifdef ALLOW_EXCH2
#include "W2_OPTIONS.h"
#endif

CBOP
C     !ROUTINE: EXCH_UV_XY_RX

C     !INTERFACE:
      SUBROUTINE EXCH_UV_XY_RX(
     U                       Uphi, Vphi,
     I                       withSigns, myThid )
      IMPLICIT NONE
C     !DESCRIPTION:
C     *==========================================================*
C     | SUBROUTINE EXCH_UV_XY_RX
C     | o Handle exchanges for _RX, two-dimensional arrays.
C     *==========================================================*
C     | Driver exchange routine which branches to cube sphere or
C     | global, simple cartesian index grid. Exchange routine is
C     | called with two arrays that are components of a vector.
C     | These components are rotated and interchanged on the
C     | rotated grid during cube exchanges.
C     *==========================================================*

C     !USES:
C     === Global data ===
#include "SIZE.h"
#include "EEPARAMS.h"
cph-mpi#include "EESUPPORT.h"
c#include "EXCH.h"

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     Uphi      :: Arrays with overlap regions are to be exchanged
C     Vphi         Note - The interface to EXCH_ assumes that
C                  the standard Fortran 77 sequence association rules
C                  apply.
C     myThid    :: My thread id.
C     withSigns :: Flag controlling whether vector is signed.
      _RX Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RX Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      LOGICAL withSigns
      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
#ifdef W2_USE_R1_ONLY
      CALL EXCH2_UV_CGRID_3D_RX(
     U                       Uphi, Vphi,
     I                       withSigns, 1, myThid )
#else
      CALL EXCH2_UV_XY_RX(
     U                       Uphi, Vphi,
     I                       withSigns, myThid )
#endif
      RETURN
#else /* ALLOW_EXCH2 */

      OLw        = OLx
      OLe        = OLx
      OLn        = OLy
      OLs        = OLy
      exchWidthX = OLx
      exchWidthY = OLy
      myNz       = 1
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_UV_RX_CUBE( Uphi, Vphi, withSigns,
     I            OLw, OLe, OLs, OLn, myNz,
     I            exchWidthX, exchWidthY,
     I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
      ELSE
c      CALL EXCH_RX( Uphi,
c    I            OLw, OLe, OLs, OLn, myNz,
c    I            exchWidthX, exchWidthY,
c    I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
c      CALL EXCH_RX( Vphi,
c    I            OLw, OLe, OLs, OLn, myNz,
c    I            exchWidthX, exchWidthY,
c    I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
c_jmc: for JAM compatibility, replace the 2 CALLs above by the 2 CPP_MACROs:
       _EXCH_XY_RX( Uphi, myThid )
       _EXCH_XY_RX( Vphi, myThid )
      ENDIF

      RETURN
#endif /* ALLOW_EXCH2 */
      END

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	jmc	1.9	C $Header: /u/gcmpack/MITgcm/eesupp/src/exch_uv_xy_rx.template,v 1.8 2007/07/28 16:07:35 heimbach Exp $
2	cnh	1.3	C $Name: $
3	adcroft	1.2
4	jmc	1.5	#include "PACKAGES_CONFIG.h"
5	adcroft	1.2	#include "CPP_EEOPTIONS.h"
6	jmc	1.7	#ifdef ALLOW_EXCH2
7			#include "W2_OPTIONS.h"
8			#endif
9	adcroft	1.2
10	cnh	1.3	CBOP
11			C !ROUTINE: EXCH_UV_XY_RX
12
13			C !INTERFACE:
14	adcroft	1.2	SUBROUTINE EXCH_UV_XY_RX(
15	jmc	1.5	U Uphi, Vphi,
16			I withSigns, myThid )
17	adcroft	1.2	IMPLICIT NONE
18	cnh	1.3	C !DESCRIPTION:
19			C ==========================================================
20	jmc	1.7	C \| SUBROUTINE EXCH_UV_XY_RX
21			C \| o Handle exchanges for _RX, two-dimensional arrays.
22	cnh	1.3	C ==========================================================
23			C \| Driver exchange routine which branches to cube sphere or
24			C \| global, simple cartesian index grid. Exchange routine is
25			C \| called with two arrays that are components of a vector.
26	jmc	1.7	C \| These components are rotated and interchanged on the
27	cnh	1.3	C \| rotated grid during cube exchanges.
28			C ==========================================================
29	adcroft	1.2
30	cnh	1.3	C !USES:
31	adcroft	1.2	C === Global data ===
32			#include "SIZE.h"
33			#include "EEPARAMS.h"
34	heimbach	1.8	cph-mpi#include "EESUPPORT.h"
35	jmc	1.9	c#include "EXCH.h"
36	adcroft	1.2
37	cnh	1.3	C !INPUT/OUTPUT PARAMETERS:
38	adcroft	1.2	C === Routine arguments ===
39	cnh	1.3	C Uphi :: Arrays with overlap regions are to be exchanged
40			C Vphi Note - The interface to EXCH_ assumes that
41			C the standard Fortran 77 sequence association rules
42			C apply.
43			C myThid :: My thread id.
44			C withSigns :: Flag controlling whether vector is signed.
45	adcroft	1.2	_RX Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46			_RX Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47			LOGICAL withSigns
48			INTEGER myThid
49
50	cnh	1.3	C !LOCAL VARIABLES:
51	jmc	1.9	#ifndef ALLOW_EXCH2
52	adcroft	1.2	C == Local variables ==
53	cnh	1.3	C OL[wens] :: Overlap extents in west, east, north, south.
54			C exchWidth[XY] :: Extent of regions that will be exchanged.
55	adcroft	1.2	INTEGER OLw, OLe, OLn, OLs, exchWidthX, exchWidthY, myNz
56	jmc	1.9	#endif
57	cnh	1.3	CEOP
58	adcroft	1.2
59	jmc	1.5	#ifdef ALLOW_EXCH2
60	jmc	1.7	#ifdef W2_USE_R1_ONLY
61			CALL EXCH2_UV_CGRID_3D_RX(
62			U Uphi, Vphi,
63			I withSigns, 1, myThid )
64			#else
65	jmc	1.5	CALL EXCH2_UV_XY_RX(
66			U Uphi, Vphi,
67			I withSigns, myThid )
68	jmc	1.7	#endif
69	jmc	1.5	RETURN
70	jmc	1.7	#else /* ALLOW_EXCH2 */
71	jmc	1.5
72	adcroft	1.2	OLw = OLx
73			OLe = OLx
74			OLn = OLy
75			OLs = OLy
76			exchWidthX = OLx
77			exchWidthY = OLy
78			myNz = 1
79			C NOTE The exchange routine we use here does not
80			C require the preceeding and following barriers.
81			C However, the slow, simple exchange interface
82			C that is calling it here is meant to ensure
83			C that threads are synchronised before exchanges
84			C begine.
85			IF (useCubedSphereExchange) THEN
86			CALL EXCH_UV_RX_CUBE( Uphi, Vphi, withSigns,
87			I OLw, OLe, OLs, OLn, myNz,
88			I exchWidthX, exchWidthY,
89			I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
90			ELSE
91			c CALL EXCH_RX( Uphi,
92			c I OLw, OLe, OLs, OLn, myNz,
93			c I exchWidthX, exchWidthY,
94			c I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
95			c CALL EXCH_RX( Vphi,
96			c I OLw, OLe, OLs, OLn, myNz,
97			c I exchWidthX, exchWidthY,
98			c I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
99			c_jmc: for JAM compatibility, replace the 2 CALLs above by the 2 CPP_MACROs:
100			_EXCH_XY_RX( Uphi, myThid )
101			_EXCH_XY_RX( Vphi, myThid )
102			ENDIF
103
104			RETURN
105	jmc	1.7	#endif /* ALLOW_EXCH2 */
106	adcroft	1.2	END
107	jmc	1.5
108			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
109
110			CEH3 ;;; Local Variables: ***
111			CEH3 ;;; mode:fortran ***
112			CEH3 ;;; End: ***