SelectedPoints(ctype, psel) with the default multiplicity == 0 creates an empty, uninitialized field that keeps psel; pass a positive multiplicity to allocate now.
qlat.qlat_gpt