Index: /trunk/psphot/src/psphot.h
===================================================================
--- /trunk/psphot/src/psphot.h	(revision 9806)
+++ /trunk/psphot/src/psphot.h	(revision 9807)
@@ -51,5 +51,5 @@
 int             psphotSortBySN (const void **a, const void **b);
 int             psphotSortByY (const void **a, const void **b);
-bool            psphotGrowthCurve (pmReadout *readout, pmPSF *psf);
+bool            psphotGrowthCurve (pmReadout *readout, pmPSF *psf, bool ignore);
 void            psphotTestArguments (int *argc, char **argv);
 bool            psphotMaskReadout (pmReadout *readout, psMetadata *recipe);
Index: /trunk/psphot/src/psphotApResid.c
===================================================================
--- /trunk/psphot/src/psphotApResid.c	(revision 9806)
+++ /trunk/psphot/src/psphotApResid.c	(revision 9807)
@@ -25,4 +25,5 @@
     // S/N limit to perform full non-linear fits
     float MAX_AP_OFFSET = psMetadataLookupF32 (&status, recipe, "MAX_AP_OFFSET");
+    bool  IGNORE_GROWTH = psMetadataLookupBool (&status, recipe, "IGNORE_GROWTH");
 
     // set limits on the aperture magnitudes
@@ -32,5 +33,6 @@
     float REF_RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_REF_RADIUS");
     psf->growth = pmGrowthCurveAlloc (3.0, 100.0, REF_RADIUS);
-    psphotGrowthCurve (readout, psf);
+
+    psphotGrowthCurve (readout, psf, IGNORE_GROWTH);
 
     psVector *mask    = psVectorAllocEmpty (300, PS_TYPE_U8);
@@ -67,5 +69,5 @@
 
         flux->data.F64[Npsf]    = pow(10.0, -0.4*source->psfMag);
-        r2rflux->data.F64[Npsf] = PS_SQR(model->radiusTMP) / flux->data.F64[Npsf];
+        r2rflux->data.F64[Npsf] = PS_SQR(model->radiusFit) / flux->data.F64[Npsf];
 
         mask->data.U8[Npsf] = 0;
Index: /trunk/psphot/src/psphotGrowthCurve.c
===================================================================
--- /trunk/psphot/src/psphotGrowthCurve.c	(revision 9806)
+++ /trunk/psphot/src/psphotGrowthCurve.c	(revision 9807)
@@ -5,5 +5,7 @@
 //     the 'center' option
 
-bool psphotGrowthCurve (pmReadout *readout, pmPSF *psf) {
+// XXX add a option to turn off the curve-of-growth (ie, make the apMag = fitMag everywhere);
+
+bool psphotGrowthCurve (pmReadout *readout, pmPSF *psf, bool ignore) {
 
     // bool status;
@@ -31,4 +33,8 @@
     pmModel *model = pmModelFromPSF (modelRef, psf);
 
+    // measure the fitMag for this model
+    pmSourcePhotometryModel (&fitMag, model);
+    psf->growth->fitMag = fitMag;
+
     // generate working image for this source
     psRegion region = {xc - dx, xc + dx, yc - dy, yc + dy};
@@ -46,15 +52,15 @@
 
         psImageKeepCircle (mask, xc, yc, radius, "OR", PM_MASK_MARK);
+        pmModelAdd (image, mask, model, false, false);
+        pmSourcePhotometryAper (&apMag, model, image, mask);
+        psImageKeepCircle (mask, xc, yc, radius, "AND", PS_NOT_U8(PM_MASK_MARK));
 
-        pmModelAdd (image, mask, model, false, false);
-
-        pmSourcePhotometryAper (&apMag, model, image, mask);
-
-        psImageKeepCircle (mask, xc, yc, radius, "AND", PS_NOT_U8(PM_MASK_MARK));
-        psf->growth->apMag->data.F32[i] = apMag;
+	// the 'ignore' mode is for testing
+	if (ignore) {
+	    psf->growth->apMag->data.F32[i] = fitMag;
+	} else {
+	    psf->growth->apMag->data.F32[i] = apMag;
+	}
     }
-
-    pmSourcePhotometryModel (&fitMag, model);
-    psf->growth->fitMag = fitMag;
 
     psf->growth->apRef = psVectorInterpolate (psf->growth->radius, psf->growth->apMag, psf->growth->refRadius);