A multicopy protocol is proposed for modeling macromolecular hydration using diffraction experimental data (X-ray or neutron) to search for a better description of delocalized water sites than that given by point water models. The model consists of one macro-molecule and several copies of each water molecule, refined simultaneously against diffraction data using molecular dynamics techniques. The protocol was applied to BPTI and an RNA tetradecamer. The sites defined by the different copies range from very ordered ones to continuous channels; they fit the density maps and agree with the diffraction amplitudes with an accuracy comparable with usual crystallographic methods. The delocalization of water in channels agrees with the high mobility observed in NMR experiments.