f-baerenbold
commented
3 years ago According to Rodi, 1987 JGR and Burchard et al., 1995 JGR, the value in Goudsmit et al., 2002 JGR is wrong. This means that the value of 0.112 in the code is presumably correct.
If stability mode is 1 (defined in config file), the turbulent diffusivity of temperature is calculated as: nu_t = cmue_2*k^2/eps, where cmue_2 = 0.072 according to Goudsmit et al., 2002.
However, in the code, cmue_2 = cmue/Prndtl = 0.09/0.8 = 0.112, so turbulent diffusivity of temperature might be overestimated by 50% (if stability mode = 1).