  We show that the numerical solution of the classical SU(3) Yang-Mills
equations of motion in the McLerran-Venugopalan model for gluon production in
central heavy ion collisions leads to a suppresion at low $p_t$ and an
enhancement at the intermediate $p_t$ region as compared to peripheral heavy
ion and pp collisions at the same energy. Our results are compared to previous,
Color Glass Condensate inspired calculations of gluon production in heavy ion
collisions. We revisit the predictions of the Color Glass Condensate model for
$pA$ ($dA$) collisions in Leading Order and show that quantum evolution--in
particular the phenomenon of geometric scaling and change of anomalous
dimensions--preserves the Cronin enhancement of $pA$ cross section (when
normalized to the leading twist term) in the Leading Order approximation even
though the $p_t$ spectrum can change. We comment on the case when gluon
radiation is included.
