/[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.4 by jmc, Tue May 12 19:44:58 2009 UTC revision 1.5 by jmc, Sun Jun 28 00:57:51 2009 UTC
# Line 35  C     !USES: Line 35  C     !USES:
35    
36  #include "SIZE.h"  #include "SIZE.h"
37  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c#include "EESUPPORT.h"  
38  #include "W2_EXCH2_SIZE.h"  #include "W2_EXCH2_SIZE.h"
39  #include "W2_EXCH2_TOPOLOGY.h"  #include "W2_EXCH2_TOPOLOGY.h"
40  #ifdef W2_FILL_NULL_REGIONS  #ifdef W2_FILL_NULL_REGIONS
# Line 108  C-     save 2 corners value (in case we Line 107  C-     save 2 corners value (in case we
107           ENDDO           ENDDO
108          ENDDO          ENDDO
109         ENDDO         ENDDO
110    C---  using or not using CubedSphereExchange: end
111          ENDIF
112    
113  C--   First call the exchanges for the two components  C--   First call the exchanges for the two components
114    
# Line 133  C- note: can substitute the low-level S/ Line 134  C- note: can substitute the low-level S/
134  c      CALL EXCH2_3D_RX( uPhi, myNz, myThid )  c      CALL EXCH2_3D_RX( uPhi, myNz, myThid )
135  c      CALL EXCH2_3D_RX( vPhi, myNz, myThid )  c      CALL EXCH2_3D_RX( vPhi, myNz, myThid )
136    
137          IF ( useCubedSphereExchange ) THEN
138    
139  C--   Then, depending on which tile we are, we may need  C--   Then, depending on which tile we are, we may need
140  C     1) to switch u and v components and also to switch the signs  C     1) to switch u and v components and also to switch the signs
141  C     2) to shift the index along the face edge.  C     2) to shift the index along the face edge.
# Line 148  C--   Loops on tile indices: Line 151  C--   Loops on tile indices:
151    
152  C-    Now choose what to do at each edge of the halo based on which face  C-    Now choose what to do at each edge of the halo based on which face
153  C    (we assume that bj is always=1)  C    (we assume that bj is always=1)
154           myTile = W2_myTileList(bi)           myTile = W2_myTileList(bi,bj)
155           myFace = exch2_myFace(myTile)           myFace = exch2_myFace(myTile)
156    
157  C--   Loops on level index:  C--   Loops on level index:
# Line 407  C--   end of Loops on tile indices (bi,b Line 410  C--   end of Loops on tile indices (bi,b
410          ENDDO          ENDDO
411         ENDDO         ENDDO
412    
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  
   
       ELSE  
 C---  not using CubedSphereExchange:  
   
 #ifndef AUTODIFF_EXCH2  
        CALL EXCH_RX( uPhi,  
      I            OLw, OLe, OLs, OLn, myNz,  
      I            exchWidthX, exchWidthY,  
      I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )  
        CALL EXCH_RX( vPhi,  
      I            OLw, OLe, OLs, OLn, myNz,  
      I            exchWidthX, exchWidthY,  
      I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )  
 #endif  
   
413  C---  using or not using CubedSphereExchange: end  C---  using or not using CubedSphereExchange: end
414        ENDIF        ENDIF
415    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22