  It has been well known for a long time that the height function of random
lozenge tilings of large domains follow a law of large number and possible
limits called dimer limit shapes are well understood. For the next order, it is
expected that fluctuations behave like version of a Gaussian Free field, at
least away from some special "frozen" regions. However despite being one of the
main questions in the domain for 20 years, only special cases have been
obtained. In this paper we show that for any specified limit shape with no
frozen region, one can construct a sequence of domains whose height functions
converge to that limit shape and where the height fluctuation converge to a
variant of the Gaussian Free Field.
