Changeset 41874

Timestamp:

Nov 3, 2021, 11:48:39 AM (5 years ago)

Author:

eugene

Message:

clean up CERROR test, add comments re: CERSTD

File:

• branches/eam_branches/ipp-dev-20210817/Ohana/src/addstar/src/ReadImageHeader.c (modified) (2 diffs)

branches/eam_branches/ipp-dev-20210817/Ohana/src/addstar/src/ReadImageHeader.c

@@ -70 +70 @@
   }
 
+  // 2021.11.03 : consider replacing or supplementing CERROR with CERSTD
+  // if (!gfits_scan (header, "CERROR",   "%lf", 1, &tmp)) tmp = 1.0;
+
   /* require Nastro > 0 unless or ACCEPT_ASTROM */
-  Nastro = 0;
-  gfits_scan (header, "NASTRO", "%d", 1, &Nastro);
-  Cerror = 0;
+  FWHM_X = FWHM_Y = Nastro = Cerror = 0;
+
+  gfits_scan (header, "NASTRO", "%d",  1, &Nastro);
   gfits_scan (header, "CERROR", "%lf", 1, &Cerror);
-  FWHM_X = FWHM_Y = 0;
   gfits_scan (header, "FWHM_X", "%lf", 1, &FWHM_X);
   gfits_scan (header, "FWHM_Y", "%lf", 1, &FWHM_Y);
+
   if (((Nastro == 0) || (Cerror > MAX_CERROR)) && !ACCEPT_ASTROM) {
     fprintf (stderr, "bad astrometric solution in header\n");
@@ -111 +114 @@
   }
 
-  /* CERROR in data file is in arcsec */
-  if (!gfits_scan (header, "CERROR",   "%lf", 1, &tmp)) tmp = 1.0;
-  image[0].cerror = tmp * 50.0;
+  /* CERROR / CERSTD in data file in arcsec; db field is units of 20 mas */
+  image[0].cerror = 50.0 * Cerror;
 
   /* get photcode from header */