A six-dimensional interaction potential for the water dimer has been fitted to ab initio interaction energies computed at 2510 dimer configurations. These energies were obtained by combining the supermolecular second-order energies extrapolated to the complete basis set limit from up to quadruple-zeta quality basis sets with the contribution from the coupled-cluster method including single, double, and noniterative triple excitations computed in a triple-zeta quality basis set. All basis sets were augmented by diffuse functions and supplemented by midbond functions. The energies have been fitted using an analytic form with the induction component represented by a polarizable term, making the potential directly transferable to clusters and the bulk phase. Geometries and energies of stationary points on the potential surface agree well with the results of high-level ab initio geometry optimizations.