Correct length scale calculation and other bugs in kEpsilonABL

[DriesAllaerts]

1. lm is calculated with k and epsilon instead of using an analytical function
2. denominator of lmax should be max() instead of min()
3. Richardson number should be -B/P such that Ri<0 corresponds to unstable conditions