/[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.10 by cnh, Tue Sep 29 18:50:56 1998 UTC revision 1.15 by adcroft, Tue May 29 14:01:36 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3  C  C
4  C     /==========================================================\  C     /==========================================================\
5  C     | GRID.h                                                   |  C     | GRID.h                                                   |
# Line 303  C     dyF    - Cell face separation in Y Line 304  C     dyF    - Cell face separation in Y
304  C     dyU    - U-point separation in Y across south-west corner of cell (m)  C     dyU    - U-point separation in Y across south-west corner of cell (m)
305  C     drC    - Cell center separation along Z axis ( units of r ).  C     drC    - Cell center separation along Z axis ( units of r ).
306  C     drF    - Cell face separation along Z axis ( units of r ).  C     drF    - Cell face separation along Z axis ( units of r ).
307  C     H      - Depth of base of fluid from upper surface f[X,Y] (m).  C     Rcolumn  -Total thickness (in r_unit) of the fluid column
308    C     R_low  - base of fluid in r_unit (Depth(m) / Pressure(Pa) at top Atmos.)
309    C     Ro_surf- surface reference (at rest) position, r_unit.
310  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
311  C              "lopped" a cell is (dimensionless scale factor).  C              "lopped" a cell is (dimensionless scale factor).
312  C              Note: The code needs terms like MIN(hFac,hFac(I+1))  C              Note: The code needs terms like MIN(hFac,hFac(I+1))
# Line 313  C     rkFac     - Vertical coordinate to Line 316  C     rkFac     - Vertical coordinate to
316  C                 ( -1 same orientation, 1 opposite orientation )  C                 ( -1 same orientation, 1 opposite orientation )
317  C                 ( vertical coord == m  -> rkFac =  1 )  C                 ( vertical coord == m  -> rkFac =  1 )
318  C                 ( vertical coord == Pa -> rkFac = -1 )  C                 ( vertical coord == Pa -> rkFac = -1 )
319    C     maskC  - cell Center land mask
320  C     maskW  - West face land mask  C     maskW  - West face land mask
321  C     maskS  - South face land mask  C     maskS  - South face land mask
322  C     recip_dxC   - Recipricol of dxC  C     recip_dxC   - Recipricol of dxC
# Line 325  C     recip_dyF   - Recipricol of dyF Line 329  C     recip_dyF   - Recipricol of dyF
329  C     recip_dyU   - Recipricol of dyU  C     recip_dyU   - Recipricol of dyU
330  C     recip_drC   - Recipricol of drC  C     recip_drC   - Recipricol of drC
331  C     recip_drF   - Recipricol of drF  C     recip_drF   - Recipricol of drF
332  C     recip_H     - Inverse of cell center open-depth  C     recip_Rcol  - Inverse of cell center column thickness (1/r_unit)
333  C     recip_hFacC - Inverse of cell open-depth f[X,Y,Z] ( dimensionless ).  C     recip_hFacC - Inverse of cell open-depth f[X,Y,Z] ( dimensionless ).
334  C     recip_hFacW   rhFacC center, rhFacW west, rhFacS south.  C     recip_hFacW   rhFacC center, rhFacW west, rhFacS south.
335  C     recip_hFacS   Note: This is precomputed here because it involves division.  C     recip_hFacS   Note: This is precomputed here because it involves division.
# Line 336  C              averaging but are just a Line 340  C              averaging but are just a
340  C              diagnostics etc.. As such xc is in m for cartesian  C              diagnostics etc.. As such xc is in m for cartesian
341  C              coordinates but degrees for spherical polar.  C              coordinates but degrees for spherical polar.
342  C     yC     - Y-coordinate of center of cell f[X,Y].  C     yC     - Y-coordinate of center of cell f[X,Y].
343    C     yG     - Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y].
344  C     xC0, yC0 - West edge x coord  ( metres or degrees )  C     xC0, yC0 - West edge x coord  ( metres or degrees )
345  C                South edge y coord ( metres or degrees )  C                South edge y coord ( metres or degrees )
346  C     rA     - R-face are f[X,Y] ( m^2 ).  C     rA     - R-face are f[X,Y] ( m^2 ).
# Line 350  C                 metric term in U equat Line 355  C                 metric term in U equat
355  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
356  C                 metric term in V equation.  C                 metric term in V equation.
357        COMMON /GRID_R/        COMMON /GRID_R/
358       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,drC,drF,       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,
359       &  H,HFacC,HFacW,HFacS,DepthInK,       &  R_low,Ro_surf,HFacC,HFacW,HFacS,
360       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,
361       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,
362       &  recip_drC,recip_drF,       &  recip_Rcol,
      &  recip_H,  
