The local, weak-equilibrium stability functions (Source File: cmue_c.F90)

INTERFACE:

    subroutine cmue_c(nlev)

DESCRIPTION:

This subroutine updates the explicit solution of (70) and (71) with shape indicated by (74). In addition to the simplifications discussed in section 4.7.38, $P_b=\epsilon_b$ is assumed in (71) to eliminate the dependency on $\overline{T}$ according to (72). As discussed in section 4.4, this implies that the last of (71) is replaced by (73). Thus, the $\Gamma$-term in (74) drops out, and the solution is characterized by $c_\mu$ and $c'_\mu$ only.

As a consequence, the numerators and the denominator appearing in (77) are of somewhat different form compared to the result in section 4.7.37. They can be written as

$$\begin{array}{rcl} D &=& d_0 + d_1 \overline{N}^2 + d_2 \over... ...{b1} \overline{N}^2 + n_{b2} \overline{S}^2 \quad . \end{array}$$ (191)

The coefficients of $D$ are given by

$$\begin{array}{rcl} d_0 &=& 36 {\cal N}^3 {\cal N}_b^2 \; , \q... ...2 - a_2^2)(a_{b2}^2 - a_{b1}^2) {\cal N} \; , \quad \end{array}$$ (192)

and the coefficients of the numerators are

$$\begin{array}{rcl} n_0 &=& 36 a_1 {\cal N}^2 {\cal N}_b^2 \; ... ...[3mm] n_2 &=& 9 a_1 (a_{b2}^2 - a_{b1}^2){\cal N}^2 \end{array}$$ (193)

and

$$\begin{array}{rcl} n_{b0} &=& 12 a_{b3} {\cal N}^3 {\cal N}_b... ... 4 (a_2^2-3 a_3^2) ) {\cal N} {\cal N}_b \; , \quad \end{array}$$ (194)

These polynomials correspond to a slightly generalized form of the solution suggested by Canuto et al. (2001) and Cheng et al. (2002). For cases with unstable stratification, the same clipping conditions on $\alpha_N$ is applied as described in section 4.7.40. For the cases of extreme shear, the limiter described in the context of (83) is active.

USES:

    use turbulence, only: an,as,at
    use turbulence, only: cmue1,cmue2
    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

DEFINED PARAMETERS:

    REALTYPE, parameter       :: asLimitFact=1.0d0
    REALTYPE, parameter       :: anLimitFact=0.5d0

REVISION HISTORY:

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