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

$\displaystyle {\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



Karsten Bolding 2012-01-24