See roman_imsim.utils.roman_utils.getPSF_Image: https://github.com/DukeCosmology/roman_imsim/blob/main/roman_imsim/utils.py#L148 We are missing this in our implementation of the galsim PSFs.