Left panel: rotation-rate profile of 16 Cyg A corresponding to the PM estimates of Ω₀ and Ω₁. The dashed line marks the bottom of the convective envelope in the assumed stellar model. Right panel: distribution of the surface rotation rate as a function of the stellar co-latitude. The red shade represents probability densities at each latitude. The black line marks the surface rotation model for the PM values of Ω₀ and Ω₁ and corresponds to the map in the right panel at r = R_⋆. The black dashed lines mark the mean surface rotation model plus (right curve) or minus (left curve) the corresponding posterior standard deviation estimated at each latitude.