  We introduce a modified form of the Kharzeev-Levin-Nardi (KLN) approach for
nuclear collisions. The new ansatz for the unintegrated gluon distribution
function preserves factorization, and the saturation scale is bound from below
by that for a single nucleon. It also reproduces the correct scaling with the
number of collisions at high transverse momentum. The corresponding Monte Carlo
implementation allows us to account for fluctuations of the hard sources
(nucleons) in the transverse plane. We compute various definitions of the
eccentricity within the new approach, which are relevant for the interpretation
of the elliptic flow. Our approach predicts breaking of the scaling of the
eccentricity with the Glauber eccentricity at the level of about 30%.
