/[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.16 by cnh, Fri Sep 21 15:13:31 2001 UTC revision 1.22 by edhill, Thu Oct 9 04:19:18 2003 UTC
# Line 104  C     | Line 104  C     |
104  C     |      "WB"                                                  C     |      "WB"                                                
105  C     |       +                                                    C     |       +                                                  
106  C     |       +                                                    C     |       +                                                  
107  C     |      \+/       /|\                       /|\                C     |      \+/       /|\                       /|\       .
108  C     |"UB"++>|-------- | -----------------etc..  | ----*---        C     |"UB"++>|-------- | -----------------etc..  | ----*---      
109  C     |       |    rVel(i=1,        |        rVel(i=Nx, *  |        C     |       |    rVel(i=1,        |        rVel(i=Nx, *  |      
110  C     |       |         j=1,        |             j=1,  *  |        C     |       |         j=1,        |             j=1,  *  |      
# Line 145  C     |"LB"++>========================== Line 145  C     |"LB"++>==========================
145  C     |                                               "PWX"        C     |                                               "PWX"      
146  C     |                                                            C     |                                                          
147  C     | Up   increasing upwards.                                    C     | Up   increasing upwards.                                  
148  C     |/|\                                                          C     |/|\                                                       .
149  C     | |                                                          C     | |                                                        
150  C     | |                                                          C     | |                                                        
151  C     | =====> E  i increasing eastwards                            C     | =====> E  i increasing eastwards                          
# Line 231  C     |  |    |                       | Line 231  C     |  |    |                       |
231  C     | \|/   |           |<---dxV(i=2,j=1,k=1)--\|/         |      C     | \|/   |           |<---dxV(i=2,j=1,k=1)--\|/         |    
232  C     |"SB"++>|___________v___________|___________v__________|_    C     |"SB"++>|___________v___________|___________v__________|_  
233  C     |       <--dxG(i=1,j=1,k=1)----->                            C     |       <--dxG(i=1,j=1,k=1)----->                          
234  C     |      /+\                                                    C     |      /+\                                              .
235  C     |       +                                                    C     |       +                                                  
236  C     |       +                                                    C     |       +                                                  
237  C     |     "WB"                                                    C     |     "WB"                                                  
# Line 287  C     |       |           |           | Line 287  C     |       |           |           |
287  C     |       |           |           |                            C     |       |           |           |                          
288  C     |       |          \|/          |                            C     |       |          \|/          |                          
289  C     |"SB"++>|___________v___________|___________v__etc...        C     |"SB"++>|___________v___________|___________v__etc...      
290  C     |      /+\                                                    C     |      /+\                                                 .
291  C     |       +                                                    C     |       +                                                  
292  C     |       +                                                    C     |       +                                                  
293  C     |     "WB"                                                    C     |     "WB"                                                  
# 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.16  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22