Random dimer models (or equivalently tiling models) have been a subject of extensive research in mathematics and physics for several decades. In this talk, we will discuss the doubly periodic Aztec diamond dimer model of growing size, with arbitrary periodicity and only mild conditions on the edge weights. In this limit, we see three types of macroscopic regions -- known as rough, smooth and frozen regions. We will discuss how the geometry of the arctic curves, the boundary of these regions, can be described in terms of an associated amoeba and an action function. In particular, we determine the number of frozen and smooth regions and the number of cusps on the arctic curves. We will also discuss the convergence of local fluctuations to the appropriate translation-invariant Gibbs measures. Joint work with Alexei Borodin.