Employing a viscosity-based survival scale for initial geometrical perturbations formed in relativistic heavy-ion collisions, we model the radial flow velocity at freeze-out. Subsequently, we use the Cooper-Frye freeze-out prescription, with viscous corrections to the distribution function, to extract the transverse momentum dependence of particle yields and flow harmonics. We fit the model parameters for central collisions, by fitting the spectra of identified particles at the Large Hadron Collider (LHC), and estimate them for other centralities using simple hydrodynamic relations. We use the results of Monte Carlo Glauber model for initial eccentricities. We demonstrate that this improved viscous blast-wave model leads to good agreement with transverse momentum distribution of elliptic and triangular flow for all centralities and estimate the shear viscosity to entropy density ratio η/s ≃ 0.24 at the LHC.