Dear Gromac users, I need to count the number of H bonds as a function of time, but also as a function of the position of the water molecule in the box. By position I mean dividing the system in slices, as g_density does. Is there any way to tell g_hbonds to do this? Thank you so much, Marcelo Carignano