The averaged error of numerical solutions on stress intensity factors for semi-elliptical surface cracks has earlier been shown to depend on the gradient of stresses applied to crack faces. In particular, the error, Epsilon, increases appreciabely with increasing the load gradient in the length direction, i.e. with concentrating KI-distribution near the surface point of a crack front. Instead of the description of the error value via exponents in a polynomial stress distribution, the variation of Epsilon is studied in terms of the dimensionless parameter, Lambda, characterizing the level of nonuniformity of the KI-distribution along a crack front. Such a transformation of data reveals the monotonic behavior of Epsilon versus Lambda, with a rather clearly defined upper bound curve. The latter, being extrapolated to the limit values of Lambda=0 and Lambda=1, yields error estimates for local stress intensity factors at the surface and deepest points of a crack front respectively