Index: trunk/Ohana/src/relastro/src/BrightCatalog.c
===================================================================
--- trunk/Ohana/src/relastro/src/BrightCatalog.c	(revision 39994)
+++ trunk/Ohana/src/relastro/src/BrightCatalog.c	(revision 39995)
@@ -294,5 +294,5 @@
       secfilt[i].Mmax   = Mmax[i];
     }
-    fprintf (stderr, "loaded data for %lld averages\n", (long long) Nrow);
+    // fprintf (stderr, "loaded data for %lld secfilt\n", (long long) Nrow);
 
     free (M    );
Index: trunk/Ohana/src/relastro/src/bcatalog.c
===================================================================
--- trunk/Ohana/src/relastro/src/bcatalog.c	(revision 39994)
+++ trunk/Ohana/src/relastro/src/bcatalog.c	(revision 39995)
@@ -365,5 +365,5 @@
   gfits_scan (&oldcatalog[0].header, "DEC1", "%lf", 1, &Dmax);
 
-  if (VERBOSE2) fprintf (stderr, "extracting from catalog covering region %f,%f to %f,%f\n", Rmin, Dmin, Rmax, Dmax);
+  if (VERBOSE) fprintf (stderr, "extracting from catalog covering region %f,%f to %f,%f (%d average)\n", Rmin, Dmin, Rmax, Dmax, (int) subcatalog[0].Naverage);
 
   float AREA = fabs(Dmax - Dmin) * fabs(Rmax - Rmin) * cos (0.5*RAD_DEG*(Dmax + Dmin));
@@ -459,6 +459,6 @@
     if (isnan(X[i])) continue;
     if (isnan(Y[i])) continue;
-    ix = MIN(0, MAX(NX - 1, (X[i] - Xmin) / 200.0));
-    iy = MIN(0, MAX(NY - 1, (Y[i] - Ymin) / 200.0));
+    ix = MAX(0, MIN(NX - 1, (X[i] - Xmin) / 200.0));
+    iy = MAX(0, MIN(NY - 1, (Y[i] - Ymin) / 200.0));
     int Nn = Nn_grid[ix][iy];
     Nm_grid[ix][iy][Nn] = subcatalog[0].average[i].Nmeasure;