104 |
C | "WB" |
C | "WB" |
105 |
C | + |
C | + |
106 |
C | + |
C | + |
107 |
C | \+/ /|\ /|\ |
C | \+/ /|\ /|\ . |
108 |
C |"UB"++>|-------- | -----------------etc.. | ----*--- |
C |"UB"++>|-------- | -----------------etc.. | ----*--- |
109 |
C | | rVel(i=1, | rVel(i=Nx, * | |
C | | rVel(i=1, | rVel(i=Nx, * | |
110 |
C | | j=1, | j=1, * | |
C | | j=1, | j=1, * | |
145 |
C | "PWX" |
C | "PWX" |
146 |
C | |
C | |
147 |
C | Up increasing upwards. |
C | Up increasing upwards. |
148 |
C |/|\ |
C |/|\ . |
149 |
C | | |
C | | |
150 |
C | | |
C | | |
151 |
C | =====> E i increasing eastwards |
C | =====> E i increasing eastwards |
231 |
C | \|/ | |<---dxV(i=2,j=1,k=1)--\|/ | |
C | \|/ | |<---dxV(i=2,j=1,k=1)--\|/ | |
232 |
C |"SB"++>|___________v___________|___________v__________|_ |
C |"SB"++>|___________v___________|___________v__________|_ |
233 |
C | <--dxG(i=1,j=1,k=1)-----> |
C | <--dxG(i=1,j=1,k=1)-----> |
234 |
C | /+\ |
C | /+\ . |
235 |
C | + |
C | + |
236 |
C | + |
C | + |
237 |
C | "WB" |
C | "WB" |
287 |
C | | | | |
C | | | | |
288 |
C | | \|/ | |
C | | \|/ | |
289 |
C |"SB"++>|___________v___________|___________v__etc... |
C |"SB"++>|___________v___________|___________v__etc... |
290 |
C | /+\ |
C | /+\ . |
291 |
C | + |
C | + |
292 |
C | + |
C | + |
293 |
C | "WB" |
C | "WB" |
319 |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
320 |
C On some platforms it may be better to precompute |
C On some platforms it may be better to precompute |
321 |
C hFacW, hFacE, ... here than do MIN on the fly. |
C hFacW, hFacE, ... here than do MIN on the fly. |
322 |
|
C gravitySign - indicates whether gravity points in the opposite |
323 |
|
C direction of R or not. |
324 |
|
C ( = +1 for R=Z (gravity points downward in Z) |
325 |
|
C ( = -1 for R=P (gravity points upward in P) |
326 |
C rkFac - Vertical coordinate to vertical index orientation. |
C rkFac - Vertical coordinate to vertical index orientation. |
327 |
C ( -1 same orientation, 1 opposite orientation ) |
C ( -1 same orientation, 1 opposite orientation ) |
328 |
C ( vertical coord == m -> rkFac = 1 ) |
C ( vertical coord == m -> rkFac = 1 ) |
329 |
C ( vertical coord == Pa -> rkFac = -1 ) |
C ( vertical coord == Pa -> rkFac = -1 ) |
330 |
|
C maskH - cell Center full-column mask (= 2D mask) |
331 |
C maskC - cell Center land mask |
C maskC - cell Center land mask |
332 |
C maskW - West face land mask |
C maskW - West face land mask |
333 |
C maskS - South face land mask |
C maskS - South face land mask |
375 |
& recip_hFacC,recip_hFacW,recip_hFacS, |
& recip_hFacC,recip_hFacW,recip_hFacS, |
376 |
& saFac, |
& saFac, |
377 |
& xC,yC,rA,rAw,rAs,rAz,xG,yG, |
& xC,yC,rA,rAw,rAs,rAz,xG,yG, |
378 |
& maskC,maskW,maskS,recip_rA,recip_rAw,recip_rAs,recip_rAz, |
& maskH, maskC,maskW,maskS, |
379 |
|
& recip_rA,recip_rAw,recip_rAs,recip_rAz, |
380 |
& tanPhiAtU, tanPhiAtV, |
& tanPhiAtU, tanPhiAtV, |
381 |
& cosfacU,cosfacV,sqcosfacU,sqcosfacV, |
& cosfacU,cosfacV,sqcosfacU,sqcosfacV, |
382 |
& drC,drF,recip_drC,recip_drF,rC,rF, |
& drC,drF,recip_drC,recip_drF,rC,rF, |
383 |
|
& gravitySign, |
384 |
& rkFac, recip_rkFac, xC0, yC0 |
& rkFac, recip_rkFac, xC0, yC0 |
385 |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
386 |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
419 |
_RS recip_rAw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS recip_rAw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
420 |
_RS recip_rAs (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS recip_rAs (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
421 |
_RS recip_rAz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS recip_rAz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
422 |
|
_RS maskH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
423 |
_RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
_RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
424 |
_RS maskW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
_RS maskW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
425 |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
436 |
_RS saFac (1:Nr) |
_RS saFac (1:Nr) |
437 |
_RS rC (1:Nr) |
_RS rC (1:Nr) |
438 |
_RS rF (1:Nr+1) |
_RS rF (1:Nr+1) |
439 |
|
_RL gravitySign |
440 |
_RS rkFac |
_RS rkFac |
441 |
_RS recip_rkFac |
_RS recip_rkFac |
442 |
_RS xC0 |
_RS xC0 |