/[MITgcm]/MITgcm/eesupp/src/exch_uv_xy_rx.template

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

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

-revision 1.7 by jmc,
Wed Jul 25 21:19:44 2007 UTC
+revision 1.12 by jmc,
Mon May 14 13:15:05 2012 UTC
 Line 12 
 C     !ROUTINE: EXCH_UV_XY_RX
  C     !INTERFACE:
        SUBROUTINE EXCH_UV_XY_RX(
-      U                       Uphi, Vphi,
+      U                          uPhi, vPhi,
-      I                       withSigns, myThid )
+      I                          withSigns, myThid )
-       IMPLICIT NONE
  C     !DESCRIPTION:
  C     *==========================================================*
  C     | SUBROUTINE EXCH_UV_XY_RX
- C     | o Handle exchanges for _RX, two-dimensional arrays.
+ C     | o Handle exchanges for _RX, 2-dimensional vector arrays.
  C     *==========================================================*
- C     | Driver exchange routine which branches to cube sphere or
+ C     | Vector arrays need to be rotated and interchaged for
- C     | global, simple cartesian index grid. Exchange routine is
+ C     | exchange operations on some grids. This driver routine
- C     | called with two arrays that are components of a vector.
+ C     | branches to support this.
- C     | These components are rotated and interchanged on the
- C     | rotated grid during cube exchanges.
  C     *==========================================================*
  C     !USES:
+       IMPLICIT NONE
  C     === Global data ===
  #include "SIZE.h"
  #include "EEPARAMS.h"
- #include "EESUPPORT.h"
- #include "EXCH.h"
  C     !INPUT/OUTPUT PARAMETERS:
  C     === Routine arguments ===
- C     Uphi      :: Arrays with overlap regions are to be exchanged
+ C     uPhi      :: 2 components of a vector field with overlap regions
- C     Vphi         Note - The interface to EXCH_ assumes that
+ C     vPhi      :: to be exchanged
- 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)
+ C     myThid    :: my Thread Id. number
-       _RX Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+       _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,
+      U                     uPhi, vPhi,
-      I                       withSigns, 1, myThid )
+      I                     withSigns, 1, myThid )
  #else
-       CALL EXCH2_UV_XY_RX(
+       CALL EXCH2_UV_3D_RX(
-      U                       Uphi, Vphi,
+      U                     uPhi, vPhi,
-      I                       withSigns, myThid )
+      I                     withSigns, 1, myThid )
  #endif
-       RETURN
  #else /* ALLOW_EXCH2 */
        OLw        = OLx
-Line 74 
 CEOP
+Line 71 
 CEOP
        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,
+        CALL EXCH1_UV_RX_CUBE( uPhi, vPhi, withSigns,
       I            OLw, OLe, OLs, OLn, myNz,
       I            exchWidthX, exchWidthY,
-      I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
+      I            EXCH_UPDATE_CORNERS, myThid )
        ELSE
- c      CALL EXCH_RX( Uphi,
+ #ifdef DISCONNECTED_TILES
- c    I            OLw, OLe, OLs, OLn, myNz,
+        CALL EXCH0_RX( uPhi,
- c    I            exchWidthX, exchWidthY,
+      I            OLw, OLe, OLs, OLn, myNz,
- c    I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
+      I            exchWidthX, exchWidthY,
- c      CALL EXCH_RX( Vphi,
+      I            EXCH_UPDATE_CORNERS, myThid )
- c    I            OLw, OLe, OLs, OLn, myNz,
+        CALL EXCH0_RX( vPhi,
- c    I            exchWidthX, exchWidthY,
+      I            OLw, OLe, OLs, OLn, myNz,
- c    I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
+      I            exchWidthX, exchWidthY,
- c_jmc: for JAM compatibility, replace the 2 CALLs above by the 2 CPP_MACROs:
+      I            EXCH_UPDATE_CORNERS, myThid )
-        _EXCH_XY_RX( Uphi, myThid )
+ #else /* DISCONNECTED_TILES */
-        _EXCH_XY_RX( Vphi, myThid )
+        CALL EXCH1_RX( uPhi,
+      I            OLw, OLe, OLs, OLn, myNz,
+      I            exchWidthX, exchWidthY,
+      I            EXCH_UPDATE_CORNERS, myThid )
+        CALL EXCH1_RX( vPhi,
+      I            OLw, OLe, OLs, OLn, myNz,
+      I            exchWidthX, exchWidthY,
+      I            EXCH_UPDATE_CORNERS, myThid )
+ #endif /* DISCONNECTED_TILES */
        ENDIF
-       RETURN
  #endif /* ALLOW_EXCH2 */
+       RETURN
        END
  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

 Legend:



Removed from v.1.7
 


changed lines


 
Added in v.1.12
 Legend:



Removed from v.1.7
 


changed lines


 
Added in v.1.12
-Removed from v.1.7
+Added in v.1.12

	ViewVC Help
Powered by ViewVC 1.1.22