###

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
- 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

Karsten Bolding
2012-12-28