wki function

Winkler index