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.
IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: namlst character(len=*), intent(in) :: fn integer, intent(in) :: nlevREVISION HISTORY:
Original author(s): Karsten Bolding, Hans Burchard, Manuel Ruiz Villarreal, Lars Umlauf