1 |
C $Header: $ |
2 |
|
3 |
#ifdef ALLOW_GMREDI |
4 |
|
5 |
C GM/Redi parameters |
6 |
_RL GM_background_K |
7 |
_RL GM_maxSlope |
8 |
_RL GM_Visbeck_alpha |
9 |
_RL GM_Visbeck_length |
10 |
_RL GM_Visbeck_depth |
11 |
_RL GM_Visbeck_maxval_K |
12 |
CHARACTER*(MAX_LEN_FNAM) GM_taper_scheme |
13 |
_RL GM_Scrit |
14 |
_RL GM_Sd |
15 |
COMMON /GM_PARAMS/ GM_background_K, |
16 |
& GM_maxSlope, |
17 |
& GM_Visbeck_alpha, |
18 |
& GM_Visbeck_length, |
19 |
& GM_Visbeck_depth, |
20 |
& GM_Visbeck_maxval_K, |
21 |
& GM_taper_scheme, |
22 |
& GM_Scrit, |
23 |
& GM_Sd |
24 |
|
25 |
|
26 |
C More GM/Redi parameters diagnosed from previous block |
27 |
C (not directly user configured) |
28 |
_RL GM_rMaxSlope |
29 |
COMMON /GM_PARAMS2/ GM_rMaxSlope |
30 |
|
31 |
|
32 |
C GM/Redi tensor elements |
33 |
|
34 |
C Bottom row of tensor corresponds to W points |
35 |
C Kwx is K_31 element, X direction at W point |
36 |
C Kwy is K_32 element, Y direction at W point |
37 |
C Kwz is K_33 element, Z direction at W point |
38 |
_RL Kwx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,MAX_NO_THREADS) |
39 |
_RL Kwy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,MAX_NO_THREADS) |
40 |
_RL Kwz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,MAX_NO_THREADS) |
41 |
COMMON /GM_Wtensor/ Kwx,Kwy,Kwz |
42 |
|
43 |
#ifdef GM_NON_UNITY_DIAGONAL |
44 |
C First/second rows of tensor corresponds to U/V points |
45 |
C Kux is K_11 element, X direction at U point |
46 |
C Kvy is K_22 element, Y direction at V point |
47 |
_RL Kux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,MAX_NO_THREADS) |
48 |
_RL Kvy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,MAX_NO_THREADS) |
49 |
COMMON /GM_UVtensor/ Kux,Kvy |
50 |
#else |
51 |
_RL Kux,Kvy |
52 |
PARAMETER(Kux=1.,Kvy=1.) |
53 |
#endif |
54 |
|
55 |
#ifdef GM_VISBECK_VARIABLE_K |
56 |
C GM mixing/stirring coefficient (spatially variable in horizontal |
57 |
C for Visbeck et al. parameterization) |
58 |
_RL VisbeckK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,MAX_NO_THREADS) |
59 |
COMMON /GM_Visbeck/ VisbeckK |
60 |
#endif |
61 |
|
62 |
#endif /* ALLOW_GMREDI */ |