We study the effects of including a running coupling constant in high-density QCD evolution. For fixed coupling constant, QCD evolution preserves the initial dependence of the saturation momentum Q s on the nuclear size A and results in an exponential dependence onFor the running coupling case, we re-derive analytical estimates for the A-and Y -dependences of the saturation scale and test them nu-where we find numerically ∆ ′ ≃ 3.2. We study the behaviour of the gluon distribution at large transverse momentum, characterizing it by an anomalous dimension 1 − γ which we define in a fixed region of small dipole sizes. In contrast to previous analytical work, we find a marked difference between the fixed coupling (γ ≃ 0.65) and running coupling (γ ∼ 0.85) results. Our numerical findings show that both a scaling function depending only on the variable r Q s and the perturbative double-leadinglogarithmic expression, provide equally good descriptions of the numerical solutions for very small r-values below the so-called scaling window.