code: def separation_cos_angle(lon0, lat0, lon1, lat1):
      
      return (np.sin(lat1) * np.sin(lat0) + np.cos(lat1) * np.cos(lat0) *
              np.cos(lon1 - lon0))
docs: Evaluate the cosine of the angular separation between two
      direction vectors.