Pair correlation functions calculation

Dear ATAT team,

I have been doing SQS calculations and have a question that has been bothering me for a long time. Your answer would be greatly appreciated. I would like to know if there is a command in ATAT to compute the pairwise correlation function that reaches the second nearest neighbour for a specific user-defined alloy structure? Can you write out the command specifically?

Thanks a lot.

Kind regards,
Shuai Li