/[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.17 by adcroft, Thu Sep 27 18:18:31 2001 UTC revision 1.22 by edhill, Thu Oct 9 04:19:18 2003 UTC
# Line 296  C     *================================= Line 296  C     *=================================
296  C     \ev  C     \ev
297  CEOP  CEOP
298    
299    #include "PACKAGES_CONFIG.h"
300    
301  C     Macros that override/modify standard definitions  C     Macros that override/modify standard definitions
302  #include "GRID_MACROS.h"  #include "GRID_MACROS.h"
303    
# Line 314  C     drF    - Cell face separation alon Line 316  C     drF    - Cell face separation alon
316  C     Rcolumn  -Total thickness (in r_unit) of the fluid column  C     Rcolumn  -Total thickness (in r_unit) of the fluid column
317  C     R_low  - base of fluid in r_unit (Depth(m) / Pressure(Pa) at top Atmos.)  C     R_low  - base of fluid in r_unit (Depth(m) / Pressure(Pa) at top Atmos.)
318  C     Ro_surf- surface reference (at rest) position, r_unit.  C     Ro_surf- surface reference (at rest) position, r_unit.
319    C     klowC  - index of the lowest "wet cell" (2D)
320  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
321  C              "lopped" a cell is (dimensionless scale factor).  C              "lopped" a cell is (dimensionless scale factor).
322  C              Note: The code needs terms like MIN(hFac,hFac(I+1))  C              Note: The code needs terms like MIN(hFac,hFac(I+1))
323  C                    On some platforms it may be better to precompute  C                    On some platforms it may be better to precompute
324  C                    hFacW, hFacE, ... here than do MIN on the fly.  C                    hFacW, hFacE, ... here than do MIN on the fly.
325    C     gravitySign - indicates whether gravity points in the opposite
326    C                   direction of R or not.
327    C                 ( = +1 for R=Z (gravity points downward in Z)
328    C                 ( = -1 for R=P (gravity points upward in P)
329  C     rkFac     - Vertical coordinate to vertical index orientation.  C     rkFac     - Vertical coordinate to vertical index orientation.
330  C                 ( -1 same orientation, 1 opposite orientation )  C                 ( -1 same orientation, 1 opposite orientation )
331  C                 ( vertical coord == m  -> rkFac =  1 )  C                 ( vertical coord == m  -> rkFac =  1 )
332  C                 ( vertical coord == Pa -> rkFac = -1 )  C                 ( vertical coord == Pa -> rkFac = -1 )
333    C     maskH  - cell Center full-column mask (= 2D mask)
334  C     maskC  - cell Center land mask  C     maskC  - cell Center land mask
335  C     maskW  - West face land mask  C     maskW  - West face land mask
336  C     maskS  - South face land mask  C     maskS  - South face land mask
# Line 361  C     tanPhiAtU - tan of the latitude at Line 369  C     tanPhiAtU - tan of the latitude at
369  C                 metric term in U equation.  C                 metric term in U equation.
370  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
371  C                 metric term in V equation.  C                 metric term in V equation.
372          COMMON /GRID_I/ klowC
373          INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
374    
375        COMMON /GRID_R/        COMMON /GRID_R/
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,
# Line 370  C                 metric term in V equat Line 381  C                 metric term in V equat
381       &  recip_hFacC,recip_hFacW,recip_hFacS,       &  recip_hFacC,recip_hFacW,recip_hFacS,
382       &  saFac,       &  saFac,
383       &  xC,yC,rA,rAw,rAs,rAz,xG,yG,       &  xC,yC,rA,rAw,rAs,rAz,xG,yG,
384       &  maskC,maskW,maskS,recip_rA,recip_rAw,recip_rAs,recip_rAz,       &  maskH, maskC,maskW,maskS,
385         &  recip_rA,recip_rAw,recip_rAs,recip_rAz,
386       &  tanPhiAtU, tanPhiAtV,       &  tanPhiAtU, tanPhiAtV,
387       &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,       &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,
388       &  drC,drF,recip_drC,recip_drF,rC,rF,       &  drC,drF,recip_drC,recip_drF,rC,rF,
389         &  gravitySign,
390       &  rkFac, recip_rkFac, xC0, yC0       &  rkFac, recip_rkFac, xC0, yC0
391        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
392        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 412  C                 metric term in V equat Line 425  C                 metric term in V equat
425        _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)
426        _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)
427        _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)
428          _RS maskH          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
429        _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)
430        _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)
431        _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)
# Line 428  C                 metric term in V equat Line 442  C                 metric term in V equat
442        _RS saFac          (1:Nr)        _RS saFac          (1:Nr)
443        _RS rC             (1:Nr)        _RS rC             (1:Nr)
444        _RS rF             (1:Nr+1)        _RS rF             (1:Nr+1)
445          _RL gravitySign
446        _RS rkFac        _RS rkFac
447        _RS recip_rkFac        _RS recip_rkFac
448        _RS xC0        _RS xC0
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22