/[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.24 by edhill, Fri Oct 31 20:35:32 2003 UTC revision 1.32 by jmc, Sun Jul 23 23:32:50 2006 UTC
# Line 313  C     drC    - Cell center separation al Line 313  C     drC    - Cell center separation al
313  C     drF    - Cell face separation along Z axis ( units of r ).  C     drF    - Cell face separation along Z axis ( units of r ).
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     xx_r_low - in TAF-sense active replacement of R_low
317  C     Ro_surf- surface reference (at rest) position, r_unit.  C     Ro_surf- surface reference (at rest) position, r_unit.
318  C     klowC  - index of the lowest "wet cell" (2D)  C     klowC  - index of the lowest "wet cell" (2D)
319  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
# Line 323  C                    hFacW, hFacE, ... h Line 324  C                    hFacW, hFacE, ... h
324  C     gravitySign - indicates the direction of gravity relative to R direction  C     gravitySign - indicates the direction of gravity relative to R direction
325  C                  (= -1 for R=Z (Z increases upward, -gravity direction  )  C                  (= -1 for R=Z (Z increases upward, -gravity direction  )
326  C                  (= +1 for R=P (P increases downward, +gravity direction)  C                  (= +1 for R=P (P increases downward, +gravity direction)
327  C     rkFac     - Vertical coordinate to vertical index orientation.  C     rkSign - Vertical coordinate to vertical index orientation.
328  C                 ( -1 same orientation, 1 opposite orientation )  C                 ( +1 same orientation, -1 opposite orientation )
329    C     globalArea :: Domain Integrated horizontal Area [m2]
330  C     maskH  - cell Center full-column mask (= 2D mask)  C     maskH  - cell Center full-column mask (= 2D mask)
331  C     maskC  - cell Center land mask  C     maskC  - cell Center land mask
332  C     maskW  - West face land mask  C     maskW  - West face land mask
# Line 350  C              averaging but are just a Line 352  C              averaging but are just a
352  C              diagnostics etc.. As such xc is in m for cartesian  C              diagnostics etc.. As such xc is in m for cartesian
353  C              coordinates but degrees for spherical polar.  C              coordinates but degrees for spherical polar.
354  C     yC     - Y-coordinate of center of cell f[X,Y].  C     yC     - Y-coordinate of center of cell f[X,Y].
355  C     yG     - Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y].  C     yG     - Y-coordinate of corner of cell ( c-grid vorticity point) f[X,Y].
356  C     xC0, yC0 - West edge x coord  ( metres or degrees )  C     xC0, yC0 - West edge x coord  ( metres or degrees )
357  C                South edge y coord ( metres or degrees )  C                South edge y coord ( metres or degrees )
358  C     rA     - R-face are f[X,Y] ( m^2 ).  C     rA     - R-face are f[X,Y] ( m^2 ).
# Line 364  C     tanPhiAtU - tan of the latitude at Line 366  C     tanPhiAtU - tan of the latitude at
366  C                 metric term in U equation.  C                 metric term in U equation.
367  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
368  C                 metric term in V equation.  C                 metric term in V equation.
369    C     angleCosC :: cosine of grid orientation angle relative to Geographic direction
370    C                  at cell center: alpha=(Eastward_dir,grid_uVel_dir)=(North_d,vVel_d)
371    C     angleSinC :: sinus  of grid orientation angle relative to Geographic direction
372    C                  at cell center: alpha=(Eastward_dir,grid_uVel_dir)=(North_d,vVel_d)
373    C     fCori     :: Coriolis parameter at grid Center point
374    C     fCoriG    :: Coriolis parameter at grid Corner point
375    C     fCoriCos  :: Coriolis Cos(phi) parameter at grid Center point (for NH)
376        COMMON /GRID_I/ klowC        COMMON /GRID_I/ klowC
377        INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
378    
379        COMMON /GRID_R/        COMMON /GRID_R/
380         &  cosfacU,cosfacV,sqCosfacU,sqCosfacV,
381         &  gravitySign, rkSign, globalArea,
382       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,       &  dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU,
383       &  R_low,Ro_surf,HFacC,HFacW,HFacS,       &  R_low,Ro_surf,hFacC,hFacW,hFacS,
384    #ifdef ALLOW_DEPTH_CONTROL
385         &  xx_r_low,
386    #endif /* ALLOW_DEPTH_CONTROL */
387       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,       &  recip_dxC,recip_dxF,recip_dxG,recip_dxV,
388       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,       &  recip_dyC,recip_dyF,recip_dyG,recip_dyU,
389       &  recip_Rcol,       &  recip_Rcol,
390       &  recip_hFacC,recip_hFacW,recip_hFacS,       &  recip_hFacC,recip_hFacW,recip_hFacS,
391       &  saFac,       &  saFac,
392       &  xC,yC,rA,rAw,rAs,rAz,xG,yG,       &  xC,yC,rA,rAw,rAs,rAz,xG,yG,
393       &  maskH, maskC,maskW,maskS,       &  maskH, maskC,maskW,maskS,
394       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,       &  recip_rA,recip_rAw,recip_rAs,recip_rAz,
395       &  tanPhiAtU, tanPhiAtV,       &  tanPhiAtU, tanPhiAtV, angleCosC, angleSinC,
      &  cosfacU,cosfacV,sqcosfacU,sqcosfacV,  
396       &  drC,drF,recip_drC,recip_drF,rC,rF,       &  drC,drF,recip_drC,recip_drF,rC,rF,
397       &  gravitySign,       &  xC0, yC0,
398       &  rkFac, recip_rkFac, xC0, yC0       &  fCori, fCoriG, fCoriCos
399          _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)
400          _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)
401          _RL sqCosfacU(1-Oly:sNy+Oly,nSx,nSy)
402          _RL sqCosfacV(1-Oly:sNy+Oly,nSx,nSy)
403          _RL gravitySign
404          _RL rkSign
405          _RL globalArea
406        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxC            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
407        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dxF            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
408        _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 393  C                 metric term in V equat Line 413  C                 metric term in V equat
413        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS dyU            (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
414        _RS R_low          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS R_low          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
415        _RS Ro_surf        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS Ro_surf        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
416        _RS HFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)  #ifdef ALLOW_DEPTH_CONTROL
417        _RS HFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)        _RL xx_r_low       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
418        _RS HFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)  #endif /* ALLOW_DEPTH_CONTROL */
419          _RS hFacC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
420          _RS hFacW          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
421          _RS hFacS          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)
422        _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)
423        _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)
424        _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 426  C                 metric term in V equat Line 449  C                 metric term in V equat
449        _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)
450        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtU      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
451        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS tanPhiAtV      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
452        _RL cosfacU(1-Oly:sNy+Oly,nSx,nSy)        _RS angleCosC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
453        _RL cosfacV(1-Oly:sNy+Oly,nSx,nSy)        _RS angleSinC      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy)  
       _RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy)  
454        _RS drC            (1:Nr)        _RS drC            (1:Nr)
455        _RS drF            (1:Nr)        _RS drF            (1:Nr)
456        _RS recip_drC      (1:Nr)        _RS recip_drC      (1:Nr)
# Line 437  C                 metric term in V equat Line 458  C                 metric term in V equat
458        _RS saFac          (1:Nr)        _RS saFac          (1:Nr)
459        _RS rC             (1:Nr)        _RS rC             (1:Nr)
460        _RS rF             (1:Nr+1)        _RS rF             (1:Nr+1)
       _RL gravitySign  
       _RS rkFac  
       _RS recip_rkFac  
461        _RS xC0        _RS xC0
462        _RS yC0        _RS yC0
463          _RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
464          _RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
465          _RS fCoriCos(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
466    
 #ifdef ALLOW_NONHYDROSTATIC  
       COMMON /GRID_NH/  
      &  recip_hFacU  
       _RS recip_hFacU    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy)  
 #endif  
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22