349 |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
350 |
C On some platforms it may be better to precompute |
C On some platforms it may be better to precompute |
351 |
C hFacW, hFacE, ... here than do MIN on the fly. |
C hFacW, hFacE, ... here than do MIN on the fly. |
352 |
C maskH :: cell Center full-column mask (= 2D mask) |
cC maskH :: cell Center full-column mask (= 2D mask) |
353 |
C maskInC :: Cell Center 2-D Interior mask (i.e., zero beyond OB) |
C maskInC :: Cell Center 2-D Interior mask (i.e., zero beyond OB) |
354 |
C maskInW :: West face 2-D Interior mask (i.e., zero on and beyond OB) |
C maskInW :: West face 2-D Interior mask (i.e., zero on and beyond OB) |
355 |
C maskInS :: South face 2-D Interior mask (i.e., zero on and beyond OB) |
C maskInS :: South face 2-D Interior mask (i.e., zero on and beyond OB) |
405 |
& recip_Rcol, |
& recip_Rcol, |
406 |
& recip_hFacC,recip_hFacW,recip_hFacS, |
& recip_hFacC,recip_hFacW,recip_hFacS, |
407 |
& xC,yC,rA,rAw,rAs,rAz,xG,yG, |
& xC,yC,rA,rAw,rAs,rAz,xG,yG, |
|
& maskH, |
|
408 |
& maskInC, maskInW, maskInS, |
& maskInC, maskInW, maskInS, |
409 |
& maskC, maskW, maskS, |
& maskC, maskW, maskS, |
410 |
& recip_rA,recip_rAw,recip_rAs,recip_rAz, |
& recip_rA,recip_rAw,recip_rAs,recip_rAz, |
411 |
& tanPhiAtU, tanPhiAtV, angleCosC, angleSinC, |
& tanPhiAtU, tanPhiAtV, angleCosC, angleSinC, |
412 |
& drC,drF,recip_drC,recip_drF,rC,rF, |
& drC,drF,recip_drC,recip_drF,rC,rF, |
413 |
& fCori, fCoriG, fCoriCos |
& fCori, fCoriG, fCoriCos |
414 |
|
c & , maskH |
415 |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
416 |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
417 |
_RS dxG (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxG (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
449 |
_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) |
450 |
_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) |
451 |
_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) |
452 |
_RS maskH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
c _RS maskH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
453 |
_RS maskInC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS maskInC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
454 |
_RS maskInW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS maskInW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
455 |
_RS maskInS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS maskInS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |