/[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.11 by jmc, Tue Jun 29 22:21:07 2004 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 /SURF_FIXED/  fixed surface arrays (Real)  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  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
# Line 28  C     phi0surf :: starting point for int Line 28  C     phi0surf :: starting point for int
28        _RS  topoZ   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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     etaHnm1 :: surface r-anomaly, etaH, at previous time level
42  C     dEtaHdt :: time derivative of total column height [r_unit/s = w unit]  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  C     PmEpR   :: keep the fresh water input (=-EmPmR) of the previous time step
44        COMMON /EXACT_ETA_LOCAL/ dEtaHdt, PmEpR        COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR
45        _RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)        _RL etaHnm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46        _RS  PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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
# Line 52  C     hFac_surfC ::  New thickness facto Line 53  C     hFac_surfC ::  New thickness facto
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     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         &     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  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  C--   COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate
74  C     rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H  C     rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H
75  C     rStarFacW :: same but for West  face  C     rStarFacW :: same but for West  face
76  C     rStarFacS :: same but for South face  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  C     rStarExpC :: column expansion factor = h^n+1/h^n , Centered
80  C     rStarExpW :: column expansion factor = h^n+1/h^n , Western  face  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  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  C     rStarDhCDt:: relative time derivative of h_Center = d.eta/dt / H
83  C     rStarDhWDt:: relative time derivative of h_West_face  (u.point)  C     rStarDhWDt:: relative time derivative of h_West_face  (u.point)
84  C     rStarDhSDt:: relative time derivative of h_South_face (v.point)  C     rStarDhSDt:: relative time derivative of h_South_face (v.point)
85        COMMON /RSTAR_CHANGE/        COMMON /RSTAR_CHANGE/
86       &     rStarFacC, rStarFacW, rStarFacS,       &     rStarFacC, rStarFacW, rStarFacS, pStarFacK,
87         &     rStarFacNm1C, rStarFacNm1W, rStarFacNm1S,
88       &     rStarExpC, rStarExpW, rStarExpS,       &     rStarExpC, rStarExpW, rStarExpS,
89       &     rStarDhCDt,rStarDhWDt,rStarDhSDt       &     rStarDhCDt,rStarDhWDt,rStarDhSDt
90        _RL  rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _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)        _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)        _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)        _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)        _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)        _RL  rStarExpS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 92  C     h0FacS :: initial (and fixed in ti Line 111  C     h0FacS :: initial (and fixed in ti
111        _RS h0FacW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _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)        _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.11  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22