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

Diff of /MITgcm/pkg/exch2/exch2_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.6 by jmc, Thu May 6 23:28:45 2010 UTC
# Line 56  CEOP Line 56  CEOP
56        exchWidthX = OLx        exchWidthX = OLx
57        exchWidthY = OLy        exchWidthY = OLy
58    
59        CALL EXCH2_RX1_CUBE( phi, 'T ',        CALL EXCH2_RX1_CUBE( phi, .FALSE., 'T ',
60       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
61       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
62       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_UPDATE_CORNERS, myThid )
63    
64  #ifdef W2_FILL_NULL_REGIONS        CALL EXCH2_RX1_CUBE( phi, .FALSE., 'T ',
       IF (useCubedSphereExchange) THEN  
        DO bj=myByLo(myThid),myByHi(myThid)  
         DO bi=myBxLo(myThid),myBxHi(myThid)  
          myTile = W2_myTileList(bi,bj)  
 C        South-east corner  
          IF ( exch2_isEedge(myTile) .EQ. 1 .AND.  
      &        exch2_isSedge(myTile) .EQ. 1 ) THEN  
           DO j=1-OLy,0  
            DO i=sNx+1,sNx+OLx  
             DO k=1,myNz  
              phi(i,j,k,bi,bj)=e2FillValue_RX  
             ENDDO  
            ENDDO  
           ENDDO  
          ENDIF  
 C        North-east corner  
          IF ( exch2_isEedge(myTile) .EQ. 1 .AND.  
      &        exch2_isNedge(myTile) .EQ. 1 ) THEN  
           DO j=sNy+1,sNy+OLy  
            DO i=sNx+1,sNx+OLx  
             DO k=1,myNz  
              phi(i,j,k,bi,bj)=e2FillValue_RX  
             ENDDO  
            ENDDO  
           ENDDO  
          ENDIF  
 C        South-west corner  
          IF ( exch2_isWedge(myTile) .EQ. 1 .AND.  
      &        exch2_isSedge(myTile) .EQ. 1 ) THEN  
           DO j=1-OLy,0  
            DO i=1-OLx,0  
             DO k=1,myNz  
              phi(i,j,k,bi,bj)=e2FillValue_RX  
             ENDDO  
            ENDDO  
           ENDDO  
          ENDIF  
 C        North-west corner  
          IF ( exch2_isWedge(myTile) .EQ. 1 .AND.  
      &        exch2_isNedge(myTile) .EQ. 1 ) THEN  
           DO j=sNy+1,sNy+OLy  
            DO i=1-OLx,0  
             DO k=1,myNz  
              phi(i,j,k,bi,bj)=e2FillValue_RX  
             ENDDO  
            ENDDO  
           ENDDO  
          ENDIF  
         ENDDO  
        ENDDO  
       ENDIF  
 #endif /* W2_FILL_NULL_REGIONS */  
   
       CALL EXCH2_RX1_CUBE( phi, 'T ',  
65       I            OLw, OLe, OLs, OLn, myNz,       I            OLw, OLe, OLs, OLn, myNz,
66       I            exchWidthX, exchWidthY,       I            exchWidthX, exchWidthY,
67       I            FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )       I            EXCH_UPDATE_CORNERS, myThid )
68    
69  #ifdef W2_FILL_NULL_REGIONS  #ifdef W2_FILL_NULL_REGIONS
70        IF (useCubedSphereExchange) THEN        IF (useCubedSphereExchange) THEN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22