*INTERFACE:*

subroutine init_turbulence(namlst,fn,nlev)

Initialises all turbulence related stuff. This routine reads a number
of namelists and allocates memory for turbulence related vectors.
The core consists of calls to the the internal functions
`generate_model()` and `analyse_model()`, discussed in
great detail in section 4.7.3 and section 4.7.4, respectively.
The former function computes the model coefficients for the generic two-equation
model (see Umlauf et al. (2003)) from physically motivated quantities
like the von Kármán constant, , the decay rate in homogeneous
turbulence, , the steady-state Richardson number, ,
and many others. The latter function does the inverse:
it computes the physically motivated quantities from the model constants
of any model currently available in GOTM.
After the call to either function, all relevant model parameters
are known to GOTM.
Then, the function `report_model()` is called, which displays all
results on the screen.

*USES:*

IMPLICIT NONE

integer, intent(in) :: namlst character(len=*), intent(in) :: fn integer, intent(in) :: nlev

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

Karsten Bolding 2012-01-24