Changeset 24805

Timestamp:

Jul 15, 2009, 2:34:48 PM (17 years ago)

Author:

eugene

Message:

test code to track down refstar precision errors

File:

• trunk/psastro/src/psastroLoadRefstars.c (modified) (2 diffs)

trunk/psastro/src/psastroLoadRefstars.c

@@ -128 +128 @@
     if (!strcmp (outformat, "PS1_DEV_0")) {
       refstars = psastroReadGetstar_PS1_DEV_0 (fits);
+
+      // XXX test
+      // FILE *outfile = fopen ("refstars.dat", "w");
+      // assert (outfile);
+      // for (int nn = 0; nn < refstars->n; nn++) {
+      //          pmAstromObj *ref = refstars->data[nn];
+      //          fprintf (outfile, "%lf %lf\n", ref->sky->r*PS_DEG_RAD, ref->sky->d*PS_DEG_RAD);
+      // }
+      // fclose (outfile);
     }
     if (refstars == NULL) {
@@ -218 +227 @@
         psMetadata *row = psFitsReadTableRow(fits, i); // Table row
 
-        ref->sky->r   = RAD_DEG*psMetadataLookupF32 (&status, row, "RA");
-        ref->sky->d   = RAD_DEG*psMetadataLookupF32 (&status, row, "DEC");
+        ref->sky->r   = RAD_DEG*psMetadataLookupF64 (&status, row, "RA");
+        ref->sky->d   = RAD_DEG*psMetadataLookupF64 (&status, row, "DEC");
         ref->Mag      = psMetadataLookupF32 (&status, row, "MAG");
         float MagC1   = psMetadataLookupF32 (&status, row, "MAG_C1");