This paper is the first attempt to build a color glass condensate/saturation model based on the next-to-leading-order (NLO) corrections to linear and nonlinear evolution in QCD. We assume that the renormalization scale is the saturation momentum and find that the scattering amplitude has geometric scaling behavior deep in the saturation domain with the explicit formula of this behavior at large τ=r2Qs2. We build a model that includes this behavior, as well as the known ingredients: (i) the behavior of the scattering amplitude in the vicinity of the saturation momentum, using the NLO Balitsky-Fadin-Kuraev-Lipatov kernel, (ii) the pre-asymptotic behavior of ln(Qs2(Y)), as a function of Y, and (iii) the impact parameter behavior of the saturation momentum, which has exponential behavior exp(-mb) at large b. We demonstrate that the model is able to describe the experimental data for the deep inelastic structure function. Despite this, our model has difficulties that are related to the small value of the QCD coupling at Qs(Y0) and the large values of the saturation momentum, which indicate the theoretical inconsistency of our description.