Fancy indexing — choose several positions at once. In a spreadsheet: “give me rows 3, 8 and 20”.
top5 = np.argsort(radius)[-5:] profiles = X[top5, :] # (5, 30) — in ONE operation
Passing an array of positions pulls all of them out together, in the order you asked for.
