Building upon recent results of Dub��dat (see arXiv:1403.6076) on the convergence of topological correlators in the double-dimer model considered on Temperleyan approximations $��^��$ to a simply connected domain $��\subset\mathbb C$ we prove the convergence of probabilities of cylindrical events for the \emph{double-dimer loop ensembles} on $��^��$ as $��\to 0$. More precisely, let $��_1,\dots,��_n\in��$ and $L$ be a macroscopic lamination on $��\setminus\{��_1,\dots,��_n\}$, i.e., a collection of disjoint simple loops surrounding at least two punctures considered up to homotopies. We show that the probabilities $P_L^��$ that one obtains $L$ after withdrawing all loops surrounding no more than one puncture from a double-dimer loop ensemble on $��^��$ converge to a conformally invariant limit $P_L$ as $��\to 0$, for each $L$. Though our primary motivation comes from 2D statistical mechanics and probability, the proofs are of a purely analytic nature. The key techniques are the analysis of entire functions on the representation variety $\mathrm{Hom}(��_1(��\setminus\{��_1,\dots,��_n\})\to\mathrm{SL}_2(\mathbb C))$ and on its (non-smooth) subvariety of locally unipotent representations. In particular, we do \emph{not} use any RSW-type arguments for double-dimers. The limits $P_L$ of the probabilities $P_L^��$ are defined as coefficients of the isomonodormic tau-function studied by Dub��dat with respect to the Fock--Goncharov lamination basis on the representation variety. The fact that $P_L$ coincides with the probability to obtain $L$ from a sample of the nested CLE(4) in $��$ requires a small additional input, namely a mild crossing estimate for this nested conformal loop ensemble., minor update following referee's comments (42 pages, 5 figures)