[2] | 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) |
---|
[4] | 83 | c vech 30/08/07 |
---|
| 84 | real hbl_0(GLOBAL_2D_ARRAY) |
---|
[2] | 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 | |
---|