/[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.18 by jmc, Thu Nov 8 16:34:57 2001 UTC revision 1.26 by heimbach, Sat Jul 24 01:26:26 2004 UTC
# Line 314  C     drF    - Cell face separation alon Line 314  C     drF    - Cell face separation alon
314  C     Rcolumn  -Total thickness (in r_unit) of the fluid column  C     Rcolumn  -Total thickness (in r_unit) of the fluid column
315  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.)
316  C     Ro_surf- surface reference (at rest) position, r_unit.  C     Ro_surf- surface reference (at rest) position, r_unit.
317    C     klowC  - index of the lowest "wet cell" (2D)
318  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
319  C              "lopped" a cell is (dimensionless scale factor).  C              "lopped" a cell is (dimensionless scale factor).
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 the direction of gravity relative to R direction
324    C                  (= -1 for R=Z (Z increases upward, -gravity direction  )
325    C                  (= +1 for R=P (P increases downward, +gravity direction)
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 362  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          COMMON /GRID_I/ klowC
368          INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
369    
370        COMMON /GRID_R/        COMMON /GRID_R/
371         &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,
372         &  gravitySign,rkFac, recip_rkFac,
373       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,
374       &  R_low,Ro_surf,HFacC,HFacW,HFacS,       &  R_low,Ro_surf,HFacC,HFacW,HFacS,
375       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,
# Line 374  C                 metric term in V equat Line 381  C                 metric term in V equat
381       &  maskH, maskC,maskW,maskS,       &  maskH, maskC,maskW,maskS,
382       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,
383       &  tanPhiAtU, tanPhiAtV,       &  tanPhiAtU, tanPhiAtV,
      &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,  
384       &  drC,drF,recip_drC,recip_drF,rC,rF,       &  drC,drF,recip_drC,recip_drF,rC,rF,
385       &  rkFac, recip_rkFac, xC0, yC0       &  xC0, yC0
386          _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)
387          _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)
388          _RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy)
389          _RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy)
390          _RL gravitySign
391          _RL rkFac
392          _RL recip_rkFac
393        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
394        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
395        _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 420  C                 metric term in V equat Line 433  C                 metric term in V equat
433        _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)
434        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
435        _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)  
436        _RS drC            (1:Nr)        _RS drC            (1:Nr)
437        _RS drF            (1:Nr)        _RS drF            (1:Nr)
438        _RS recip_drC      (1:Nr)        _RS recip_drC      (1:Nr)
# Line 431  C                 metric term in V equat Line 440  C                 metric term in V equat
440        _RS saFac          (1:Nr)        _RS saFac          (1:Nr)
441        _RS rC             (1:Nr)        _RS rC             (1:Nr)
442        _RS rF             (1:Nr+1)        _RS rF             (1:Nr+1)
       _RS rkFac  
       _RS recip_rkFac  
443        _RS xC0        _RS xC0
444        _RS yC0        _RS yC0
445    
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22