We present a recently introduced strategy to study non-perturbatively thermal QCD up to temperatures of the order of the electro-weak scale, combining step scaling techniques and shifted boundary conditions. The former allow to renormalize the theory for a range of scales which spans several orders of magnitude with a moderate computational cost. Shifted boundary conditions avoid the need for the zero temperature subtraction in the Equation of State. As a consequence, the simulated lattices do not have to accommodate two very different scales, the pion mass and the temperature. Effective field theory arguments guarantee that finite volume effects can be kept under control safely. As a first application of this strategy, we present the results of the computation of the hadronic screening spectrum in QCD with $N_f=3$ flavours of massless quarks for temperatures from $T\sim 1$ GeV up to $\sim 160$ GeV.