The non-local, approximate weak-equilibrium stability function


INTERFACE:

   subroutine cmue_b(nlev)
DESCRIPTION:

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

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


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



Karsten Bolding 2012-12-28