Using all available data on the deep-inelastic cross-sections at HERA at
x<0.01, we look for geometric scaling of the form \sigma^{\gamma^*p}(\tau)
where the scaling variable \tau behaves alternatively like \log(Q^2)-\lambda Y,
as in the original definition, or \log(Q^2)-\lambda \sqrt{Y}, which is
suggested by the asymptotic properties of the Balitsky-Kovchegov (BK) equation
with running QCD coupling constant. A ``Quality Factor'' (QF) is defined,
quantifying the phenomenological validity of the scaling and the uncertainty on
the intercept \lambda. Both choices have a good QF, showing that the second
choice is as valid as the first one, predicted for fixed coupling constant. A
comparison between the QCD asymptotic predictions and data is made and the QF
analysis shows that the agreement can be reached, provided going beyond leading
logarithmic accuracy for the BK equation.Comment: 4 pages, 4 figure