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.