Use separate variables for eddy diffusion coefficient and drag coefficient in procedure coefkz (following LMDZ). coefkzmin only computes eddy diffusion coefficient, not drag coefficient.
trunk/Sources/phylmd/clmain.f