/[MITgcm]/MITgcm/model/inc/SURFACE.h
ViewVC logotype

Diff of /MITgcm/model/inc/SURFACE.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by jmc, Tue Dec 10 02:47:26 2002 UTC revision 1.13 by heimbach, Thu Dec 8 15:44:33 2005 UTC
# Line 21  C--   COMMON /SURF_FIXED/  fixed surface Line 21  C--   COMMON /SURF_FIXED/  fixed surface
21  C     Bo_surf  :: Boyancy|1/rho [ocean|atmos] at surface level [=g|alpha(p_o)]  C     Bo_surf  :: Boyancy|1/rho [ocean|atmos] at surface level [=g|alpha(p_o)]
22  C     recip_Bo :: 1/Bo_surf  C     recip_Bo :: 1/Bo_surf
23  C     topoZ    :: topographic height [m] (used mainly for atmosphere)  C     topoZ    :: topographic height [m] (used mainly for atmosphere)
24  C     phi0surf :: starting point for integrating phi'_Hyd  C     phi0surf :: starting point for integrating phi_Hyd
25        COMMON /SURF_FIXED/ Bo_surf, recip_Bo, topoZ, phi0surf        COMMON /SURF_FIXED/ Bo_surf, recip_Bo, topoZ, phi0surf
26        _RL  Bo_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  Bo_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
27        _RL  recip_Bo(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  recip_Bo(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 39  C IMPORTANT:  ksurfC,W,S = Nr+1  where t Line 39  C IMPORTANT:  ksurfC,W,S = Nr+1  where t
39        INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40    
41  #ifdef EXACT_CONSERV  #ifdef EXACT_CONSERV
42  C     hDivFlow :: Div. Barotropic Flow at current time [transport unit m3/s]  C     etaHnm1 :: surface r-anomaly, etaH, at previous time level
43        COMMON /EXACT_ETA_LOCAL/ hDivFlow  C     dEtaHdt :: time derivative of total column height [r_unit/s = w unit]
44        _RL hDivFlow(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)  C     PmEpR   :: keep the fresh water input (=-EmPmR) of the previous time step
45          COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR
46          _RL etaHnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
47          _RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
48          _RS  PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
49  #endif  #endif
50    
51  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
52  C--   COMMON /SURF_CHANGE/ transient varariable used for Non-Lin Free-Surf  C--   COMMON /SURF_CHANGE/ transient variables used for Non-Lin Free-Surf
53  C     hFac_surfC ::  New thickness factor of the surface level  C     hFac_surfC ::  New thickness factor of the surface level
54  C                        center (Tracer point)  C                        center (Tracer point)
55  C     hFac_surfW ::  idem, West  interface (U point)  C     hFac_surfW ::  idem, West  interface (U point)
56  C     hFac_surfS ::  idem, South interface (V point)  C     hFac_surfS ::  idem, South interface (V point)
 C     PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step  
57        COMMON /SURF_CHANGE/        COMMON /SURF_CHANGE/
58       &     hFac_surfC, hFac_surfW, hFac_surfS,       &     hFac_surfC, hFac_surfW, hFac_surfS
      &     PmEpR  
59        _RS  hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
60        _RS  hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
61        _RS  hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
62        _RS  PmEpR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
63    C     Local variables in common block
64    C     Rmin_surf :: minimum r_value of the free surface position
65    C                  that satisfy  the hFacInf criteria
66          COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
67          _RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
68    
69    C--   COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate
70    C     rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H
71    C     rStarFacW :: same but for West  face
72    C     rStarFacS :: same but for South face
73    C     rStarExpC :: column expansion factor = h^n+1/h^n , Centered
74    C     rStarExpW :: column expansion factor = h^n+1/h^n , Western  face
75    C     rStarExpS :: column expansion factor = h^n+1/h^n , Southern face
76    C     rStarDhCDt:: relative time derivative of h_Center = d.eta/dt / H
77    C     rStarDhWDt:: relative time derivative of h_West_face  (u.point)
78    C     rStarDhSDt:: relative time derivative of h_South_face (v.point)
79          COMMON /RSTAR_CHANGE/
80         &     rStarFacC, rStarFacW, rStarFacS,
81         &     rStarExpC, rStarExpW, rStarExpS,
82         &     rStarDhCDt,rStarDhWDt,rStarDhSDt
83          _RL  rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84          _RL  rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85          _RL  rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86          _RL  rStarExpC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87          _RL  rStarExpW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88          _RL  rStarExpS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89          _RL  rStarDhCDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90          _RL  rStarDhWDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91          _RL  rStarDhSDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92    
93    C--   COMMON /RSTAR_FIXED/ fixed thickness ratio ( r* discretization )
94    C     h0FacC :: initial (and fixed in time) hFacC factor
95    C     h0FacW :: initial (and fixed in time) hFacW factor
96    C     h0FacS :: initial (and fixed in time) hFacS factor
97          COMMON /RSTAR_FIXED/
98         & h0FacC, h0FacW, h0FacS
99          _RS h0FacC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
100          _RS h0FacW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
101          _RS h0FacS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
102    
103  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22