upd

teuben · Oct 31, 2023 · b1fa6f9 · b1fa6f9
1 parent f891ef4
commit b1fa6f9
Showing 1 changed file with 5 additions and 8 deletions.
diff --git a/src/image/fits/sdinfo.c b/src/image/fits/sdinfo.c
@@ -6,7 +6,7 @@
  * 11-dec-2019   PJT       mdarray reduction example
  * 29-feb-2020   PJT       verbose, raw I/O
  * 28-sep-2021   PJT       better cfitsio usage, report more SDFITS properties
- *    apr-2023   PJT       some mods for the DYSH work
+ *    apr-2023   PJT       some mods/clarifications for the DYSH work
  *
  * Benchmark 6 N2347 files:  2.4"  (this is with mom=0 stats)
  * dims=5 for NGC5291:       31-35ms (depending in 1 or 3 levels)
@@ -41,7 +41,7 @@ string defv[] = {
     "bench=1\n           How many times to run benchmark",
     "mode=-1\n           Mode how much to process the SDFITS files (-1 means all)",
     "datadim=2\n         1: DATA[nrows*nchan]   2: DATA[nrows][nchan]",
-    "VERSION=0.9g\n      19-apr-2023 PJT",
+    "VERSION=0.9h\n      2-may-2023 PJT",
     NULL,
 };
 
@@ -440,6 +440,7 @@ void nemo_main(void)
       dprintf(1,"MODE=3\n");
 
       if (Qraw) ndims = -1;
+
       if (ndims == 0) {  
 
 	for (k=0; k<ncolcheck; ++k) {
@@ -487,15 +488,14 @@ void nemo_main(void)
 	  else
 	    dprintf(0,"DATA2 %g %g ... %g (only 1 row)\n",data2[0][0], data2[0][1], data2[0][nchan-1]);
 	  if (blfit >= 0) {
-	    dprintf(0,"BASELINE %d\n",blfit);
+	    dprintf(0,"BASELINE %d (full row)\n",blfit);
 	    for (ii=0; ii<nrows; ii++)
 	      baseline(nchan, NULL, data2[ii], blfit, coeffs, errors);
 	  }
 	}
 
 	if (mode >= 0 && mode < 5) continue;
 	dprintf(1,"MODE=5\n");
-
 
 	if (mom >= 0) {
 	  dprintf(0,"Stats\n");
@@ -536,8 +536,6 @@ void nemo_main(void)
 		   mean_moment(&m), -1.0,
 		   min_moment(&m), max_moment(&m), n_moment(&m));
 	} // stats
-
-
 
       } else if (ndims > 0) {     // special processing
 
@@ -658,8 +656,7 @@ void nemo_main(void)
       } // ndims > 0
 
       fits_close_file(fptr, &status);            /* close the file */
-      fits_report_error(stderr, status);     /* print out any error messages */
-
+      fits_report_error(stderr, status);         /* print out any error messages */
 
     } //for(j) loop over files
 }