In this paper we review and update the Amaral-Gay Ducati-Betemps-Soyez saturation model, by testing it against the recent H1-ZEUS combined data on deep inelastic scattering, including heavy quarks in the dipole amplitude. We obtain that this model, which is based on traveling wave solutions of the Balitsky-Kovchegov equation and built in the momentum space framework, yields very accurate descriptions of the reduced cross section, $��_{r}(x,y,Q^{2})$, as well as DIS structure functions such as $F_{2}(x,Q^{2})$ and $F_{L}(x,Q^{2})$, all measured at HERA. Additionally, it provides good descriptions of heavy quark structure functions, $F_{2}^{cc}$ and $F_{2}^{bb}$ at small-$x$ and $Q^{2}\lesssim 60$ GeV$^{2}$. We also use the improved model to make predictions for structure functions to be measured in the near future at LHeC.

13 pages, 7 figures