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.