INTERFACE:
subroutine cmue_a(nlev)DESCRIPTION:
The solution of (72) and (73) has the shape indicated by (76). This subroutine is used to update the quantities , and , defined in (76), from which all turbulent fluxes can be computed. The non-linear terms and are updated by evaluating the right hand side of (67) at the old time step.
The numerators and the denominator appearing in (79) are polynomials of the form
The coefficients of are given by
USES:
use turbulence, only: eps use turbulence, only: P,B,Pb,epsb use turbulence, only: an,as,at,r 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 NONEINPUT PARAMETERS:
number of vertical layers integer, intent(in) :: nlevBUGS:
Test stage. Do not yet use.REVISION HISTORY:
Original author(s): Lars Umlauf
Karsten Bolding 2012-12-28