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 |