Adding function pmSubtractionMatch to do the PSF-matching. It calls all the other functions, doing what ppSub has done for a long time. However, since we want to use the same functionality in ppStack, I moved it in here to avoid code repitition. Currently untested.