Global energy confinement scaling for neutral-beam-heated tokamaks