| [2089] | 1 | |
|---|
| 2 | module Mod_PHY_AT_ctr |
|---|
| 3 | |
|---|
| 4 | |
|---|
| 5 | !--------------------------------------------------------------------------+ |
|---|
| 6 | ! Sat 25-May-2013 MAR | |
|---|
| 7 | ! module Mod_PHY_AT_ctr contains the logical variables used to set up | |
|---|
| 8 | ! MAR PHYsics | |
|---|
| 9 | ! | |
|---|
| 10 | ! version 3.p.4.1 created by H. Gallee, Mon 1-Apr-2013 | |
|---|
| 11 | ! as Mod_PHY____ctr | |
|---|
| 12 | ! Last Modification by H. Gallee, Sat 25-May-2013 | |
|---|
| 13 | ! | |
|---|
| 14 | !--------------------------------------------------------------------------+ |
|---|
| 15 | |
|---|
| 16 | |
|---|
| 17 | use Mod_Real |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | IMPLICIT NONE |
|---|
| 21 | |
|---|
| 22 | |
|---|
| 23 | ! E-e and K-l models switches |
|---|
| 24 | ! --------------------------- |
|---|
| 25 | |
|---|
| 26 | logical :: AT_TKE ! Set-Up of Turbulent Kinetic Energy Computation |
|---|
| 27 | |
|---|
| 28 | logical :: Ee_Duynkerke ! Dunkerke (1988) E-epsilon model of turbulence |
|---|
| 29 | logical :: Ee_Kitada ! Kitada (1987) E-epsilon model of turbulence |
|---|
| 30 | logical :: Ee_HuangRamn ! Huang and Raman (1991) E-epsilon model of turbulence |
|---|
| 31 | logical :: Kl_TherryLac ! Therry & Lacarrere (1983) K-l model of turbulence |
|---|
| 32 | |
|---|
| 33 | |
|---|
| 34 | ! Box weighted vertical moving Averages of TKE and eps |
|---|
| 35 | ! ---------------------------------------------------- |
|---|
| 36 | |
|---|
| 37 | logical :: AT_vav = .FALSE. ! Box weighted vertical moving Averages of TKE & eps |
|---|
| 38 | integer, SAVE :: mz__KA = 1 ! Level No below which moving Averages of TKE and eps is effective |
|---|
| 39 | real(kind=real8), SAVE :: zz__KA = 20. ! Height below which moving Averages of TKE and eps is effective |
|---|
| 40 | |
|---|
| 41 | |
|---|
| 42 | ! Flux Limitor |
|---|
| 43 | ! ------------ |
|---|
| 44 | |
|---|
| 45 | logical :: AT_LIM = .FALSE. ! pkt Turbulent Flux Limitor |
|---|
| 46 | real(kind=real8), SAVE :: DThMAX = 5. ! Maximum allowed d(PT ) in 1 d(t) [K/h] |
|---|
| 47 | real(kind=real8), SAVE :: Dpkt_X ! Maximum allowed d(pkt) in 1 d(t) = DThMAX /100.**(R/Cp) /3600. [K/s] |
|---|
| 48 | |
|---|
| 49 | |
|---|
| 50 | end module Mod_PHY_AT_ctr |
|---|
