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.
IMPLICIT NONEINPUT 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