Changeset 8649

Timestamp:

Aug 27, 2006, 1:38:02 PM (20 years ago)

Author:

eugene

Message:

temporary test for inconsistent average/measure coords

File:

• trunk/Ohana/src/opihi/dvo/pmeasure.c (modified) (2 diffs)

trunk/Ohana/src/opihi/dvo/pmeasure.c

@@ -7 +7 @@
   
   FILE *f;
-  int i, j, k, m, N, VERBOSE;
+  int i, j, k, m, N, status, VERBOSE;
   double Mz, Mr, mag;
-  double Radius, Rmin, Rmax;
+  double Radius, Rmin, Rmax, R, D;
   unsigned IDclip, IDchoice, LimExclude;
   unsigned FlagChoice, FlagClip, flags;
@@ -122 +122 @@
         if (LimExclude && (Zvec.elements[N] > 0.99)) continue;
         if (Zvec.elements[N] < 0.011) continue;
-        fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], catalog.average[i].R - catalog.measure[m+k].dR_PS/3600.0, catalog.average[i].D - catalog.measure[m+k].dD_PS/3600.0, &graphmode.coords);
+        R = catalog.average[i].R - catalog.measure[m+k].dR_PS/3600.0;
+        D = catalog.average[i].D - catalog.measure[m+k].dD_PS/3600.0;
+        // XXX drop this check
+        if ((R < Rmin) || (R > Rmax) || (D < -90.0) || (D > 90.0)) {
+          char *date;
+          date = sec_to_date (catalog.measure[m+k].t);
+          gprint (GP_LOG, "out: %f, %f : %s : (%f, %f) + (%f, %f)\n", R, D, date, catalog.average[i].R, catalog.average[i].D, catalog.measure[m+k].dR_PS/3600.0, catalog.measure[m+k].dD_PS/3600.0);
+          free (date);
+        }
+        status = fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], R, D, &graphmode.coords);
+        if (!status) {
+          fprintf (stderr, ".");
+          continue;
+        }
         N ++;
       }