  We present an extraction of unpolarized partonic transverse momentum
distributions (TMDs) from a simultaneous fit of available data measured in
semi-inclusive deep-inelastic scattering, Drell-Yan and Z-boson production. To
connect data at different scales, we use TMD evolution at next-to-leading
logarithmic accuracy. The analysis is restricted to the low-transverse-momentum
region, with no matching to fixed-order calculations at high transverse
momentum. We introduce specific choices to deal with TMD evolution at low
scales, of the order of 1 GeV$^2$. This could be considered as a first attempt
at a global fit of TMDs.
