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.