Changeset 37655 for branches/eam_branches/ipp-20140904/Ohana/src/relastro/src/ICRF.c

Timestamp:

Nov 20, 2014, 8:34:19 PM (12 years ago)

Author:

eugene

Message:

fix parallel relastro load catalog, move load icrf out of bcatalog

File:

• branches/eam_branches/ipp-20140904/Ohana/src/relastro/src/ICRF.c (modified) (1 diff)

branches/eam_branches/ipp-20140904/Ohana/src/relastro/src/ICRF.c

@@ -48 +48 @@
   return TRUE;
 }
+
+int select_catalog_ICRF (Catalog *catalog, int Ncatalog) {
+
+  if (!USE_ICRF_CORRECT) return TRUE;
+
+  int N = 0;
+
+  char filename[1024];
+  snprintf (filename, 1024, "%s/test.icrf.dat", CATDIR);
+  FILE *f = fopen (filename, "w");
+
+  int c, i, j;
+  for (c = 0; c < Ncatalog; c++) {
+
+    for (i = 0; i < catalog[c].Naverage; i++) {
+      if (!(catalog[c].average[i].flags & ID_ICRF_QSO)) continue;
+
+      // only save a single value
+      int savedICRF = FALSE;
+      for (j = 0; j < catalog[c].average[i].Nmeasure; j++) {
+          
+        int offset = catalog[c].average[i].measureOffset + j;
+          
+        if (!(catalog[c].measureT[offset].dbFlags & ID_MEAS_ICRF_QSO)) continue;
+
+        fprintf (f, "%d %d %d : %f %f : %f %f\n", c, i, offset,
+                 catalog[c].average[i].R, catalog[c].average[i].D,
+                 catalog[c].measureT[offset].R, catalog[c].measureT[offset].D);
+
+        if (savedICRF) continue;
+
+        ICRFsave (c, i, offset);
+        N ++;
+        savedICRF = TRUE;
+      }
+    }
+  }
+  fclose (f);
+
+  fprintf (stderr, "added %d ICRF QSO\n", N);
+
+  return TRUE;
+}