Index: trunk/ppImage/src/ppImageCheckNoise.c
===================================================================
--- trunk/ppImage/src/ppImageCheckNoise.c	(revision 33590)
+++ trunk/ppImage/src/ppImageCheckNoise.c	(revision 33667)
@@ -77,5 +77,5 @@
             *stats = *statsDefaults;
             statsDefaults->tmpData = NULL;
-
+	    
             psImageStats (stats, subset, submask, 0xffff);
 
@@ -84,5 +84,23 @@
 	    // float normVariance = PS_SQR(stats->sampleStdev) / cellMedian;
 	    // if (!isfinite(normVariance)) fprintf (stderr, "** normVariance is nan **\n");
-
+	    psTrace("ppImage.noise",3,
+		    "PPIM.N: %s %s %d %d (%d %d) %f\n",
+		    psMetadataLookupStr(NULL,inReadout->parent->parent->concepts,"CHIP.NAME"),
+		    psMetadataLookupStr(NULL,inReadout->parent->concepts,"CELL.NAME"),
+		    ix,iy,subset->numCols,subset->numRows,normVariance);
+		    
+	    for (int dy = 0; dy < subset->numRows; dy++) {
+	      for (int dx = 0; dx < subset->numCols; dx++) {
+		psTrace("ppImage.noise",5,
+			"PPIN: %s %s %d %d %d %d %f %d %f\n",
+			psMetadataLookupStr(NULL,inReadout->parent->parent->concepts,"CHIP.NAME"),
+			psMetadataLookupStr(NULL,inReadout->parent->concepts,"CELL.NAME"),
+			ix,iy,dx,dy,
+			subset->data.F32[dy][dx],
+			(int) submask->data.PS_TYPE_VECTOR_MASK_DATA[dy][dx],
+			normVariance);
+	      }
+	    }
+	    //	    binning->data.F32[iy][ix] = normVariance;
 	    // apply resulting value to the input pixels
 	    for (int jy = fineRegion.y0; jy < fineRegion.y1; jy++) {
@@ -100,5 +118,5 @@
         }
     }
-
+    
     psFree (rng);
     psFree (binning);