Update the turbulent kinetic energy

INTERFACE:

   subroutine do_tke(nlev,dt,u_taus,u_taub,z0s,z0b,h,NN,SS)

DESCRIPTION:

Based on user input, this routine calls the appropriate routines for calculating the turbulent kinetic energy. The user has the choice between an algebraic equation described in section 4.17, and two versions of the dynamic transport equation of the TKE described in section 4.11 and section 4.12. The former uses $k$-$\epsilon$ notation, the latter the notation of Mellor and Yamada (1982). Apart from this, both equations are identical and update the vectors tke and tkeo, which is the value of the tke at the old time step.

USES:

   IMPLICIT NONE

INPUT PARAMETERS:

   integer,  intent(in)                :: nlev
   REALTYPE, intent(in)                :: dt,u_taus,u_taub,z0s,z0b
   REALTYPE, intent(in)                :: h(0:nlev)
   REALTYPE, intent(in)                :: NN(0:nlev),SS(0:nlev)

REVISION HISTORY:

   Original author(s): Karsten Bolding, Hans Burchard,
                       Manuel Ruiz Villarreal, Lars Umlauf