81 |
C Note: In the spherical polar case delX and delY are given in |
C Note: In the spherical polar case delX and delY are given in |
82 |
C degrees and are relative to some starting latitude and |
C degrees and are relative to some starting latitude and |
83 |
C longitude - phiMin and thetaMin. |
C longitude - phiMin and thetaMin. |
84 |
|
xC0 = thetaMin |
85 |
|
yC0 = phiMin |
86 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
87 |
jG = myYGlobalLo + (bj-1)*sNy |
jG = myYGlobalLo + (bj-1)*sNy |
88 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
199 |
ENDDO |
ENDDO |
200 |
ENDDO |
ENDDO |
201 |
ENDDO |
ENDDO |
202 |
_EXCH_XY_R4(zA, myThid ) |
|
203 |
|
DO bj = myByLo(myThid), myByHi(myThid) |
204 |
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
205 |
|
DO K=1,Nz |
206 |
|
DO J=1,sNy |
207 |
|
DO I=1,sNx |
208 |
|
IF (HFacC(I,J,K,bi,bj) .NE. 0. D0 ) THEN |
209 |
|
rHFacC(I,J,K,bi,bj) = 1. D0 / HFacC(I,J,K,bi,bj) |
210 |
|
ELSE |
211 |
|
rHFacC(I,J,K,bi,bj) = 0. D0 |
212 |
|
ENDIF |
213 |
|
IF (HFacW(I,J,K,bi,bj) .NE. 0. D0 ) THEN |
214 |
|
rHFacW(I,J,K,bi,bj) = 1. D0 / HFacW(I,J,K,bi,bj) |
215 |
|
maskW(I,J,K,bi,bj) = 1. D0 |
216 |
|
ELSE |
217 |
|
rHFacW(I,J,K,bi,bj) = 0. D0 |
218 |
|
maskW(I,J,K,bi,bj) = 0.0 D0 |
219 |
|
ENDIF |
220 |
|
IF (HFacS(I,J,K,bi,bj) .NE. 0. D0 ) THEN |
221 |
|
rHFacS(I,J,K,bi,bj) = 1. D0 / HFacS(I,J,K,bi,bj) |
222 |
|
maskS(I,J,K,bi,bj) = 1. D0 |
223 |
|
ELSE |
224 |
|
rHFacS(I,J,K,bi,bj) = 0. D0 |
225 |
|
maskS(I,J,K,bi,bj) = 0. D0 |
226 |
|
ENDIF |
227 |
|
ENDDO |
228 |
|
ENDDO |
229 |
|
ENDDO |
230 |
|
ENDDO |
231 |
|
ENDDO |
232 |
|
C Now sync. and get/send edge regions that are shared with |
233 |
|
C other threads. |
234 |
|
_EXCH_XYZ_R4(rHFacC , myThid ) |
235 |
|
_EXCH_XYZ_R4(rHFacW , myThid ) |
236 |
|
_EXCH_XYZ_R4(rHFacS , myThid ) |
237 |
|
_EXCH_XYZ_R4(maskW , myThid ) |
238 |
|
_EXCH_XYZ_R4(maskS , myThid ) |
239 |
|
_EXCH_XY_R4 (zA , myThid ) |
240 |
|
|
241 |
|
CcnhDebugStarts |
242 |
|
tanPhiAtU = 0. _d 0 |
243 |
|
tanPhiAtV = 0. _d 0 |
244 |
|
_EXCH_XY_R4 (tanPhiAtU , myThid ) |
245 |
|
_EXCH_XY_R4 (tanPhiAtV , myThid ) |
246 |
|
CcnhDebugEnds |
247 |
|
|
248 |
|
|
249 |
C |
C |
250 |
RETURN |
RETURN |
251 |
END |
END |