Cleaned up some messy stuff that was only there for my debugging interests. Also, giving sampled in any array works fine. This comes at the price of looping over arguments, but as long as you are not simulating a lot (not sure exactly what number) of galaxies this is not a big deal. Fine for the number of galaxies in our size coadds, so it should be okay. The only reason I mention it is python looping is slow.