We have performed studies of the three-dimensional random-field XY model on 32 samples of L×L×L simple cubic lattices with periodic boundary conditions, with a random field strength of h_{r} = 1.5, for L= 128, using a parallelized Monte Carlo algorithm. We present results for the sample-averaged magnetic structure factor S(k[over ⃗]) over a range of temperature, using both random hot start and ferromagnetic cold start initial states, and k[over ⃗] along the [1,0,0] and [1,1,1] directions. At T= 1.875, S(k[over ⃗]) shows a broad peak near |k[over ⃗]|=0, with a correlation length which is limited by thermal fluctuations, rather than the lattice size. As T is lowered, this peak grows and sharpens. By T= 1.5, it is clear that the correlation length is larger than L= 128. The lowest temperature for which S(k[over ⃗]) was calculated is T= 1.421875, where the hot start and cold start initial conditions usually do not find the same local minimum in the phase space. Our results are consistent with the idea that there is a finite value of T below which S(k[over ⃗]) diverges slowly as |k[over ⃗]| goes to zero. This divergence would imply that the relaxation time of the spins is also diverging. That is the signature of an ergodicity-breaking phase transition.