The algebraic velocity variances (Source File: variances.F90)

INTERFACE:

    subroutine variances(nlev,SSU,SSV)

DESCRIPTION:

Using (70) and the solution shown in (74) and the variances of the turbulent velocity fluctations can be evaluated according to

$$\begin{array}{rcl} \dfrac{\langle u'^2 \rangle}{k} &=& \dfrac... ...{3}-a_3\right) P +\dfrac83 a_5 G \right) \; , \quad \end{array}$$ (180)

where the diffusivities are computed according to (44) (also see section 4.7.39 and section 4.7.40), and the buoyancy production, $G$, follows from (147).

USES:

   use turbulence,  only:     uu,vv,ww
   use turbulence,  only:     tke,eps,P,B,num
   use turbulence,  only:     cc1,ct1,a2,a3,a5
   IMPLICIT NONE

INPUT PARAMETERS:

 
   number of vertical layers
   integer,  intent(in)                 :: nlev
 
   square of shear frequency (1/s^2)
   (from u- and v-component)
   REALTYPE, intent(in)                :: SSU(0:nlev),SSV(0:nlev)

REVISION HISTORY:

    Original author(s): Lars Umlauf
    $Log: variances.F90,v $
    Revision 1.1  2005-06-27 10:54:33  kbk
    new files needed