I wrote this function to minimize the distances between orders assigned by the model by penalizing any deviations from the centroid of each route.