Did you plot columns 2 and 3? Did any of them look like a g(r)?
Column 2 is g(r) and 3 is basically its integral, from which you can get the "solvation number". Coluimn 1 is not the number of frames, but it is quite easy to figure out what it is.
Since you are not doing a solvent-solvent g(r), but rather a site (your reference point) to solvent g(r) it may be that you should specify iSdist instead of iGdist - I don't remember exactly how the code works for this combination.
