We examine the behavior of a function sampled from the invariant measure associated to the focusing discrete Non Linear Schrödinger equation, defined on a discrete torus of dimension $d \geq 3$, and nonlinearity parameter $p>4$, in the infinite volume limit. The Gibbs measure has two parameters, the inverse temperature and the strength of the non linearity, and the scaling is such that the non linear and linear parts of the Hamiltonian contribute on the same scale. It was shown by Dey, Kirkpatrick and the first author that this measure undergoes a phase transition, concerning concentration of mass of a typical function {at the level of partition functions}. We prove that the three regions of the phase diagram yield three distinct local limits, a massive Gaussian free field, a massless Gaussian free field plus a random constant, and finally a (possibly trivial) mixture of massive Gaussian free fields, where some mass is ``lost'' to the region of concentration. Our proof relies on the analysis of the spherical model of a ferromagnet. The measure under consideration is an exponential tilt of the spherical model, and the additional tilt can be seen to induce an additional phase transition. Although our motivations come from the NLS equation and soliton resolution conjecture, our proofs are completely probabilistic, and can be read without any knowledge of PDE theory.