363       &  recip_hFacC,recip_hFacW,recip_hFacS,       &  recip_hFacC,recip_hFacW,recip_hFacS,
      &  rkFac, recip_rkFac,  
364       &  saFac,       &  saFac,
365       &  xC,yC,rA,rC,rF,yC0,xC0,       &  xC,yC,rA,rAw,rAs,rAz,xG,yG,
366       &  maskW,maskS,       &  maskC,maskW,maskS,recip_rA,recip_rAw,recip_rAs,recip_rAz,
367       &  tanPhiAtU, tanPhiAtV       &  tanPhiAtU, tanPhiAtV,
368         &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,
369         &  drC,drF,recip_drC,recip_drF,rC,rF,
370         &  rkFac, recip_rkFac, xC0, yC0
371        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
372        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
373        _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 370  C                 metric term in V equat Line 376  C                 metric term in V equat
376        _RS dyF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
377        _RS dyG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
378        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
379        _RS drC            (1:Nr)        _RS R_low          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
380        _RS drF            (1:Nr)        _RS Ro_surf        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RS rkFac  
       _RS recip_rkFac  
       _RS DepthInK       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
       _RS H              (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
381        _RS HFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS HFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
382        _RS HFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS HFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
383        _RS HFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS HFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
# Line 387  C                 metric term in V equat Line 389  C                 metric term in V equat
389        _RS recip_dyF      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dyF      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
390        _RS recip_dyG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dyG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
391        _RS recip_dyU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dyU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
392        _RS recip_drC      (1:Nr)        _RS recip_Rcol     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RS recip_drF      (1:Nr)  
       _RS recip_h        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
393        _RS recip_hFacC    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS recip_hFacC    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
394        _RS recip_hFacW    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS recip_hFacW    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
395        _RS recip_hFacS    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RS recip_hFacS    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
       _RS saFac          (1:Nr)  
396        _RS xC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS xC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
397        _RS xC0        _RS xG             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
398        _RS yC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS yC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
399        _RS yC0        _RS yG             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
400        _RS rA             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS rA             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
401        _RS rC             (1:Nr)        _RS rAw            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
402        _RS rF             (1:Nr+1)        _RS rAs            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
403          _RS rAz            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
404          _RS recip_rA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
405          _RS recip_rAw      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
406          _RS recip_rAs      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
407          _RS recip_rAz      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
408          _RS maskC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
409        _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)
410        _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)
411        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
412        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
413          _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)
414          _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)
415          _RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy)
416          _RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy)
417          _RS drC            (1:Nr)
418          _RS drF            (1:Nr)
419          _RS recip_drC      (1:Nr)
420          _RS recip_drF      (1:Nr)
421          _RS saFac          (1:Nr)
422          _RS rC             (1:Nr)
423          _RS rF             (1:Nr+1)
424          _RS rkFac
425          _RS recip_rkFac
426          _RS xC0
427          _RS yC0
428    
429    #ifdef ALLOW_NONHYDROSTATIC
430          COMMON /GRID_NH/
431         &  recip_hFacU
432          _RS recip_hFacU    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
433    #endif
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22