We report on results for the Landau gauge gluon propagator computed from large statistical ensembles and look at the compatibility of the results with the Gribov-Zwanziger tree level prediction for its refined and very refine