315 |
& deepFacC, deepFac2C, recip_deepFacC, recip_deepFac2C, |
& deepFacC, deepFac2C, recip_deepFacC, recip_deepFac2C, |
316 |
& deepFacF, deepFac2F, recip_deepFacF, recip_deepFac2F, |
& deepFacF, deepFac2F, recip_deepFacF, recip_deepFac2F, |
317 |
& gravitySign, rkSign, globalArea |
& gravitySign, rkSign, globalArea |
318 |
_RL cosFacU (1-Oly:sNy+Oly,nSx,nSy) |
_RL cosFacU (1-OLy:sNy+OLy,nSx,nSy) |
319 |
_RL cosFacV (1-Oly:sNy+Oly,nSx,nSy) |
_RL cosFacV (1-OLy:sNy+OLy,nSx,nSy) |
320 |
_RL sqCosFacU (1-Oly:sNy+Oly,nSx,nSy) |
_RL sqCosFacU (1-OLy:sNy+OLy,nSx,nSy) |
321 |
_RL sqCosFacV (1-Oly:sNy+Oly,nSx,nSy) |
_RL sqCosFacV (1-OLy:sNy+OLy,nSx,nSy) |
322 |
_RL deepFacC (Nr) |
_RL deepFacC (Nr) |
323 |
_RL deepFac2C (Nr) |
_RL deepFac2C (Nr) |
324 |
_RL deepFacF (Nr+1) |
_RL deepFacF (Nr+1) |
350 |
C rSurfS :: surface reference position at Southern edge location [r_unit]. |
C rSurfS :: surface reference position at Southern edge location [r_unit]. |
351 |
C hFac :: Fraction of cell in vertical which is open i.e how |
C hFac :: Fraction of cell in vertical which is open i.e how |
352 |
C "lopped" a cell is (dimensionless scale factor). |
C "lopped" a cell is (dimensionless scale factor). |
353 |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
C Note: The code needs terms like MIN(hFac,hFac(I-1)) |
354 |
C On some platforms it may be better to precompute |
C On some platforms it may be better to precompute |
355 |
C hFacW, hFacE, ... here than do MIN on the fly. |
C hFacW, hFacS, ... here than do MIN on the fly. |
356 |
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) |
357 |
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) |
358 |
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) |
381 |
C yC :: Y-coordinate of center of cell f[X,Y]. |
C yC :: Y-coordinate of center of cell f[X,Y]. |
382 |
C yG :: Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y]. |
C yG :: Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y]. |
383 |
C rA :: R-face are f[X,Y] ( m^2 ). |
C rA :: R-face are f[X,Y] ( m^2 ). |
384 |
C Note: In a cartesian framework zA is simply dx*dy, |
C Note: In a cartesian framework rA is simply dx*dy, |
385 |
C however we use zA to allow for non-globally |
C however we use rA to allow for non-globally |
386 |
C orthogonal coordinate frames (with appropriate |
C orthogonal coordinate frames (with appropriate |
387 |
C metric terms). |
C metric terms). |
388 |
C rC :: R-coordinate of center of cell f[Z] (units of r). |
C rC :: R-coordinate of center of cell f[Z] (units of r). |
476 |
_RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
477 |
_RS maskW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RS maskW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
478 |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
479 |
_RS drC (Nr) |
_RS drC (Nr+1) |
480 |
_RS drF (Nr) |
_RS drF (Nr) |
481 |
_RS recip_drC (Nr) |
_RS recip_drC (Nr+1) |
482 |
_RS recip_drF (Nr) |
_RS recip_drF (Nr) |
483 |
_RS rC (Nr) |
_RS rC (Nr) |
484 |
_RS rF (Nr+1) |
_RS rF (Nr+1) |
513 |
C kSurfW :: vertical index of the surface U point |
C kSurfW :: vertical index of the surface U point |
514 |
C kSurfS :: vertical index of the surface V point |
C kSurfS :: vertical index of the surface V point |
515 |
C kLowC :: index of the r-lowest "wet cell" (2D) |
C kLowC :: index of the r-lowest "wet cell" (2D) |
516 |
C IMPORTANT: ksurfC,W,S = Nr+1 and kLowC = 0 where the fluid column |
C IMPORTANT: kSurfC,W,S = Nr+1 and kLowC = 0 where the fluid column |
517 |
C is empty (continent) |
C is empty (continent) |
518 |
COMMON /GRID_I/ |
COMMON /GRID_I/ |
519 |
& kSurfC, kSurfW, kSurfS, |
& kSurfC, kSurfW, kSurfS, |