ini_cturbs.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
SUBROUTINE ini_cturbs
!!
!!****     *INI_CTURBS*  - routine to initialize the turbulence scheme 
!!                        constants.
!!
!!      PURPOSE
!!      -------
!         The purpose of this routine is to initialize the turbulence 
!       scheme constants that are stored in module MODD_CTURB
!
!!      METHOD
!!      ------
!!        The constants are set to their numerical values
!!
!!      AUTHOR
!!      ------
!!        V. Masson         * Meteo-France *
!!
!!      MODIFICATIONS
!!      -------------
!! --------------------------------------------------------------------------
!
!*        0. DECLARATIONS
!            ------------
!
USE modd_csts,        ONLY : xkarman
USE modd_canopy_turb, ONLY : xtkemin, xcep, xced, xalpsbl, xa0, xcmfs, xasbl,&
                               xcshf, xctp  
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!  ---------------------------------------------------------------------------
!
!         1. SETTING THE NUMERICAL VALUES
!            ----------------------------
!
!         1.1 Constant for dissipation of Tke
!

REAL(KIND=JPRB) :: zhook_handle

IF (lhook) CALL dr_hook('INI_CTURBS',0,zhook_handle)
xced  = 0.845
!       Redelsperger-Sommeria (1981) = 0.70
!       Schmidt-Schumann      (1989) = 0.845
!       Cheng-Canuto-Howard   (2002) = 0.845
!
!
!         1.2 Constant for wind pressure-correlations
!
xcep  = 2.11
!       Redelsperger-Sommeria (1981) = 4.
!       Schmidt-Schumann      (1989) = 3.5
!       Cheng-Canuto-Howard   (2002) = 2.11
!
!
!         1.3 Constant a0 for wind pressure-correlations
!
xa0   = 0.6
!       Redelsperger-Sommeria (1981) = 0.6
!       Schmidt-Schumann      (1989) = 0.55
!       Cheng-Canuto-Howard   (2002) = 0.6
!
!
!         1.4 Constant for temperature and vapor pressure-correlations
!
xctp  = 4.65
!       Redelsperger-Sommeria (1981) = 4.
!       Schmidt-Schumann      (1989) = 3.25
!       Cheng-Canuto-Howard   (2002) = 4.65
!
!
!         1.5 Value related to the TKE universal function within SBL
!
xalpsbl = 4.63
!       Redelsperger et al 2001     = 4.63
!       Wyngaard et al. 1974        = 3.75
!       Stull 1988                  = 4.75
!
!
!
!         2. Derivated constants
!            -------------------
!
!         2.1 Constant in fluxes equations
!
xcmfs= 2./3./xcep*(1.-xa0)   !Constant for the momentum flux due to shear (RS)
!
! Redelsperger-Sommeria (1981) ......... 0.066
! Schmidt-Schumann      (1989) ......... 0.086
! Cheng-Canuto-Howard   (2002) ......... 0.126
!
xcshf= 2./3./xctp            !Constant for the sensible heat flux(RS)
!
! Redelsperger-Sommeria (1981) ......... 0.167
! Schmidt-Schumann      (1989) ......... 0.204
! Cheng-Canuto-Howard   (2002) ......... 0.143
!
!
!         2.2 Value related to the TKE universal function within SBL
!
!
xasbl   = 0.5*( xalpsbl**(3./2.)*xkarman*xced + xkarman/sqrt(xalpsbl)/xcmfs )
!       Redelsperger et al 2001
!
!
!
!         3. MINIMUM VALUES 
!            --------------
!
xtkemin=1.e-4
IF (lhook) CALL dr_hook('INI_CTURBS',1,zhook_handle)
!

END SUBROUTINE ini_cturbs