Index: trunk/psphot/src/psphotReadout.c
===================================================================
--- trunk/psphot/src/psphotReadout.c	(revision 34317)
+++ trunk/psphot/src/psphotReadout.c	(revision 34404)
@@ -128,10 +128,4 @@
     }
 
-    // find blended neighbors of very saturated stars (detections->newSources)
-    // if (!psphotDeblendSatstars (config, view, filerule)) {
-    // psError (PSPHOT_ERR_UNKNOWN, false, "failed on satstar deblend analysis");
-    // return psphotReadoutCleanup (config, view, filerule);
-    // }
-
     // mark blended peaks PS_SOURCE_BLEND (detections->newSources)
     // XXX I've deactivated this because it was preventing galaxies close to stars from being
@@ -149,4 +143,10 @@
     }
 
+    // find and subtract radial profile models for saturated stars (XXX change name eventually)
+    if (!psphotDeblendSatstars (config, view, filerule)) {
+	psError (PSPHOT_ERR_UNKNOWN, false, "failed on satstar deblend analysis");
+	return psphotReadoutCleanup (config, view, filerule);
+    }
+
     // if we were not supplied a PSF model, determine the IQ stats here (detections->newSources)
     if (!psphotImageQuality (config, view, filerule)) { // pass 1
@@ -221,5 +221,5 @@
     // NOTE: this block performs the 2nd pass low-significance PSF detection stage
     { 
-	// add noise for subtracted objects
+	// add noise for subtracted objects & subtracted saturated stars
 	psphotAddNoise (config, view, filerule); // pass 1 (detections->allSources)
 
@@ -253,5 +253,4 @@
 	// merge the newly selected sources into the existing list
 	// NOTE: merge OLD and NEW
-	// XXX check on free of sources...
 	psphotMergeSources (config, view, filerule); // (detections->newSources + detections->allSources -> detections->allSources)
 
@@ -294,5 +293,4 @@
 	// merge the newly selected sources into the existing list
 	// NOTE: merge OLD and NEW
-	// XXX check on free of sources...
 	psphotMergeSources (config, view, filerule); // (detections->newSources + detections->allSources -> detections->allSources)