Index: /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/ReadImageHeader.c
===================================================================
--- /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/ReadImageHeader.c	(revision 37004)
+++ /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/ReadImageHeader.c	(revision 37005)
@@ -45,5 +45,6 @@
   image[0].NY = Ny;
 
-  image[0].refColor = NAN;
+  image[0].refColorBlue = NAN;
+  image[0].refColorRed = NAN;
 
   if (!gfits_scan (header, "TZERO",   "%d",  1, &image[0].tzero) && !ACCEPT_TIME) {
@@ -135,5 +136,6 @@
 
   tmp = 0;
-  gfits_scan (header, "FLIMIT",   "%lf", 1, &tmp);
+  // gfits_scan (header, "FLIMIT",   "%lf", 1, &tmp);
+  gfits_scan (header, "DETEFF.MAGREF", "%lf", 1, &tmp);
   image[0].detection_limit = tmp * 10.0;
 
Index: /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches.c
===================================================================
--- /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches.c	(revision 37004)
+++ /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches.c	(revision 37005)
@@ -36,5 +36,5 @@
   ALLOCATE (Y2, double, NAVE);
   ALLOCATE (N2, off_t,  NAVE);
-  ALLOCATE (catalog[0].found, off_t, NAVE);
+  ALLOCATE (catalog[0].found_t, off_t, NAVE);
   /* for secfilt j and star i, secfilt[i*Nsecfilt+j] */
 
@@ -77,5 +77,4 @@
   if (Nstars < 1) {
     if (VERBOSE) fprintf (stderr, "skipping %s, no overlapping stars\n", catalog[0].filename);
-    free (catalog[0].found);
     free (X1);
     free (Y1);
@@ -92,5 +91,5 @@
     RD_to_XY (&X2[i], &Y2[i], catalog[0].average[i].R, catalog[0].average[i].D, &tcoords);
     N2[i] = i;
-    catalog[0].found[N2[i]] = -1;
+    catalog[0].found_t[N2[i]] = -1;
   }
   if (Nave > 1) sort_coords_index (X2, Y2, N2, Nave);
@@ -219,9 +218,9 @@
       }
       /* this catalog star matches more than one image star */
