/[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.5 by cnh, Sat Aug 15 16:55:47 1998 UTC revision 1.13 by adcroft, Fri Feb 2 21:04:47 2001 UTC
# Line 129  C     |       |        /|\          | Line 129  C     |       |        /|\          |
129  C     |       |-------- | -----------------etc..  | ----*---     |  C     |       |-------- | -----------------etc..  | ----*---     |
130  C     |       |    rVel(i=1,        |        rVel(i=Nx, *  |     |  C     |       |    rVel(i=1,        |        rVel(i=Nx, *  |     |
131  C     |       |         j=1,        |             j=1,  *  |     |  C     |       |         j=1,        |             j=1,  *  |     |
132  C     |       |         k=Nz)       |             k=Nz) *  |     |  C     |       |         k=Nr)       |             k=Nr) *  |     |
133  C     |U(i=1, ==>       x         ==>U(i=2,       x     *==>U    |  C     |U(i=1, ==>       x         ==>U(i=2,       x     *==>U    |
134  C     |  j=1, |      T(i=1,         |  j=1,    T(i=Nx,  *(i=Nx+1,|  C     |  j=1, |      T(i=1,         |  j=1,    T(i=Nx,  *(i=Nx+1,|
135  C     |  k=Nz)|        j=1,         |  k=Nz)     j=1,   *  |j=1, |  C     |  k=Nr)|        j=1,         |  k=Nr)     j=1,   *  |j=1, |
136  C     |       |        k=Nz)        |            k=Nz)  *  |k=Nz)|  C     |       |        k=Nr)        |            k=Nr)  *  |k=Nr)|
137  C     |       |                     |                   *  |     |  C     |       |                     |                   *  |     |
138  C     |"LB"++>==============================================     |  C     |"LB"++>==============================================     |
139  C     |                                               "PWX"      |  C     |                                               "PWX"      |
# Line 309  C              "lopped" a cell is (dimen Line 309  C              "lopped" a cell is (dimen
309  C              Note: The code needs terms like MIN(hFac,hFac(I+1))  C              Note: The code needs terms like MIN(hFac,hFac(I+1))
310  C                    On some platforms it may be better to precompute  C                    On some platforms it may be better to precompute
311  C                    hFacW, hFacE, ... here than do MIN on the fly.  C                    hFacW, hFacE, ... here than do MIN on the fly.
312    C     rkFac     - Vertical coordinate to vertical index orientation.
313    C                 ( -1 same orientation, 1 opposite orientation )
314    C                 ( vertical coord == m  -> rkFac =  1 )
315    C                 ( vertical coord == Pa -> rkFac = -1 )
316  C     maskW  - West face land mask  C     maskW  - West face land mask
317  C     maskS  - South face land mask  C     maskS  - South face land mask
318  C     recip_dxC   - Recipricol of dxC  C     recip_dxC   - Recipricol of dxC
# Line 332  C              averaging but are just a Line 336  C              averaging but are just a
336  C              diagnostics etc.. As such xc is in m for cartesian  C              diagnostics etc.. As such xc is in m for cartesian
337  C              coordinates but degrees for spherical polar.  C              coordinates but degrees for spherical polar.
338  C     yC     - Y-coordinate of center of cell f[X,Y].  C     yC     - Y-coordinate of center of cell f[X,Y].
339    C     yG     - Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y].
340  C     xC0, yC0 - West edge x coord  ( metres or degrees )  C     xC0, yC0 - West edge x coord  ( metres or degrees )
341  C                South edge y coord ( metres or degrees )  C                South edge y coord ( metres or degrees )
342  C     rA     - R-face are f[X,Y] ( m^2 ).  C     rA     - R-face are f[X,Y] ( m^2 ).
# Line 346  C                 metric term in U equat Line 351  C                 metric term in U equat
351  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
352  C                 metric term in V equation.  C                 metric term in V equation.
353        COMMON /GRID_R/        COMMON /GRID_R/
354       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,dzC,drF,       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,drC,drF,
355       &  H,HFacC,HFacW,HFacS,       &  H,HFacC,HFacW,HFacS,DepthInK,
356       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,
357       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,
358       &  recip_drC,recip_drF,       &  recip_drC,recip_drF,
359       &  recip_H,       &  recip_H,
360       &  recip_hFacC,recip_hFacW,recip_hFacS,       &  recip_hFacC,recip_hFacW,recip_hFacS,
361         &  rkFac, recip_rkFac,
362       &  saFac,       &  saFac,
363       &  xC,yC,rA,rC,rF,yC0,xC0,       &  xC,yC,rA,rAw,rAs,rAz,rC,rF,yC0,xC0,xG,yG,
364       &  maskW,maskS,       &  maskW,maskS,recip_rA,recip_rAw,recip_rAs,recip_rAz,
365       &  tanPhiAtU, tanPhiAtV       &  tanPhiAtU, tanPhiAtV
366        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
367        _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 365  C                 metric term in V equat Line 371  C                 metric term in V equat
371        _RS dyF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
372        _RS dyG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyG            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
373        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
374        _RS drC            (1:Nz)        _RS drC            (1:Nr)
375        _RS drF            (1:Nz)        _RS drF            (1:Nr)
376          _RS rkFac
377          _RS recip_rkFac
378          _RS DepthInK       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
379        _RS H              (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS H              (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
380        _RS HFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS HFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
381        _RS HFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS HFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
382        _RS HFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS HFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
383        _RS recip_dxC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dxC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
384        _RS recip_dxF      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dxF      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
385        _RS recip_dxG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_dxG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 379  C                 metric term in V equat Line 388  C                 metric term in V equat
388        _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)
389        _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)
390        _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)
391        _RS recip_dzC      (1:Nz)        _RS recip_drC      (1:Nr)
392        _RS recip_drF      (1:Nz)        _RS recip_drF      (1:Nr)
393        _RS recip_h        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS recip_h        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
394        _RS recip_hFacC    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS recip_hFacC    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
395        _RS recip_hFacW    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS recip_hFacW    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
396        _RS recip_hFacS    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS recip_hFacS    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
397        _RS saFac          (1:Nz)        _RS saFac          (1:Nr)
398        _RS xC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS xC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
399          _RS xG             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
400        _RS xC0        _RS xC0
401        _RS yC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS yC             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
402          _RS yG             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
403        _RS yC0        _RS yC0
404        _RS rA             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS rA             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
405        _RS rC             (1:Nz)        _RS rAw            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
406        _RS rF             (1:Nz+1)        _RS rAs            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
407        _RS maskW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS rAz            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
408        _RS maskS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nz,nSx,nSy)        _RS recip_rA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
409          _RS recip_rAw      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
410          _RS recip_rAs      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
411          _RS recip_rAz      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
412          _RS rC             (1:Nr)
413          _RS rF             (1:Nr+1)
414          _RS maskW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
415          _RS maskS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
416        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
417        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
418    
419    #ifdef ALLOW_NONHYDROSTATIC
420  C--   COMMON /GRID_I/ Integer valued grid defining variables        COMMON /GRID_NH/
421  C     iDep  - Index of last "non-land" cell for each column f[X,Y].       &  recip_hFacU
422        COMMON /GRID_I/        _RS recip_hFacU    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
423       &      iDep  #endif
       INTEGER iDep(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22