  We consider a version of the McLerran-Venugopalan model by Lam and Mahlon
where confinement is implemented via colored noise in the infrared. This model
does not assume an infinite momentum frame, hence the boosted nuclei are not
infinitely thin. Instead, the nuclei have a finite extension in the
longitudinal direction and therefore depend on the longitudinal coordinate. In
this fully three dimensional framework an x dependence of the gluon
distribution function emerges naturally. In order to fix the parameters of the
model, we calculate the gluon distribution function and compare it with the
JR09 parametrization of the data. We explore the parameter space of the model
to attain a working framework that can be used to calculate the initial
conditions in heavy ion collisions.
