Index: /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/ImageOps.c
===================================================================
--- /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/ImageOps.c	(revision 36799)
+++ /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/ImageOps.c	(revision 36800)
@@ -859,5 +859,5 @@
     N = 0;
     for (k = 0; k < catalog[0].average[j].Nmeasure; k++) {
-      //reset flag on each invocation
+      // reset flag on each invocation
       measure[k].dbFlags &= ~ID_MEAS_POOR_ASTROM;
 
@@ -898,4 +898,7 @@
   double *R, *D, *dR, *dD, *d2;
   StatType statsR, statsD;
+
+  // XXX we are not going to use this for now
+  return;
 
   Nsecfilt = GetPhotcodeNsecfilt();
Index: /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/UpdateObjects.c
===================================================================
--- /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/UpdateObjects.c	(revision 36799)
+++ /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/UpdateObjects.c	(revision 36800)
@@ -140,4 +140,5 @@
 	// does the measurement pass the supplied filtering constraints?
 	// MeasFilterTestTiny does not test psfQF
+	// exclude bad detections based on: photcodes, psfQF, time range, photflags & astromBadMask, mag_inst
 	int keepMeasure = measureBig ? MeasFilterTest(&measureBig[k], FALSE) : MeasFilterTestTiny(&measure[k], FALSE);
 	if (!keepMeasure) {
@@ -147,6 +148,26 @@
 	}
 
+	// mark (as POOR) any measurements which are deviant from the mean by > ExcludeBogusRadius
+	if (ExcludeBogus) {
+	  double Ri = getMeanR (&measure[k], &catalog[i].average[j], &catalog[i].secfilt[j*Nsecfilt]);
+	  double Di = getMeanD (&measure[k], &catalog[i].average[j], &catalog[i].secfilt[j*Nsecfilt]);
+	  coords.crval1 = catalog[i].average[j].R;
+	  coords.crval2 = catalog[i].average[j].D;
+	  double Xi, Yi;
+	  RD_to_XY (&Xi, &Yi, Ri, Di, &coords);
+	  double radius = hypot(Xi, Yi);
+	  if (radius > ExcludeBogusRadius) {
+	      measure[k].dbFlags &= ~ID_MEAS_USED_OBJ;
+	      if (measureBig) { measureBig[k].dbFlags &= ~ID_MEAS_USED_OBJ; }
+	      measure[k].dbFlags |= ID_MEAS_POOR_ASTROM;
+	      if (measureBig) { measureBig[k].dbFlags |= ID_MEAS_POOR_ASTROM; }
+	      continue;
+	  }
+	  measure[k].dbFlags &= ~ID_MEAS_POOR_ASTROM;
+	  if (measureBig) { measureBig[k].dbFlags &= ~ID_MEAS_POOR_ASTROM; }
+	}
+
 	// outlier rejection
-	if (FlagOutlier && (measure[k].dbFlags & ID_MEAS_POOR_ASTROM)) {
+	if (FALSE && FlagOutlier && (measure[k].dbFlags & ID_MEAS_POOR_ASTROM)) {
 	  measure[k].dbFlags &= ~ID_MEAS_USED_OBJ;
 	  if (measureBig) { measureBig[k].dbFlags &= ~ID_MEAS_USED_OBJ; }
@@ -271,6 +292,4 @@
 
       if (mode == FIT_PM_AND_PAR) {
-	// fprintf (stderr, "parallax fitting is still untested (%s, %d)\n", __FILE__, __LINE__);
-
 	float pXmin = +2.0;
 	float pXmax = -2.0;
Index: /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/args.c
===================================================================
--- /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/args.c	(revision 36799)
+++ /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/args.c	(revision 36800)
@@ -179,4 +179,5 @@
   FlagOutlier = FALSE;
   if ((N = get_argument (argc, argv, "-clip"))) {
+    fprintf (stderr, "-clip is currently disabled\n");
     remove_argument (N, &argc, argv);
     CLIP_THRESH = atof (argv[N]);
@@ -598,4 +599,5 @@
   FlagOutlier = FALSE;
   if ((N = get_argument (argc, argv, "-clip"))) {
+    fprintf (stderr, "-clip is currently disabled\n");
     remove_argument (N, &argc, argv);
     CLIP_THRESH = atof (argv[N]);
Index: /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/bcatalog.c
===================================================================
--- /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/bcatalog.c	(revision 36799)
+++ /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/bcatalog.c	(revision 36800)
@@ -104,5 +104,6 @@
 
       // filter out outliers - these are detections inconsistent with the offset distribution
-      if (FlagOutlier && (catalog[0].measure[offset].dbFlags & ID_MEAS_POOR_ASTROM)) {
+      // XXX disable this for now
+      if (FALSE && FlagOutlier && (catalog[0].measure[offset].dbFlags & ID_MEAS_POOR_ASTROM)) {
 	catalog[0].measure[offset].dbFlags &= ~ID_MEAS_USED_CHIP;
 	if (FALSE && (abs(catalog[0].measure[offset].t - Tref) < 10) && (catalog[0].measure[offset].photcode == Cref)) {
@@ -122,5 +123,4 @@
 	  RD_to_XY (&Xi, &Yi, Ri, Di, &coords);
 	  double radius = hypot(Xi, Yi);
-	  // XXX this radius needs to be configurable
 	  if (radius > ExcludeBogusRadius) {
 	      NskipBogus ++;
Index: /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/load_catalogs.c
===================================================================
--- /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/load_catalogs.c	(revision 36799)
+++ /branches/eam_branches/ipp-20140423/Ohana/src/relastro/src/load_catalogs.c	(revision 36800)
@@ -6,5 +6,6 @@
 Catalog *load_catalogs (SkyList *skylist, int *Ncatalog, int subselect, int hostID, char *hostpath) {
 
-  int i, j, k, m, Nstar;
+  int i, j, Nstar;
+  // int k, m;
   Catalog *catalog, *pcatalog, tcatalog;
 
@@ -73,4 +74,5 @@
       if (RESET) {
 	for (j = 0; j < catalog[Ncat].Naverage; j++) {
+# if (0)	  
 	  catalog[Ncat].average[j].flags = 0;
 	  m = catalog[Ncat].average[j].measureOffset;
@@ -78,4 +80,5 @@
 	    catalog[Ncat].measure[m+k].dbFlags = 0;
 	  }
+# endif
 	}
       }