/[MITgcm]/MITgcm/model/inc/GRID.h
ViewVC logotype

Diff of /MITgcm/model/inc/GRID.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.22 by edhill, Thu Oct 9 04:19:18 2003 UTC revision 1.27 by jmc, Wed Apr 6 18:22:11 2005 UTC
# Line 296  C     *================================= Line 296  C     *=================================
296  C     \ev  C     \ev
297  CEOP  CEOP
298    
 #include "PACKAGES_CONFIG.h"  
   
299  C     Macros that override/modify standard definitions  C     Macros that override/modify standard definitions
300  #include "GRID_MACROS.h"  #include "GRID_MACROS.h"
301    
# Line 322  C              "lopped" a cell is (dimen Line 320  C              "lopped" a cell is (dimen
320  C              Note: The code needs terms like MIN(hFac,hFac(I+1))  C              Note: The code needs terms like MIN(hFac,hFac(I+1))
321  C                    On some platforms it may be better to precompute  C                    On some platforms it may be better to precompute
322  C                    hFacW, hFacE, ... here than do MIN on the fly.  C                    hFacW, hFacE, ... here than do MIN on the fly.
323  C     gravitySign - indicates whether gravity points in the opposite  C     gravitySign - indicates the direction of gravity relative to R direction
324  C                   direction of R or not.  C                  (= -1 for R=Z (Z increases upward, -gravity direction  )
325  C                 ( = +1 for R=Z (gravity points downward in Z)  C                  (= +1 for R=P (P increases downward, +gravity direction)
 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 )
 C                 ( vertical coord == m  -> rkFac =  1 )  
 C                 ( vertical coord == Pa -> rkFac = -1 )  
328  C     maskH  - cell Center full-column mask (= 2D mask)  C     maskH  - cell Center full-column mask (= 2D mask)
329  C     maskC  - cell Center land mask  C     maskC  - cell Center land mask
330  C     maskW  - West face land mask  C     maskW  - West face land mask
# Line 369  C     tanPhiAtU - tan of the latitude at Line 364  C     tanPhiAtU - tan of the latitude at
364  C                 metric term in U equation.  C                 metric term in U equation.
365  C     tanPhiAtV - tan of the latitude at V point. Used for spherical polar  C     tanPhiAtV - tan of the latitude at V point. Used for spherical polar
366  C                 metric term in V equation.  C                 metric term in V equation.
367    C     fCori     :: Coriolis parameter at grid Center point
368    C     fCoriG    :: Coriolis parameter at grid Corner point
369    C     fCoriCos  :: Coriolis Cos(phi) parameter at grid Center point (for NH)
370        COMMON /GRID_I/ klowC        COMMON /GRID_I/ klowC
371        INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
372    
373        COMMON /GRID_R/        COMMON /GRID_R/
374         &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,
375         &  gravitySign,rkFac, recip_rkFac,
376       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,
377       &  R_low,Ro_surf,HFacC,HFacW,HFacS,       &  R_low,Ro_surf,HFacC,HFacW,HFacS,
378       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,
# Line 384  C                 metric term in V equat Line 384  C                 metric term in V equat
384       &  maskH, maskC,maskW,maskS,       &  maskH, maskC,maskW,maskS,
385       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,
386       &  tanPhiAtU, tanPhiAtV,       &  tanPhiAtU, tanPhiAtV,
      &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,  
387       &  drC,drF,recip_drC,recip_drF,rC,rF,       &  drC,drF,recip_drC,recip_drF,rC,rF,
388       &  gravitySign,       &  xC0, yC0,
389       &  rkFac, recip_rkFac, xC0, yC0       &  fCori, fCoriG, fCoriCos
390          _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)
391          _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)
392          _RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy)
393          _RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy)
394          _RL gravitySign
395          _RL rkFac
396          _RL recip_rkFac
397        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
398        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
399        _RS dxG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 431  C                 metric term in V equat Line 437  C                 metric term in V equat
437        _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)
438        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
439        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)  
       _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)  
       _RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy)  
       _RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy)  
440        _RS drC            (1:Nr)        _RS drC            (1:Nr)
441        _RS drF            (1:Nr)        _RS drF            (1:Nr)
442        _RS recip_drC      (1:Nr)        _RS recip_drC      (1:Nr)
# Line 442  C                 metric term in V equat Line 444  C                 metric term in V equat
444        _RS saFac          (1:Nr)        _RS saFac          (1:Nr)
445        _RS rC             (1:Nr)        _RS rC             (1:Nr)
446        _RS rF             (1:Nr+1)        _RS rF             (1:Nr+1)
       _RL gravitySign  
       _RS rkFac  
       _RS recip_rkFac  
447        _RS xC0        _RS xC0
448        _RS yC0        _RS yC0
449          _RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
450          _RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
451          _RS fCoriCos(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
452    
453    
454  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
455        COMMON /GRID_NH/        COMMON /GRID_NH/
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22