In this article, we prove that the height function associated with the square-ice model (i.e.~the six-vertex model with a=b=c=1 on the square lattice), or, equivalently, of the uniform random homomorphisms from \mathbb Z^2 to \mathbb Z, has logarithmic variance. This establishes a strong form of roughness of this height function.