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 |
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 |
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 |
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 |
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 |