/[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.2 by adcroft, Tue May 29 14:01:36 2001 UTC revision 1.16 by jmc, Sat Sep 11 21:23:09 2010 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  C  C
4  C     /==========================================================\  CBOP
5  C     | SURFACE.h                                                |  C    !ROUTINE: SURFACE.h
6  C     | o Header file defining surface-related model varaibles   |  C    !INTERFACE:
7  C     |==========================================================|  C    include SURFACE.h
8  C     | Contains variables relative to the surface position      |  C    !DESCRIPTION: \bv
9  C     | that are held fixed in linear free-surface formulation   |  C     *==========================================================*
10  C     | but can vary with time with a non-linear free-surface.   |  C     | SURFACE.h
11  C     \==========================================================/  C     | o Header file defining surface-related model variables
12    C     *==========================================================*
13  C--   COMMON /SOLVE_BAROT/  Barotropic variables common block  C     | Contains variables relative to the surface position
14  C     Bo_surf  -Boyancy|1/rho [ocean|atmos] at surface level [= g | alpha(p_o)]  C     | that are held fixed in linear free-surface formulation
15  C     recip_Bo     = 1/Bo_surf  C     | but can vary with time with a non-linear free-surface.
16        COMMON /SOLVE_BAROT/ Bo_surf, recip_Bo  C     *==========================================================*
17        _RL  Bo_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  C     \ev
18    CEOP
19    
20    C--   COMMON /SURF_FIXED/  fixed surface arrays (Real)
21    C     Bo_surf  :: Boyancy|1/rho [ocean|atmos] at surface level [=g|alpha(p_o)]
22    C     recip_Bo :: 1/Bo_surf
23    C     topoZ    :: topographic height [m] (used mainly for atmosphere)
24    C     phi0surf :: starting point for integrating phi_Hyd
25          COMMON /SURF_FIXED/ Bo_surf, recip_Bo, topoZ, phi0surf
26          _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)
28          _RS  topoZ   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
29          _RS  phi0surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
30    
31    C--   COMMON /SURF_CORREC/ Common block for correction of source/sink of
32    C--                        Tracer due to W at the surface with Linear
33    C--                        Free Surface
34    C     TsurfCor :: Pot.Temp Linear-Free-Surface correction term [K.r_Unit/s]
35    C     SsurfCor :: Salinity Linear-Free-Surface correction term [psu.r_Unit/s]
36          COMMON /SURF_CORREC/ TsurfCor, SsurfCor
37          _RL TsurfCor
38          _RL SsurfCor
39    
40    #ifdef EXACT_CONSERV
41    C     etaHnm1 :: surface r-anomaly, etaH, at previous time level
42    C     dEtaHdt :: time derivative of total column height [r_unit/s = w unit]
43    C     PmEpR   :: keep the fresh water input (=-EmPmR) of the previous time step
44          COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR
45          _RL etaHnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
46          _RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
47          _RS  PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48    #endif
49    
50    #ifdef NONLIN_FRSURF
51    C--   COMMON /SURF_CHANGE/ transient variables used for Non-Lin Free-Surf
52    C     hFac_surfC ::  New thickness factor of the surface level
53    C                        center (Tracer point)
54    C     hFac_surfW ::  idem, West  interface (U point)
55    C     hFac_surfS ::  idem, South interface (V point)
56          COMMON /SURF_CHANGE/
57         &     hFac_surfC, hFac_surfW, hFac_surfS
58          _RS  hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
59          _RS  hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
60          _RS  hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
61    
62    C     Local variables in common block
63    C     Rmin_surf :: minimum r_value of the free surface position
64    C                  that satisfy  the hFacInf criteria
65          COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
66          _RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
67    
68    C--   COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate
69    C     rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H
70    C     rStarFacW :: same but for West  face
71    C     rStarFacS :: same but for South face
72    C     rStarExpC :: column expansion factor = h^n+1/h^n , Centered
73    C     rStarExpW :: column expansion factor = h^n+1/h^n , Western  face
74    C     rStarExpS :: column expansion factor = h^n+1/h^n , Southern face
75    C     rStarDhCDt:: relative time derivative of h_Center = d.eta/dt / H
76    C     rStarDhWDt:: relative time derivative of h_West_face  (u.point)
77    C     rStarDhSDt:: relative time derivative of h_South_face (v.point)
78          COMMON /RSTAR_CHANGE/
79         &     rStarFacC, rStarFacW, rStarFacS,
80         &     rStarExpC, rStarExpW, rStarExpS,
81         &     rStarDhCDt,rStarDhWDt,rStarDhSDt
82          _RL  rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
83          _RL  rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84          _RL  rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85          _RL  rStarExpC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86          _RL  rStarExpW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87          _RL  rStarExpS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88          _RL  rStarDhCDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89          _RL  rStarDhWDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90          _RL  rStarDhSDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91    
92    C--   COMMON /RSTAR_FIXED/ fixed thickness ratio ( r* discretization )
93    C     h0FacC :: initial (and fixed in time) hFacC factor
94    C     h0FacW :: initial (and fixed in time) hFacW factor
95    C     h0FacS :: initial (and fixed in time) hFacS factor
96          COMMON /RSTAR_FIXED/
97         & h0FacC, h0FacW, h0FacS
98          _RS h0FacC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
99          _RS h0FacW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
100          _RS h0FacS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
101    
102  C--   COMMON /SURF_INDEX/ Common block for surface related index  C--   COMMON /SIGMA_CHANGE/ transient variables used with r* Coordinate
103  C     k_surf - vertical index of the surface tracer cell  C     etaHw    :: surface r-anomaly (etaH) at Western  edge (U location)
104        COMMON /SURF_INDEX/ k_surf  C     etaHs    :: surface r-anomaly (etaH) at Southern edge (V location)
105        INTEGER k_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  C     dEtaWdt  :: time derivative of etaH at Western  edge (U location)
106    C     dEtaSdt  :: time derivative of etaH at Southern edge (V location)
107          COMMON /SIGMA_CHANGE/
108         &  etaHw, etaHs,
109         &  dEtaWdt, dEtaSdt
110          _RL  etaHw  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
111          _RL  etaHs  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
112          _RL  dEtaWdt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
113          _RL  dEtaSdt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
114    #endif /* NONLIN_FRSURF */
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22