Normalized integrated spectra of the NaD absorption line profile after the subtraction of the stellar contribution in galaxies with medium and high quality stellar modeling of the integrated spectra (Sect. 3.1). Each spectrum shown covers a rest-frame range of ~46 Å (i.e., ~2300 km s^-1). For each source, in the upper panels, the different kinematic components, are shown in blue and green, i.e., blueshifted and systemic or redshifted component, according to our classification (Table 1). The red curve shows the total contribution coming from the NaD Gaussian fit. In pink, we indicate the model of the HeI emission line. This line was successfully modeled with one (e.g., F10257-4339) or two components (e.g., F04315-0840) in 21 cases. In two cases, F07027-6011 (N) and F10015-0614, the line modeling of the (prominent) HeI emission cannot reproduce the observed line profile well. Therefore, in these cases, the HeI is simply excluded when fitting the NaD absorption. In the lower panels, the residuals (i.e., data – model) are presented with the rest-frame NaD wavelengths indicated in orange. The galaxy IDs follow that of Table 1, although here they are shortened for better visualization.