The local, weak-equilibrium stability functions

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.25, $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.24. 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.27. 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