source: trunk/LMDZ.MARS/libf/phymars/turb_mod.F90 @ 1635

Last change on this file since 1635 was 1242, checked in by aslmd, 11 years ago

LMDZ.MARS. MESOSCALE. Corrected problems affecting LES runs with new physics (sensibFlux calculated in vdifc is just used as an input to thermals, it does not have the right sign and apparently is not W/m2). Further reduction or factorisation of mesoscale includes. Also made ustar and tstar available in turb_mod from the calculations in vdifc.

File size: 994 bytes
Line 
1module turb_mod
2
3  !! variables
4  REAL,SAVE,ALLOCATABLE :: q2(:,:)    ! Turbulent Kinetic Energy
5  REAL,allocatable,SAVE :: l0(:)
6  REAL,SAVE,ALLOCATABLE :: ustar(:)
7  REAL,SAVE,ALLOCATABLE :: wstar(:)
8  REAL,SAVE,ALLOCATABLE :: tstar(:)
9  REAL,SAVE,ALLOCATABLE :: hfmax_th(:)
10  REAL,SAVE,ALLOCATABLE :: zmax_th(:)
11  REAL,SAVE,ALLOCATABLE :: sensibFlux(:)
12  LOGICAL :: turb_resolved = .false.
13      ! this is a flag to say 'turbulence is resolved'
14      ! mostly for LES use. default is FALSE (for GCM and mesoscale)
15
16contains
17
18  subroutine ini_turb_mod(ngrid,nlayer)
19 
20  implicit none
21  integer,intent(in) :: ngrid ! number of atmospheric columns
22  integer,intent(in) :: nlayer ! number of atmospheric layers
23
24    allocate(q2(ngrid,nlayer+1))
25    allocate(l0(ngrid))
26    allocate(wstar(ngrid))
27    allocate(ustar(ngrid))
28    allocate(tstar(ngrid))
29    allocate(hfmax_th(ngrid))
30    allocate(zmax_th(ngrid))
31    allocate(sensibFlux(ngrid))
32
33  end subroutine ini_turb_mod
34
35end module turb_mod
Note: See TracBrowser for help on using the repository browser.