1 |
jmc |
1.5 |
C $Header: /u/gcmpack/MITgcm/verification/tutorial_deep_convection/code/MOM_VISC.h,v 1.1 2013/09/30 18:19:41 jmc Exp $ |
2 |
dfer |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
jmc |
1.3 |
C- Common file for length scales |
5 |
dfer |
1.1 |
|
6 |
|
|
#ifdef ALLOW_MOM_COMMON |
7 |
|
|
|
8 |
jmc |
1.4 |
C-- COMMON /MOM_VISC_PAR_L/ logical-type parameters for Momemtum viscosity |
9 |
|
|
C useHarmonicVisc :: harmonic horizontal viscosity is used |
10 |
|
|
C useBiharmonicVisc :: biharmonic horizontal viscosity is used |
11 |
|
|
C useVariableVisc :: variable (in space or time) viscosity is used |
12 |
|
|
COMMON /MOM_VISC_PAR_L/ |
13 |
|
|
& useHarmonicVisc, useBiharmonicVisc, useVariableVisc |
14 |
|
|
LOGICAL useHarmonicVisc, useBiharmonicVisc, useVariableVisc |
15 |
|
|
|
16 |
|
|
COMMON /MOM_VISC_LENGTH/ L2_D, L2_Z, |
17 |
|
|
& L3_D, L3_Z, |
18 |
|
|
& L4rdt_D, L4rdt_Z |
19 |
dfer |
1.1 |
_RL L2_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
20 |
|
|
_RL L2_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
21 |
|
|
_RL L3_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
22 |
|
|
_RL L3_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
23 |
|
|
_RL L4rdt_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
24 |
|
|
_RL L4rdt_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
25 |
|
|
|
26 |
jmc |
1.5 |
#ifdef ALLOW_SMAG_3D |
27 |
|
|
C smag3D_hLsC :: horiz. grid length scale (power 2/3) at grid cell center |
28 |
|
|
C smag3D_hLsW :: horiz. grid length scale (power 2/3) at western edge |
29 |
|
|
C smag3D_hLsS :: horiz. grid length scale (power 2/3) at southern egde |
30 |
|
|
C smag3D_hLsZ :: horiz. grid length scale (power 2/3) at grid cell corner |
31 |
|
|
COMMON /MOM_SMAG_3D_LENGTH/ |
32 |
|
|
& smag3D_hLsC, smag3D_hLsW, |
33 |
|
|
& smag3D_hLsS, smag3D_hLsZ |
34 |
|
|
_RS smag3D_hLsC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
35 |
|
|
_RS smag3D_hLsW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
36 |
|
|
_RS smag3D_hLsS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
37 |
|
|
_RS smag3D_hLsZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
38 |
|
|
#endif /* ALLOW_SMAG_3D */ |
39 |
|
|
|
40 |
gforget |
1.2 |
#ifdef ALLOW_3D_VISCAH |
41 |
|
|
C viscAhDfld, viscAhZfld :: full 3D specification of Laplacian Viscosity |
42 |
jmc |
1.3 |
C coeff. for mixing of momentum horizontally ( units of m^2/s ) |
43 |
jmc |
1.4 |
COMMON /MOM_VISC_3D_VISCAH/ |
44 |
gforget |
1.2 |
& viscAhDfld, viscAhZfld |
45 |
|
|
_RL viscAhDfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
46 |
|
|
_RL viscAhZfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
47 |
|
|
#endif |
48 |
|
|
|
49 |
|
|
#ifdef ALLOW_3D_VISCA4 |
50 |
|
|
C viscA4Dfld, viscA4Zfld :: full 3D specification of Bi-harmonic Viscosity |
51 |
jmc |
1.3 |
C coeff. for mixing of momentum horizontally ( units of m^4/s ) |
52 |
jmc |
1.4 |
COMMON /MOM_VISC_3D_VISCA4/ |
53 |
gforget |
1.2 |
& viscA4Dfld, viscA4Zfld |
54 |
|
|
_RL viscA4Dfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
55 |
|
|
_RL viscA4Zfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
56 |
|
|
#endif |
57 |
|
|
|
58 |
jmc |
1.4 |
#ifdef ALLOW_NONHYDROSTATIC |
59 |
|
|
C viscAh_W :: Horizontal harmonic viscosity for vertical momentum |
60 |
|
|
C viscA4_W :: Horizontal biharmonic viscosity for vertical momentum |
61 |
|
|
COMMON /MOM_VISC_NH/ |
62 |
|
|
& viscAh_W, viscA4_W |
63 |
|
|
_RL viscAh_W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
64 |
|
|
_RL viscA4_W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
65 |
|
|
#endif /* ALLOW_NONHYDROSTATIC */ |
66 |
|
|
|
67 |
dfer |
1.1 |
#endif /* ALLOW_MOM_COMMON */ |