Changeset 25636

Timestamp:

Sep 28, 2009, 5:48:38 PM (17 years ago)

Author:

eugene

Message:

get the 180 valid range right

File:

• branches/eam_branches/20090715/Ohana/src/opihi/dvo/skycoverage.c (modified) (6 diffs)

branches/eam_branches/20090715/Ohana/src/opihi/dvo/skycoverage.c

@@ -16 +16 @@
   Coords coords;
   int typehash;
+  int PhotcodeSelect;
+  PhotCode *PhotcodeValue;
 
   WITH_MOSAIC = FALSE;
@@ -43 +45 @@
     remove_argument (N, &argc, argv);
     ByName = TRUE;
+  }
+
+  PhotcodeValue = NULL;
+  PhotcodeSelect = FALSE;
+  if ((N = get_argument (argc, argv, "-photcode"))) {
+    if (!InitPhotcodes ()) return (FALSE);
+    PhotcodeSelect = TRUE;
+    remove_argument (N, &argc, argv);
+    PhotcodeValue = GetPhotcodebyName (argv[N]);
+    if (PhotcodeValue == NULL) {
+      gprint (GP_ERR, "photcode not found in photcode table\n");
+      return (FALSE);
+    }
+    remove_argument (N, &argc, argv);
   }
 
@@ -85 +101 @@
  
   if (argc != 4) {
-    gprint (GP_ERR, "USAGE: skycoverage (buffer) (pixscale) (Npts) [-time start range] [-name name]\n");
+    gprint (GP_ERR, "USAGE: skycoverage (buffer) (pixscale) (Npts)\n");
+    gprint (GP_ERR, "  options: [-time start range] [-trange start stop] [-name name] [-photcode name] [+mosaic] [-mosaic] [-ra-center RA]\n");
     gprint (GP_ERR, "       (buffer) saves bitmapped AIT plot\n");
     gprint (GP_ERR, "       (pixscale) specifies the pixel size in degrees\n");
@@ -128 +145 @@
     for (xs = 0; xs < Nx; xs++) {
       status = XY_to_RD (&r, &d, (double)(xs), (double)(ys), &coords);
-      status &= (r > 0);
-      status &= (r < 360);
+      status &= (r >= 0);
+      status &= (r <= 360);
       if (status) {
         V[ys*Nx + xs] = 2;
@@ -141 +158 @@
     if (ByName && strcmp (name, image[i].name)) continue;
     if (TimeSelect && ((image[i].tzero < tzero) || (image[i].tzero+image[i].trate*image[i].NY > tzero + trange))) continue;
+
+    if (PhotcodeSelect) {
+      if (PhotcodeValue[0].type == PHOT_DEP) {
+        if (PhotcodeValue[0].code != image[i].photcode) continue;
+      } else {
+        if (PhotcodeValue[0].code != GetPhotcodeEquivCodebyCode (image[i].photcode)) continue;
+      }
+    }
+
     if (!FindMosaicForImage (image, Nimage, i)) continue;
 
@@ -167 +193 @@
         XY_to_RD (&r, &d, Xi, Yi, &image[i].coords);
         r = ohana_normalize_angle (r);
+        if (r - RaCenter > +180.0) r -= 360.0;
+        if (r - RaCenter < -180.0) r += 360.0;
         status = RD_to_XY (&Xs, &Ys, r, d, &coords);
         if (Xs < 0) continue;