Changeset 4844 for trunk/Ohana/src/opihi/dvo/cmatch.c

Timestamp:

Aug 20, 2005, 5:40:29 AM (21 years ago)

Author:

eugene

Message:

removing pointless Fread/char

File:

• trunk/Ohana/src/opihi/dvo/cmatch.c (modified) (3 diffs)

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

@@ -1 +1 @@
 # include "dvo1.h"
-# define BYTES_STAR 23
-# define BLOCK 1000
-# define DNSTARS 1000
 
 int cmatch (int argc, char **argv) {
@@ -15 +12 @@
   Vector *rvec, *dvec, *mvec, *drvec, *ddvec, *dmvec;
 
+  VarConfig ("GSCDIR", "%s", gscdir);
   VarConfig ("CATDIR", "%s", catdir);
-  VarConfig ("GSCDIR", "%s", gscdir);
+  if (VarConfig ("CATMODE", "%s", catalog.catmode) == NULL) 
+    strcpy (catalog.catmode, "RAW");
       
   if (argc != 9) {
     fprintf (stderr, "USAGE: cmatch file radius (RA) (DEC) (Mag) (dRA) (dDEC) (dMag)\n");
+    fprintf (stderr, "       match a set of object coordinates with a DVO db table\n");
     return (FALSE);
   }
+
+  /*** this function is not well-defined.  re-assess it and re-code it ***/
+  fprintf (stderr, "disabled for now\n");
+  return (FALSE);
 
   radius = atof (argv[2]);
@@ -34 +38 @@
   /* load data from the photometry database file */
   sprintf (filename, "%s/%s", catdir, argv[1]);
-  if (!fits_read_header (filename, &catalog1.header)) {
-    fprintf (stderr, "no stars in %s\n", filename);
+  catalog1.filename = filename;
+  switch (lock_catalog (&catalog, LCK_SOFT)) {
+  case 2:
+    unlock_catalog (&catalog);
+  case 0:
     return (FALSE);
   }
-  f = fopen (filename, "r");
-  if (f == (FILE *) NULL) {
-    fprintf (stderr, "no stars in %s\n", filename);
+  if (!load_catalog (&catalog, LOAD_AVES | LOAD_MEAS | LOAD_SECF, TRUE)) {
+    unlock_catalog (&catalog);
     return (FALSE);
   }
-  fseek (f, catalog1.header.size, SEEK_SET); 
-  
-  /** find number of stars, measurements **/
-  catalog1.Naverage = 0;
-  fits_scan (&catalog1.header, "NSTARS", "%d", 1, &catalog1.Naverage);
-  if (catalog1.Naverage == 0) {
-    fprintf (stderr, "no stars in catalog %s\n", filename);
-    fclose (f);
-    return (TRUE);
-  }
-  ALLOCATE (catalog1.average, Average, catalog1.Naverage);
-  Nitems = catalog1.Naverage;
-  nitems = Fread (catalog1.average, sizeof(Average), Nitems, f, "average");
-  if (nitems != Nitems) {
-    fprintf (stderr, "ERROR: failed to read data from catalog file %s (1)\n", filename);
-    fclose (f);
-    free (catalog1.average);
-    return (FALSE);
-  }
-  fclose (f);
+  unlock_catalog (&catalog);
   fprintf (stderr, "read %d stars from phot catalog file %s\n", catalog1.Naverage, filename);
 
-  /* load data from the GSC file */
-  sprintf (filename, "%s/%s", gscdir, argv[1]);
-  f = fopen (filename, "r");
-  if (f == (FILE *) NULL) {
-    fprintf (stderr, "no stars in %s\n", filename);
-    free (catalog1.average);
-    return (FALSE);
-    /* return (FALSE); */
-  }
-  
   nstar = 0;
   NSTARS = DNSTARS;