How to obtain the average force per contact by diameter

Hi all,

I want to analyze the average force per contact by diameter.

Using the command "compute pair/gran/local force id pos", I can obtain the id of the two contacted particles. Then using the command "compute property/atom diameter", I can obtain the diameter of each id. Based on some coding in Matlab, I can obtain the average force per contact by diameter. This procedure may not be that direct and I want to know is there any better method to output what I need?

As I see, we may cannot output the diameter of the two contacted particles using command "compute pair/gran/local" directly. Do I have to output two files and then do some post-coding? Any suggestions are greatly appreciated.

Best,
Rachel