We present a global analysis of available data on inclusive structure functions measured in electron-proton scattering at small values of Bjorken-x, including the latest data from the combined HERA analysis on reduced cross sections. Our approach relies on the dipole formulation of DIS together with the use of the non-linear running coupling BK equation for the description of the small-x dynamics. With the resulting parametrization we are able to describe the latest FL data measured by the H1 collaboration. Further, we discuss the kinematical domain where significant deviations from NLO-DGLAP should be expected and the ability of non-linnear physics to account for such deviations.