Index: trunk/psphot/src/fit_galaxies.c
===================================================================
--- trunk/psphot/src/fit_galaxies.c	(revision 4114)
+++ trunk/psphot/src/fit_galaxies.c	(revision 4251)
@@ -8,9 +8,4 @@
     int   Nfit = 0;
     int   Niter = 0;
-
-    float XBORDER  = psMetadataLookupF32 (&status, config, "XBORDER");
-    float YBORDER  = psMetadataLookupF32 (&status, config, "YBORDER");
-    psRegion *keep = psRegionAlloc (XBORDER, -XBORDER, YBORDER, -YBORDER);
-    keep           = psRegionForImage (keep, image, keep);
 
     float RADIUS   = psMetadataLookupF32 (&status, config, "FIT_RADIUS");
@@ -25,6 +20,11 @@
 
 	// need a better model guess and a better radius choice
-	pmSourceSetPixelCircle (source, image, RADIUS);
-	pmSourceMaskRegion (source, keep);
+	// when radius is not fixed, we will need to check if new radius fits on image
+# if (0)
+	if (model->radius > OUTER) {
+	  // allocate image, noise, mask arrays for each peak (square of radius OUTER)
+	  pmSourceDefinePixels (source, imdata, x, y, OUTER);
+	}
+# endif
 	psModel  *model  = pmSourceModelGuess (source, modelType);