We present the extraction of unpolarized quark transverse momentum dependent parton distribution functions (TMDPDFs) and the non-perturbative part of TMD evolution kernel from the global analysis of Drell-Yan and $Z$-boson production data. The analysis is performed at the next-to-next-to-leading order (NNLO) in perturbative QCD, using the $zeta$-prescription. The estimation of the error-propagation from the experimental uncertainties to non-perturbative function is made by the replica method. The importance of the inclusion of the precise LHC data and its influence on the determination of non-perturbative functions is discussed.