  The JIMWLK equation with a "daughter dipole" running coupling is solved
numerically starting from an initial condition given by the
McLerran-Venugopalan model. The resulting Wilson line configurations are then
used to compute the spectrum of gluons comprising the glasma inital state of a
high energy heavy ion collision. The development of a geometrical scaling
region makes the spectrum of produced gluons harder. Thus the ratio of the mean
gluon transverse momentum to the saturation scale grows with energy. Also the
total gluon multiplicity increases with energy slightly faster than the
saturation scale squared.
