A multimode-like scheme for selecting the centers of Gaussian basis functions when computing vibrational spectra