  We evaluate the unintegrated gluon distribution of the proton starting from a
parametrization of the color dipole cross section including DGLAP evolution and
saturation effects. To this end, we perform the Fourier-Bessel transform of
$\sigma(x,r)/\alpha(r)$. At large transverse momentum of gluons we match the
so-obtained distribution to the logarithmic derivative of the collinear gluon
distribution. We check our approach by calculating the proton structure
function $F_L(x,Q^2)$ finding good agreement with HERA data.
