Changeset 35066

Timestamp:

Jan 28, 2013, 9:27:49 AM (13 years ago)

Author:

mhuber

Message:

bug fix so max-density not pull all/most sources from catalog

File:

• trunk/Ohana/src/getstar/src/select_by_region.c (modified) (3 diffs)

trunk/Ohana/src/getstar/src/select_by_region.c

@@ -53 +53 @@
     gfits_scan (&catalog[0].header, "DEC1", "%lf", 1, &Dmax);
 
-    if (VERBOSE) fprintf (stderr, "extracting from catalog covering region %f,%f to %f,%f\n", Rmin, Dmin, Rmax, Dmax);
-
     AREA = fabs(Dmax - Dmin) * fabs(Rmax - Rmin) * cos (0.5*RAD_DEG*(Dmax + Dmin));
+    if (VERBOSE) fprintf (stderr, "extracting from catalog covering region %f,%f to %f,%f for area %f\n", Rmin, Dmin, Rmax, Dmax, AREA);
     assert (AREA > 0);
 
@@ -83 +82 @@
       if (isnan(mag)) continue;
       if (mag < MinMagValue) continue;
-      if (mag > MagLimitValue) continue;
 
       bin = (mag - MagMin) / dMag;
@@ -101 +99 @@
     }
     MagLimitValue = i*dMag + MagMin;
-    if (VERBOSE) fprintf (stderr, "using %d stars in mag range %f - %f\n", Nsum, MinMagValue, MagLimitValue);
+    if (VERBOSE) fprintf (stderr, "using %d (Nmax %d) stars in mag range %f - %f\n", Nsum, Nmax, MinMagValue, MagLimitValue);
   }