I wanted to count the number of water molecules around the protein for a very huge microsecond long trajectory. Is there any specific way to do that with the help of gromacs tools? Using VMD with tickle its and easy task but the size of the trajectory make this a difficult one.