Calculating Non-polar contacts

GROMACS version:
GROMACS modification: Yes/No
Hello! I would like to calculate hydrophobic contacts within a specific domain of my protein. To do that, I created an index file where the hydrophobic residues of this domain are selected. Which command should I use to calculate the possible contacts between these residues?