How can I calculate the charge distribution of a water molecule?
Two suggestions: MolCalc: Calculate Properties > Polarity and Solvation gives you a 3D surface plot with the electrostatic potential superimposed.

Another suggestion is Avogadro/GAMESS/MacMolPlt. See for example here