/[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.4.2.2 by heimbach, Mon Apr 8 20:10:36 2002 UTC revision 1.20 by jmc, Tue Apr 29 21:03:00 2014 UTC
# Line 7  C    !INTERFACE: Line 7  C    !INTERFACE:
7  C    include SURFACE.h  C    include SURFACE.h
8  C    !DESCRIPTION: \bv  C    !DESCRIPTION: \bv
9  C     *==========================================================*  C     *==========================================================*
10  C     | SURFACE.h                                                  C     | SURFACE.h
11  C     | o Header file defining surface-related model varaibles      C     | o Header file defining surface-related model variables
12  C     *==========================================================*  C     *==========================================================*
13  C     | Contains variables relative to the surface position        C     | Contains variables relative to the surface position
14  C     | that are held fixed in linear free-surface formulation      C     | that are held fixed in linear free-surface formulation
15  C     | but can vary with time with a non-linear free-surface.      C     | but can vary with time with a non-linear free-surface.
16  C     *==========================================================*  C     *==========================================================*
17  C     \ev  C     \ev
18  CEOP  CEOP
19    
20  C--   COMMON /SOLVE_BAROT/  Barotropic variables common block  C--   COMMON /SURF_FIXED/  fixed surface arrays (Real)
21  C     Bo_surf  :: Boyancy|1/rho [ocean|atmos] at surface level [= g | alpha(p_o)]  C     Bo_surf  :: Buoyancy|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        COMMON /SOLVE_BAROT/ Bo_surf, recip_Bo  C     topoZ    :: topographic height [m] (used mainly for atmosphere)
24        _RL  Bo_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  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_INDEX/ Common block for surface related index  C--   COMMON /SURF_CORREC/ Common block for correction of source/sink of
32  C     ksurfC ::  vertical index of the surface tracer cell  C--                        Tracer due to W at the surface with Linear
33  C     ksurfW ::  vertical index of the surface U point  C--                        Free Surface
34  C     ksurfS ::  vertical index of the surface V point  C     TsurfCor :: Pot.Temp Linear-Free-Surface correction term [K.r_Unit/s]
35  C IMPORTANT:  ksurfC,W,S = Nr+1  where the fluid column is empty (continent)  C     SsurfCor :: Salinity Linear-Free-Surface correction term [psu.r_Unit/s]
36        COMMON /SURF_INDEX/ ksurfC, ksurfW, ksurfS        COMMON /SURF_CORREC/ TsurfCor, SsurfCor
37        INTEGER ksurfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL TsurfCor
38        INTEGER ksurfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL SsurfCor
       INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
39    
40  #ifdef EXACT_CONSERV  #ifdef EXACT_CONSERV
41  C     hDivFlow :: Div. Barotropic Flow at current time [transport unit m3/s]  C     etaHnm1 :: surface r-anomaly, etaH, at previous time level
42        COMMON /EXACT_ETA_LOCAL/ hDivFlow  C     dEtaHdt :: time derivative of total column height [r_unit/s = w unit]
43        _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
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          _RL  PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48  #endif  #endif
49    
50  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
51  C--   COMMON /SURF_CHANGE/ transient varariable used for Non-Lin Free-Surf  C--   COMMON /SURF_CHANGE/ transient variables used for Non-Lin Free-Surf
52  C     hFac_surfC ::  New thickness factor of the surface level  C     hFac_surfC ::  New thickness factor of the surface level
53  C                        center (Tracer point)  C                        center (Tracer point)
54  C     hFac_surfW ::  idem, West  interface (U point)  C     hFac_surfW ::  idem, West  interface (U point)
55  C     hFac_surfS ::  idem, South interface (V point)  C     hFac_surfS ::  idem, South interface (V point)
56  C     PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step  C     hFac_surfNm1C, etc. :: prior values
57        COMMON /SURF_CHANGE/        COMMON /SURF_CHANGE/
58       &     hFac_surfC, hFac_surfW, hFac_surfS,       &     hFac_surfC, hFac_surfW, hFac_surfS,
59       &     PmEpR       &     hFac_surfNm1C, hFac_surfNm1W, hFac_surfNm1S
60        _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)
61        _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)
62        _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)
63        _RS  PmEpR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  hFac_surfNm1C(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
64          _RS  hFac_surfNm1W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
65          _RS  hFac_surfNm1S(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
66    
67    C     Local variables in common block
68    C     Rmin_surf :: minimum r_value of the free surface position
69    C                  that satisfy  the hFacInf criteria
70          COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
71          _RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
72    
73    C--   COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate
74    C     rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H
75    C     rStarFacW :: same but for West  face
76    C     rStarFacS :: same but for South face
77    C     pStarFacK :: rStarFacC**atm_kappa (for atmosphere in p* coords)
78    C     rStarFacNm1C, etc. :: prior values
79    C     rStarExpC :: column expansion factor = h^n+1/h^n , Centered
80    C     rStarExpW :: column expansion factor = h^n+1/h^n , Western  face
81    C     rStarExpS :: column expansion factor = h^n+1/h^n , Southern face
82    C     rStarDhCDt:: relative time derivative of h_Center = d.eta/dt / H
83    C     rStarDhWDt:: relative time derivative of h_West_face  (u.point)
84    C     rStarDhSDt:: relative time derivative of h_South_face (v.point)
85          COMMON /RSTAR_CHANGE/
86         &     rStarFacC, rStarFacW, rStarFacS, pStarFacK,
87         &     rStarFacNm1C, rStarFacNm1W, rStarFacNm1S,
88         &     rStarExpC, rStarExpW, rStarExpS,
89         &     rStarDhCDt,rStarDhWDt,rStarDhSDt
90          _RL  rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91          _RL  rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92          _RL  rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93          _RL  pStarFacK (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94          _RL  rStarFacNm1C (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
95          _RL  rStarFacNm1W (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
96          _RL  rStarFacNm1S (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
97          _RL  rStarExpC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98          _RL  rStarExpW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
99          _RL  rStarExpS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100          _RL  rStarDhCDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101          _RL  rStarDhWDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
102          _RL  rStarDhSDt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
103    
104    C--   COMMON /RSTAR_FIXED/ fixed thickness ratio ( r* discretization )
105    C     h0FacC :: initial (and fixed in time) hFacC factor
106    C     h0FacW :: initial (and fixed in time) hFacW factor
107    C     h0FacS :: initial (and fixed in time) hFacS factor
108          COMMON /RSTAR_FIXED/
109         & h0FacC, h0FacW, h0FacS
110          _RS h0FacC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
111          _RS h0FacW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
112          _RS h0FacS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
113    
114    C--   COMMON /SIGMA_CHANGE/ transient variables used with r* Coordinate
115    C     etaHw    :: surface r-anomaly (etaH) at Western  edge (U location)
116    C     etaHs    :: surface r-anomaly (etaH) at Southern edge (V location)
117    C     dEtaWdt  :: time derivative of etaH at Western  edge (U location)
118    C     dEtaSdt  :: time derivative of etaH at Southern edge (V location)
119          COMMON /SIGMA_CHANGE/
120         &  etaHw, etaHs,
121         &  dEtaWdt, dEtaSdt
122          _RL  etaHw  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
123          _RL  etaHs  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
124          _RL  dEtaWdt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
125          _RL  dEtaSdt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
126  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
Legend:
Removed from v.1.4.2.2  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22