Non gaussianity parameter in gromacs

Dear all,

Is there a module in Gromacs to get non-Gaussian parameters to study dynamical heterogeneity ? If not? how can I calculate it? is there any script for it?

Thank you