haversine function

Haversine great circle distance