/[MITgcm]/MITgcm/pkg/exch2/exch2_uv_bgrid_3d_rx.template
ViewVC logotype

Diff of /MITgcm/pkg/exch2/exch2_uv_bgrid_3d_rx.template

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

revision 1.5 by jmc, Sun Jun 28 00:57:51 2009 UTC revision 1.8 by jmc, Mon Mar 26 19:43:10 2012 UTC
# Line 58  C     exchWidth[XY] :: Extent of regions Line 58  C     exchWidth[XY] :: Extent of regions
58  C     uPhiNW,uPhiSE :: temporary array to hold corner value (CS grid)  C     uPhiNW,uPhiSE :: temporary array to hold corner value (CS grid)
59  C     vPhiNW,vPhiSE :: temporary array to hold corner value (CS grid)  C     vPhiNW,vPhiSE :: temporary array to hold corner value (CS grid)
60  C     uLoc,vLoc     :: local copy of the vector components with haloes filled.  C     uLoc,vLoc     :: local copy of the vector components with haloes filled.
61  C     msgBuf        :: Informational/error meesage buffer  C     msgBuf        :: Informational/error message buffer
62    
63        INTEGER local_maxDim        INTEGER local_maxDim
64        PARAMETER( local_maxDim = 8*Nr )        PARAMETER( local_maxDim = 8*Nr )
# Line 112  C---  using or not using CubedSphereExch Line 112  C---  using or not using CubedSphereExch
112    
113  C--   First call the exchanges for the two components  C--   First call the exchanges for the two components
114    
115         CALL EXCH2_RX1_CUBE( uPhi, 'T ',         CALL EXCH2_RX1_CUBE( uPhi, .FALSE., 'T ',
116       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
117       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
118       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_IGNORE_CORNERS, myThid )
119         CALL EXCH2_RX1_CUBE( uPhi, 'T ',         CALL EXCH2_RX1_CUBE( uPhi, .FALSE., 'T ',
120       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
121       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
122       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_UPDATE_CORNERS, myThid )
123    
124         CALL EXCH2_RX1_CUBE( vPhi, 'T ',         CALL EXCH2_RX1_CUBE( vPhi, .FALSE., 'T ',
125       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
126       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
127       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_IGNORE_CORNERS, myThid )
128         CALL EXCH2_RX1_CUBE( vPhi, 'T ',         CALL EXCH2_RX1_CUBE( vPhi, .FALSE., 'T ',
129       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
130       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
131       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_UPDATE_CORNERS, myThid )
132    
133  C- note: can substitute the low-level S/R calls above with:  C- note: can substitute the low-level S/R calls above with:
134  c      CALL EXCH2_3D_RX( uPhi, myNz, myThid )  c      CALL EXCH2_3D_RX( uPhi, myNz, myThid )
# Line 149  C--   Loops on tile indices: Line 149  C--   Loops on tile indices:
149         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
150          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
151    
152  C-    Now choose what to do at each edge of the halo based on which face  C-    Choose what to do at each edge of the halo based on which face we are
 C    (we assume that bj is always=1)  
153           myTile = W2_myTileList(bi,bj)           myTile = W2_myTileList(bi,bj)
154           myFace = exch2_myFace(myTile)           myFace = exch2_myFace(myTile)
155    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22