I like the use of argsort. However, a transpose won't be needed unless ```coord.ndim > 1```, so you could avoid the call to argsort and transpose for 1d coords as well as scalars. I appreciate that neither of these calls will be computationally expensive.