1 |
C $Header: /u/gcmpack/MITgcm/model/inc/NH_VARS.h,v 1.2 2006/06/20 20:57:37 baylor Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
CBOP |
5 |
C !ROUTINE: NH_VARS.h |
6 |
C !INTERFACE: |
7 |
C include "NH_VARS.h" |
8 |
C !DESCRIPTION: |
9 |
C \bv |
10 |
C *==========================================================* |
11 |
C | NH_VARS.h |
12 |
C | o Additional state variables for non-hydrostatic model |
13 |
C *==========================================================* |
14 |
C | In N-H mode, wVel becomes a prognostic variable: need |
15 |
C | to hold two levels of time tendency for w (for AB) |
16 |
C *==========================================================* |
17 |
C \ev |
18 |
CEOP |
19 |
|
20 |
#ifdef ALLOW_NONHYDROSTATIC |
21 |
C-- COMMON /NH_VARS_R/ REAL NH state variables |
22 |
C phi_nh :: Non-hydrostatic potential (=NH-Pressure/rhoConst) |
23 |
C dPhiNH :: Surface-Hydrostatic correction to Non-hydrostatic Phi |
24 |
C gX, gxNm1 :: Time tendencies at current and previous time levels. |
25 |
C viscA?_W :: Horizontal variable viscosities |
26 |
|
27 |
#ifdef ALLOW_ADAMSBASHFORTH_3 |
28 |
COMMON /NH_VARS_R/ |
29 |
& phi_nh, dPhiNH, |
30 |
& gW, gwNm, |
31 |
& viscAh_W, viscA4_W |
32 |
#else /* ALLOW_ADAMSBASHFORTH_3 */ |
33 |
COMMON /NH_VARS_R/ |
34 |
& phi_nh, dPhiNH, |
35 |
& gW, gwNm1, |
36 |
& viscAh_W, viscA4_W |
37 |
#endif /* ALLOW_ADAMSBASHFORTH_3 */ |
38 |
_RL phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
39 |
_RL dPhiNH(1-OLx:sNx+OLx,1-OLy:sNy+OLy, nSx,nSy) |
40 |
_RL gW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
41 |
#ifdef ALLOW_ADAMSBASHFORTH_3 |
42 |
_RL gwNm (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,2) |
43 |
#else /* ALLOW_ADAMSBASHFORTH_3 */ |
44 |
_RL gwNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
45 |
#endif /* ALLOW_ADAMSBASHFORTH_3 */ |
46 |
_RL viscAh_W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
47 |
_RL viscA4_W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
48 |
|
49 |
#endif /* ALLOW_NONHYDROSTATIC */ |