In the area of statistical mechanics, an important open question is to show that the height function associated with the square ice model (i.e., planar six vertex model with uniform weights), or equivalently the uniform graph homeomorphisms, converges to a continuum Gaussian free field In the scaling limit, I will review some recent results about this model, including that the single point height function, upon renormalization, converges to a Gaussian random variable.