Changeset 35600 for branches/eam_branches/ipp-20130509

Timestamp:

May 27, 2013, 2:32:48 PM (13 years ago)

Author:

eugene

Message:

kapa name cannot have spaces (treated as separate arguments); add some warning messages in some failure cases; set tcoords.Npolyterms to 0 for image region matches

Location:

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src

Files:

• ImageOps.c (modified) (1 diff)
• load_images.c (modified) (1 diff)
• plot_scatter.c (modified) (2 diffs)
• plotstuff.c (modified) (1 diff)
• select_images.c (modified) (1 diff)

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src/ImageOps.c

@@ -691 +691 @@
   if (!PLOTSTUFF) return;
 
+  if (Npts == 0) {
+    fprintf (stderr, "cannot set image Mcal values yet (nan Mave)\n");
+    return;
+  }
+
   ALLOCATE (xlist, double, Npts);
 

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src/load_images.c

@@ -95 +95 @@
 # endif
 
+  if (Nimage == 0) {
+    fprintf (stderr, "no images selected for analysis : problem with region or photcode?\n");
+    exit (4);
+  }
+
   // match chips to mosaics (if applicable)
   initMosaics (subset, Nsubset, image, inSubset, Nimage);

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src/plot_scatter.c

@@ -48 +48 @@
 
                 Mrel = catalog[i].secfilt[Nsecfilt*j+Nsec].M;
+                if (isnan(Mrel)) continue;
+
                 xlist[N] = Mrel;
                 ylist[N] = PhotSysTiny (&catalog[i].measureT[m], &catalog[i].averageT[j], &catalog[i].secfilt[j*Nsecfilt]) - Mcal - Mmos - Mgrid - Mrel;
@@ -54 +56 @@
             }
         }
+    }
+
+    if (N == 0) {
+      fprintf (stderr, "no valid average values yet\n");
+      continue;
     }
 

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src/plotstuff.c

@@ -23 +23 @@
   char name[100];
 
-  sprintf (name, "gastro [%d]", N);
+  sprintf (name, "gastro:%d", N);
 
   // if -plot is supplied, the plots are shown on the screen; 

branches/eam_branches/ipp-20130509/Ohana/src/relphot/src/select_images.c

@@ -407 +407 @@
   tcoords.pc1_1  = tcoords.pc2_2 = 1.0;
   tcoords.pc1_2  = tcoords.pc2_1 = 0.0;
+
+  tcoords.Npolyterms = 0;
+  memset (tcoords.polyterms, 0, 14*sizeof(float));
   strcpy (tcoords.ctype, "RA---TAN");