Use separate variables for eddy diffusion coefficient and drag coefficient in procedure coefkz (following LMDZ). coefkzmin only computes eddy diffusion coefficient, not drag coefficient.