Changeset 8643

Timestamp:

Aug 27, 2006, 1:29:25 PM (20 years ago)

Author:

eugene

Message:

sidereal time from Keyword or from jd and longitude

File:

• trunk/Ohana/src/addstar/src/gstars.c (modified) (3 diffs)

trunk/Ohana/src/addstar/src/gstars.c

@@ -46 +46 @@
   }
     
-  /* get ST (used for airmass calculation) */
-  gfits_scan (&header, STKeyword, "%s", 1, line);
-  /* remove ':' characters */
-  for (c = strchr (line, ':'); c != (char *) NULL; c = strchr (line, ':')) { *c = ' '; }
-  sscanf (line, "%d %d %lf", &hour, &min, &sec);
-  image[0].sidtime = hour + min/60.0 + sec/3600.0;
-  image[0].latitude = Latitude;
-
   /* CERROR in data file is in arcsec */
   if (!gfits_scan (&header, "CERROR",   "%lf", 1, &tmp)) tmp = 1.0;
@@ -111 +103 @@
     image[0].tzero = parse_time (&header);
   }
+
+  if (STKeyword[0]) {
+    /* get ST (used for airmass calculation) */
+    gfits_scan (&header, STKeyword, "%s", 1, line);
+    /* remove ':' characters */
+    for (c = strchr (line, ':'); c != (char *) NULL; c = strchr (line, ':')) { *c = ' '; }
+    sscanf (line, "%d %d %lf", &hour, &min, &sec);
+    image[0].sidtime = hour + min/60.0 + sec/3600.0;
+  } else {
+    double jd;
+    jd = sec_to_jd (image[0].tzero);
+    image[0].sidtime  = ohana_lst (jd, Longitude);
+  }
+  image[0].latitude = Latitude;
 
   tmp = 0;
@@ -213 +219 @@
 
     XY_to_RD (&stars[N].R, &stars[N].D, stars[N].X, stars[N].Y, &image[0].coords);
-    while (stars[N].R < 0.0) stars[N].R += 360.0;
+    while (stars[N].R <    0.0) stars[N].R += 360.0;
     while (stars[N].R >= 360.0) stars[N].R -= 360.0;
     stars[N].found = -1;