I mean the number of contacts per protein residue with different different parts of the lipids. It more or less can be done in GROMACS, but you need to create many indexes for each trajectory, so it is quite a long analysis. I was wondering if there is any tcl script to use with VMD to do so.