1 | ! |
---|
2 | ! $Id: mixing.h,v 1.6 2005/10/10 13:40:18 pmarches Exp $ |
---|
3 | ! |
---|
4 | ! This is include file "mixing.h" |
---|
5 | ! ==== == ======= ==== ========== |
---|
6 | ! |
---|
7 | #ifdef UV_VIS2 |
---|
8 | real visc2_r(GLOBAL_2D_ARRAY) |
---|
9 | real visc2_p(GLOBAL_2D_ARRAY) |
---|
10 | common /mixing_visc2_r/visc2_r /mixing_visc2_p/visc2_p |
---|
11 | #endif |
---|
12 | #ifdef UV_VIS4 |
---|
13 | real visc4_r(GLOBAL_2D_ARRAY) |
---|
14 | real visc4_p(GLOBAL_2D_ARRAY) |
---|
15 | common /mixing_visc4_r/visc4_r /mixing_visc4_p/visc4_p |
---|
16 | #endif |
---|
17 | #ifdef TS_DIF2 |
---|
18 | real diff2(GLOBAL_2D_ARRAY,NT) |
---|
19 | common /mixing_diff2/diff2 |
---|
20 | #endif |
---|
21 | #ifdef TS_DIF4 |
---|
22 | real diff4(GLOBAL_2D_ARRAY,NT) |
---|
23 | common /mixing_diff4/diff4 |
---|
24 | #endif |
---|
25 | #if defined UV_VIS2 || defined TS_DIF2 |
---|
26 | # ifdef SMAGORINSKY |
---|
27 | real visc3d_r(GLOBAL_2D_ARRAY,N) |
---|
28 | real visc3d_p(GLOBAL_2D_ARRAY,N) |
---|
29 | common /mixing_visc3d_r/visc3d_r /mixing_visc3d_p/visc3d_p |
---|
30 | # endif |
---|
31 | #endif |
---|
32 | |
---|
33 | #ifdef SOLVE3D |
---|
34 | real Akv(GLOBAL_2D_ARRAY,0:N) |
---|
35 | real Akt(GLOBAL_2D_ARRAY,0:N,2) |
---|
36 | common /mixing_Akv/Akv /mixing_Akt/Akt |
---|
37 | #ifdef RANDOM_WALK |
---|
38 | real dAktdz(GLOBAL_2D_ARRAY,0:N) |
---|
39 | common /mixing_dAktdz/dAktdz |
---|
40 | #endif |
---|
41 | |
---|
42 | # if defined ANA_VMIX || defined BVF_MIXING \ |
---|
43 | || defined LMD_MIXING || defined LMD_SKPP || defined LMD_BKPP \ |
---|
44 | || defined MY2_MIXING || defined MY25_MIXING |
---|
45 | real bvf(GLOBAL_2D_ARRAY,0:N) |
---|
46 | common /mixing_bvf/ bvf |
---|
47 | # endif |
---|
48 | |
---|
49 | # ifdef MY25_MIXING |
---|
50 | ! |
---|
51 | ! Mellor-Yamada (1982) Level 2.5 vertical mixing variables. |
---|
52 | ! Akq Vertical mixing coefficient (m^2/s) for turbulent energy. |
---|
53 | ! Lscale Turbulent length scale (m). |
---|
54 | ! q2 Turbulent energy squared (m^2/s^2) at horizontal |
---|
55 | ! RHO-points and vertical W-points. |
---|
56 | ! q2l Turbulent energy squared times turbulent length scale |
---|
57 | ! [m^3/s^2] at horizontal RHO-points and vertical |
---|
58 | ! W-points. |
---|
59 | ! |
---|
60 | real Akq(GLOBAL_2D_ARRAY,0:N) |
---|
61 | real Lscale(GLOBAL_2D_ARRAY,N) |
---|
62 | real q2(GLOBAL_2D_ARRAY,0:N,2) |
---|
63 | real q2l(GLOBAL_2D_ARRAY,0:N,2) |
---|
64 | common /my25_mix_Akq/Akq /my25_mix_Lscale/Lscale |
---|
65 | & /my25_mix_q2/q2 /my25_mix_q2l/q2l |
---|
66 | #endif /* MY25_MIXING */ |
---|
67 | |
---|
68 | # if defined LMD_SKPP || defined LMD_BKPP |
---|
69 | ! |
---|
70 | ! Large/McWilliams/Doney oceanic planetary boundary layer variables. |
---|
71 | ! ghats Boundary layer nonlocal transport (m/s^2). |
---|
72 | ! hbl Depth of oceanic surface boundary layer (m). |
---|
73 | ! hbbl Depth of oceanic bottom boundary layer (m). |
---|
74 | ! kbl Index of first grid level below "hbl". |
---|
75 | ! ustar Turbulent friction velocity (m/s). |
---|
76 | ! |
---|
77 | integer kbl(GLOBAL_2D_ARRAY) |
---|
78 | # ifdef BIOLOGY |
---|
79 | real hel(GLOBAL_2D_ARRAY) |
---|
80 | common /lmd_hel/hel |
---|
81 | # endif |
---|
82 | real hbl(GLOBAL_2D_ARRAY) |
---|
83 | c vech 30/08/07 |
---|
84 | real hbl_0(GLOBAL_2D_ARRAY) |
---|
85 | real hbbl(GLOBAL_2D_ARRAY) |
---|
86 | common /lmd_kpp_kbl/kbl /lmd_kpp_hbl/hbl |
---|
87 | & /lmd_kpp_hbbl/hbbl |
---|
88 | # ifdef LMD_NONLOCAL |
---|
89 | real ghats(GLOBAL_2D_ARRAY,0:N) |
---|
90 | common /lmd_kpp_ghats/ghats |
---|
91 | # endif |
---|
92 | # endif /* LMD_SKPP || LMD_BKPP */ |
---|
93 | |
---|
94 | # ifdef LMD_MIXING |
---|
95 | real ustar(GLOBAL_2D_ARRAY) |
---|
96 | common /lmd_kpp_ustar/ustar |
---|
97 | # endif /* LMD_MIXING */ |
---|
98 | |
---|
99 | #endif /* SOLVE3D */ |
---|
100 | |
---|