Index: branches/pap/ppImage/src/ppImageCheckCTE.c
===================================================================
--- branches/pap/ppImage/src/ppImageCheckCTE.c	(revision 23948)
+++ branches/pap/ppImage/src/ppImageCheckCTE.c	(revision 25027)
@@ -60,8 +60,6 @@
             fineRegion = psImageBinningSetFineRegion (binning, ruffRegion);
             fineRegion = psRegionForImage (image, fineRegion);
-            if (fineRegion.x0 >= image->numCols || fineRegion.x1 >= image->numCols ||
-                fineRegion.y0 >= image->numRows || fineRegion.y1 >= image->numRows) {
-                continue;
-            }
+            if (fineRegion.x0 >= image->numCols) continue;
+	    if (fineRegion.y0 >= image->numRows) continue;
 
             psImage *subset  = psImageSubset (image, fineRegion);
@@ -85,8 +83,13 @@
 	    float normVariance = PS_SQR(stats->robustStdev) / cellMedian;
 	    // float normVariance = PS_SQR(stats->sampleStdev) / cellMedian;
+	    // if (!isfinite(normVariance)) fprintf (stderr, "** normVariance is nan **\n");
 
 	    // apply resulting value to the input pixels
 	    for (int jy = fineRegion.y0; jy < fineRegion.y1; jy++) {
+	      if (jy < 0) continue;
+	      if (jy >= image->numRows) continue;
 	      for (int jx = fineRegion.x0; jx < fineRegion.x1; jx++) {
+		if (jx < 0) continue;
+		if (jx >= image->numCols) continue;
 		image->data.F32[jy][jx] = normVariance;
 	      }