function sdssphot_get_sersicfit, objid, filter n = n_elements(objid) ; exists = bytarr(n) sersic_fits = replicate({ROMSersicFit}, n) sersic_fit = {ROMSersicFit} filters = ['i', 'r', 'g'] channel = where(filters eq filter, count) if count eq 0 then begin print, 'Cannot find filter' return, -1 endif for i = 0, n-1 do begin ; check to see if the file exists filepath = sdss_rom_sersic_path(objid, (*pclass).filters[channel], /full) if file_test(filepath) then begin ; exists[i] = 1B struct_assign, mrdfits(filepath, 1), sersic_fit sersic_fits[i] = sersic_fit endif else begin ; exists[i] = 0B endelse endfor return, sersic_fits end