I was getting negative uncertainties that were messing up the specombine step. These lines get the error spectrum from frame 2
866: specerr = hdu[0].data[2, ap]
1017: specerrs[i] = interp(lamgrid, lam, hdu[0].data[2][ap] ** 2.0) ** 0.5
Looking at the spectra though, I think the uncertainty is in frame 3?