-      if (catalog[0].found[n] > -1) {
-	catalog[0].measure[catalog[0].found[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
+      if (catalog[0].found_t[n] > -1) {
+	catalog[0].measure[catalog[0].found_t[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
 	catalog[0].measure[Nmeas].dbFlags |= ID_MEAS_BLEND_OBJ;
       } else {
-	catalog[0].found[n] = Nmeas;
+	catalog[0].found_t[n] = Nmeas;
       }
       /* Nm is updated, but not written out in -update mode (for existing entries)
@@ -329,5 +328,4 @@
   if (VERBOSE) fprintf (stderr, "Nstars, Nave, Nmeas: "OFF_T_FMT" "OFF_T_FMT" "OFF_T_FMT", ("OFF_T_FMT" matches)\n",  Nstars,  Nave,  Nmeas,  Nmatch);
 
-  free (catalog[0].found);
   free (X1);
   free (Y1);
Index: /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest.c
===================================================================
--- /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest.c	(revision 37004)
+++ /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest.c	(revision 37005)
@@ -36,5 +36,5 @@
   ALLOCATE (Y2, double, NAVE);
   ALLOCATE (N2, off_t,  NAVE);
-  ALLOCATE (catalog[0].found, off_t, NAVE);
+  ALLOCATE (catalog[0].found_t, off_t, NAVE);
   /* for secfilt j and star i, secfilt[i*Nsecfilt+j] */
 
@@ -77,5 +77,4 @@
   if (Nstars < 1) {
     if (VERBOSE) fprintf (stderr, "skipping %s, no overlapping stars\n", catalog[0].filename);
-    free (catalog[0].found);
     free (X1);
     free (Y1);
@@ -92,5 +91,5 @@
     RD_to_XY (&X2[i], &Y2[i], catalog[0].average[i].R, catalog[0].average[i].D, &tcoords);
     N2[i] = i;
-    catalog[0].found[N2[i]] = -1;
+    catalog[0].found_t[N2[i]] = -1;
   }
   if (Nave > 1) sort_coords_index (X2, Y2, N2, Nave);
@@ -261,5 +260,5 @@
        Nm is recalculated in build_meas_links if loaded table is not sorted */
     stars[N].found = Nmeas;
-    catalog[0].found[n] = Nmeas;
+    catalog[0].found_t[n] = Nmeas;
     catalog[0].average[n].Nmeasure ++;
     Nmeas ++;
@@ -390,5 +389,4 @@
   if (VERBOSE) fprintf (stderr, "Nstars, Nave, Nmeas, Nlens: "OFF_T_FMT" "OFF_T_FMT" "OFF_T_FMT" "OFF_T_FMT", ("OFF_T_FMT" matches)\n",  Nstars,  Nave,  Nmeas,  Nlens, Nmatch);
 
-  free (catalog[0].found);
   free (X1);
   free (Y1);
Index: /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest_refstars.c
===================================================================
--- /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest_refstars.c	(revision 37004)
+++ /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_closest_refstars.c	(revision 37005)
@@ -31,5 +31,5 @@
   ALLOCATE (Y2, double, NAVE);
   ALLOCATE (N2, off_t,  NAVE);
-  ALLOCATE (catalog[0].found, off_t, NAVE);
+  ALLOCATE (catalog[0].found_t, off_t, NAVE);
   /* for secfilt j and star i, secfilt[i*Nsecfilt+j] */
 
@@ -71,5 +71,4 @@
   if (Nstars < 1) {
     if (VERBOSE) fprintf (stderr, "skipping %s, no overlapping stars\n", catalog[0].filename);
-    free (catalog[0].found);
     free (X1);
     free (Y1);
@@ -86,5 +85,5 @@
     RD_to_XY (&X2[i], &Y2[i], catalog[0].average[i].R, catalog[0].average[i].D, &tcoords);
     N2[i] = i;
-    catalog[0].found[N2[i]] = -1;
+    catalog[0].found_t[N2[i]] = -1;
   }
   if (Nave > 1) sort_coords_index (X2, Y2, N2, Nave);
@@ -220,9 +219,9 @@
     /* this catalog star matches more than one image star */
     // XXX should this be an average flag?
-    if (catalog[0].found[n] > -1) {
-      catalog[0].measure[catalog[0].found[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
+    if (catalog[0].found_t[n] > -1) {
+      catalog[0].measure[catalog[0].found_t[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
       catalog[0].measure[Nmeas].dbFlags |= ID_MEAS_BLEND_OBJ;
     } else {
-      catalog[0].found[n] = Nmeas;
+      catalog[0].found_t[n] = Nmeas;
     }
 
@@ -331,5 +330,4 @@
   if (VERBOSE) fprintf (stderr, "Nstars, Nave, Nmeas: "OFF_T_FMT" "OFF_T_FMT" "OFF_T_FMT", ("OFF_T_FMT" matches)\n",  Nstars,  Nave,  Nmeas,  Nmatch);
 
-  free (catalog[0].found);
   free (X1);
   free (Y1);
Index: /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_refstars.c
===================================================================
--- /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_refstars.c	(revision 37004)
+++ /branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/find_matches_refstars.c	(revision 37005)
@@ -37,5 +37,5 @@
   ALLOCATE (Y2, double, NAVE);
   ALLOCATE (N2, off_t,  NAVE);
-  ALLOCATE (catalog[0].found, off_t, NAVE);
+  ALLOCATE (catalog[0].found_t, off_t, NAVE);
   REALLOCATE (catalog[0].average, Average, NAVE);
   REALLOCATE (catalog[0].secfilt, SecFilt, NAVE*catalog[0].Nsecfilt);
@@ -69,5 +69,5 @@
     RD_to_XY (&X2[i], &Y2[i], catalog[0].average[i].R, catalog[0].average[i].D, &tcoords);
     N2[i] = i;
-    catalog[0].found[N2[i]] = -1;
+    catalog[0].found_t[N2[i]] = -1;
   }
   if (Nave > 1) sort_coords_index (X2, Y2, N2, Nave);
@@ -190,9 +190,9 @@
       }
       /* this catalog star matches more than one image star */
-      if (catalog[0].found[n] > -1) {
-	catalog[0].measure[catalog[0].found[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
+      if (catalog[0].found_t[n] > -1) {
+	catalog[0].measure[catalog[0].found_t[n]].dbFlags |= ID_MEAS_BLEND_OBJ;
 	catalog[0].measure[Nmeas].dbFlags |= ID_MEAS_BLEND_OBJ;
       } else {
-	catalog[0].found[n] = Nmeas;
+	catalog[0].found_t[n] = Nmeas;
       }
 
@@ -305,5 +305,4 @@
   if (VERBOSE) fprintf (stderr, "Nstars, Nave, Nmeas: %d "OFF_T_FMT" "OFF_T_FMT", ("OFF_T_FMT" matches)\n",  Nstars,  Nave,  Nmeas,  Nmatch);
 
-  free (catalog[0].found);
   free (X1);
   free (Y1);