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

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

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

revision 1.3 by jmc, Tue May 12 19:44:58 2009 UTC revision 1.4 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    
# Line 67  CEOP Line 66  CEOP
66        negOne = 1.        negOne = 1.
67        IF (withSigns) negOne = -1.        IF (withSigns) negOne = -1.
68    
       IF ( useCubedSphereExchange ) THEN  
 C---  using CubedSphereExchange:  
   
69  C--   First call the exchanges  C--   First call the exchanges
70    
71         CALL EXCH2_RX1_CUBE( phi, 'T ',         CALL EXCH2_RX1_CUBE( phi, 'T ',
# Line 84  C--   First call the exchanges Line 80  C--   First call the exchanges
80  C- note: can substitute the low-level S/R call above with:  C- note: can substitute the low-level S/R call above with:
81  c      CALL EXCH2_3D_RX( phi, myNz, myThid )  c      CALL EXCH2_3D_RX( phi, myNz, myThid )
82    
83          IF ( useCubedSphereExchange ) THEN
84    C---  using CubedSphereExchange:
85    
86  C--   Then we may need to switch u and v components  C--   Then we may need to switch u and v components
87  C     and/or the signs depending on which cube face we are located.  C     and/or the signs depending on which cube face we are located.
88    
# Line 94  C--   Loops on tile indices: Line 93  C--   Loops on tile indices:
93  C--   Then we may need to switch the signs depending on which cube face  C--   Then we may need to switch the signs depending on which cube face
94  C      we are located (we assume that bj is always=1).  C      we are located (we assume that bj is always=1).
95  C     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
96           myTile = W2_myTileList(bi)           myTile = W2_myTileList(bi,bj)
97           myFace = exch2_myFace(myTile)           myFace = exch2_myFace(myTile)
98    
99  C--   Loops on level index:  C--   Loops on level index:
# Line 181  C--    end of Loops on tile and level in Line 180  C--    end of Loops on tile and level in
180          ENDDO          ENDDO
181         ENDDO         ENDDO
182    
       ELSE  
 C---  not using CubedSphereExchange:  
   
 #ifndef AUTODIFF_EXCH2  
        CALL EXCH_RX( phi,  
      I            OLw, OLe, OLs, OLn, myNz,  
      I            exchWidthX, exchWidthY,  
      I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )  
 #endif  
   
183  C---  using or not using CubedSphereExchange: end  C---  using or not using CubedSphereExchange: end
184        ENDIF        ENDIF
185    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22