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

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

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

revision 1.3 by edhill, Mon Apr 5 15:27:06 2004 UTC revision 1.4 by jmc, Tue Sep 21 21:10:45 2004 UTC
# Line 103  C          Uphi(snx+1,sny+1,bi,bj)= vPhi Line 103  C          Uphi(snx+1,sny+1,bi,bj)= vPhi
103           ENDIF           ENDIF
104    
105  C        Now zero out the null areas that should not be used in the numerics  C        Now zero out the null areas that should not be used in the numerics
106    C-       Also add one valid u,v value next to the corner, that allows
107    C         to compute vorticity on a wider stencil (e.g., vort3(0,1) & (1,0))
108           IF ( exch2_isWedge(myTile) .EQ. 1 .AND.           IF ( exch2_isWedge(myTile) .EQ. 1 .AND.
109       &        exch2_isSedge(myTile) .EQ. 1 ) THEN       &        exch2_isSedge(myTile) .EQ. 1 ) THEN
110  C         Zero SW corner points  C         Zero SW corner points
# Line 116  C         Zero SW corner points Line 118  C         Zero SW corner points
118              vPhi(I,J,bi,bj)=0.              vPhi(I,J,bi,bj)=0.
119             ENDDO             ENDDO
120            ENDDO            ENDDO
121                uPhi(0,0,bi,bj)=vPhi(1,0,bi,bj)
122                vPhi(0,0,bi,bj)=uPhi(0,1,bi,bj)
123           ENDIF           ENDIF
124           IF ( exch2_isWedge(myTile) .EQ. 1 .AND.           IF ( exch2_isWedge(myTile) .EQ. 1 .AND.
125       &        exch2_isNedge(myTile) .EQ. 1 ) THEN       &        exch2_isNedge(myTile) .EQ. 1 ) THEN
# Line 130  C         Zero NW corner points Line 134  C         Zero NW corner points
134              vPhi(I,J,bi,bj)=0.              vPhi(I,J,bi,bj)=0.
135             ENDDO             ENDDO
136            ENDDO            ENDDO
137              IF ( withSigns ) THEN
138                uPhi(0,sNy+1,bi,bj)=-vPhi(1,sNy+2,bi,bj)
139                vPhi(0,sNy+2,bi,bj)=-uPhi(0,sNy,bi,bj)
140              ELSE
141                uPhi(0,sNy+1,bi,bj)= vPhi(1,sNy+2,bi,bj)
142                vPhi(0,sNy+2,bi,bj)= uPhi(0,sNy,bi,bj)
143              ENDIF
144           ENDIF           ENDIF
145           IF ( exch2_isEedge(myTile) .EQ. 1 .AND.           IF ( exch2_isEedge(myTile) .EQ. 1 .AND.
146       &        exch2_isSedge(myTile) .EQ. 1 ) THEN       &        exch2_isSedge(myTile) .EQ. 1 ) THEN
# Line 144  C         Zero SE corner points Line 155  C         Zero SE corner points
155              vPhi(I,J,bi,bj)=0.              vPhi(I,J,bi,bj)=0.
156             ENDDO             ENDDO
157            ENDDO            ENDDO
158              IF ( withSigns ) THEN
159                uPhi(sNx+2,0,bi,bj)=-vPhi(sNx,0,bi,bj)
160                vPhi(sNx+1,0,bi,bj)=-uPhi(sNx+2,1,bi,bj)
161              ELSE
162                uPhi(sNx+2,0,bi,bj)= vPhi(sNx,0,bi,bj)
163                vPhi(sNx+1,0,bi,bj)= uPhi(sNx+2,1,bi,bj)
164              ENDIF
165           ENDIF           ENDIF
166           IF ( exch2_isEedge(myTile) .EQ. 1 .AND.           IF ( exch2_isEedge(myTile) .EQ. 1 .AND.
167       &        exch2_isNedge(myTile) .EQ. 1 ) THEN       &        exch2_isNedge(myTile) .EQ. 1 ) THEN
# Line 158  C         Zero NE corner points Line 176  C         Zero NE corner points
176              vPhi(I,J,bi,bj)=0.              vPhi(I,J,bi,bj)=0.
177             ENDDO             ENDDO
178            ENDDO            ENDDO
179                uPhi(sNx+2,sNy+1,bi,bj)=vPhi(sNx,sNy+2,bi,bj)
180                vPhi(sNx+1,sNy+2,bi,bj)=uPhi(sNx+2,sNy,bi,bj)
181           ENDIF           ENDIF
182          ENDDO          ENDDO
183         ENDDO         ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22