Index: /trunk/psphot/src/psphotChooseAnalysisOptions.c
===================================================================
--- /trunk/psphot/src/psphotChooseAnalysisOptions.c	(revision 36452)
+++ /trunk/psphot/src/psphotChooseAnalysisOptions.c	(revision 36453)
@@ -250,7 +250,7 @@
 	    psSphere ptGal, ptSky;
 	    GetGalacticCoords (&ptGal, &ptSky, toGal, chip, source->peak->xf, source->peak->yf);
-            float b = ptGal.r;
-            float limit = GAL_LIMIT + GAL_LIMIT_BULGE * exp(-0.5*(b*b/GAL_LIMIT_SIGMA2));
-            if (fabs(ptGal.d) < limit) continue;
+            float l = ptGal.r;
+            float b_min = GAL_LIMIT + GAL_LIMIT_BULGE * exp(-0.5*(l*l/GAL_LIMIT_SIGMA2));
+            if (fabs(ptGal.d) < b_min) continue;
 	    // include an exception for low density skycells below the limit?
 	}
@@ -533,7 +533,7 @@
 		psSphere ptGal, ptSky;
 		GetGalacticCoords (&ptGal, &ptSky, toGal, chips->data[imageID], source->peak->xf, source->peak->yf);
-                float b = ptGal.r;
-                float limit = GAL_LIMIT + GAL_LIMIT_BULGE * exp(-0.5*(b*b/GAL_LIMIT_SIGMA2));
-		if (fabs(ptGal.d) < limit) continue;
+                float l = ptGal.r;
+                float b_min = GAL_LIMIT + GAL_LIMIT_BULGE * exp(-0.5*(l*l/GAL_LIMIT_SIGMA2));
+		if (fabs(ptGal.d) < b_min) continue;
 	    }