/[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.1.2.1 by jmc, Fri Mar 30 22:45:06 2001 UTC revision 1.20 by jmc, Tue Apr 29 21:03:00 2014 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  :: Buoyancy|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          _RL  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    C     hFac_surfNm1C, etc. :: prior values
57          COMMON /SURF_CHANGE/
58         &     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)
61          _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)
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
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 /SURF_INDEX/ Common block for surface related index  C--   COMMON /SIGMA_CHANGE/ transient variables used with r* Coordinate
115  C     k_surf - vertical index of the surface tracer cell  C     etaHw    :: surface r-anomaly (etaH) at Western  edge (U location)
116        COMMON /SURF_INDEX/ k_surf  C     etaHs    :: surface r-anomaly (etaH) at Southern edge (V location)
117        INTEGER k_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  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 */
Legend:
Removed from v.1.1.2.1  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22