The non-local, approximate weak-equilibrium stability function (Source File: cmue_b.F90)

INTERFACE:

    subroutine cmue_b(nlev)

DESCRIPTION:

This subroutine is used to update the quantities $c_\mu$, $c'_\mu$ and $\Gamma$, defined in (74), from which all turbulent fluxes can be computed. This done exactly as described in section 4.7.37, with the exception that equilibrium $P+G=\epsilon$ and $P_b=\epsilon_b$ is assumed in computing the non-linear terms in (65), leading to the particularly simple expressions

$${\cal N} = \dfrac{c_1}{2} \; , \quad {\cal N}_b = c_{b1} \quad .$$ (190)

USES:

    use turbulence, only: an,as,at
    use turbulence, only: cmue1,cmue2,gam
    use turbulence, only: cm0
    use turbulence, only: cc1
    use turbulence, only: ct1,ctt
    use turbulence, only: a1,a2,a3,a4,a5
    use turbulence, only: at1,at2,at3,at4,at5
 
    IMPLICIT NONE

INPUT PARAMETERS:

    number of vertical layers
    integer, intent(in)       :: nlev

BUGS:

   Test stage. Do not yet use.

REVISION HISTORY:

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