Finding Atlas regions for top X% of sources

[Geeta]

I want to do the above operation using scripting.