Index: /branches/eam_branch_20080324/psphot/src/psphotSourceSize.c
===================================================================
--- /branches/eam_branch_20080324/psphot/src/psphotSourceSize.c	(revision 17117)
+++ /branches/eam_branch_20080324/psphot/src/psphotSourceSize.c	(revision 17118)
@@ -47,4 +47,33 @@
 	// XXX for now, just skip any sources with masked pixels
 	source->crNsigma  = -5.0;
+	bool keep = true;
+	for (int iy = -1; (iy <= +1) && keep; iy++) {
+	    for (int ix = -1; (ix <= +1) && keep; ix++) {
+		if (mask[yPeak+iy][xPeak+ix]) { keep &= false; }
+	    }
+	}
+	if (!keep) continue;
+
+	// need a reasonably fast way to follow the psf R = 1 annulus
+
+	// save the PSF model from the Ensemble fit
+	pmModel *PSF = source->modelPSF;
+
+	// convert model to polar terms
+	psEllipseShape shape;
+	psF32 *PAR = PSF->params->data.F32;
+	shape.sx   = PAR[PM_PAR_SXX] / M_SQRT2;
+	shape.sy   = PAR[PM_PAR_SYY] / M_SQRT2;
+	shape.sxy  = PAR[PM_PAR_SXY];
+
+	psEllipseShapeToAxes (shape, 20.0);
+	
+
+	for (theta = 0.0; theta < 360.0; theta += 10.0) {
+	  
+
+
+	// XXX for now, just skip any masked pixels
+	extSum = 0.0;
 	bool keep = true;
 	for (int iy = -1; (iy <= +1) && keep; iy++) {