Index: /branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtraction.c
===================================================================
--- /branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtraction.c	(revision 30343)
+++ /branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtraction.c	(revision 30344)
@@ -25,4 +25,6 @@
 #include "pmSubtractionVisual.h"
 #include "pmSubtractionThreads.h"
+
+bool psFitsWriteImageSimple (char *filename, psImage *image, psMetadata *header);
 
 #include "pmSubtraction.h"
@@ -1437,4 +1439,7 @@
     psFree(kernelErr2);
 
+    static int nOut1 = 0;
+    static int nOut2 = 0;
+
     // Calculate covariances
     // This can be fairly involved, so we only do it for a small number of instances
@@ -1450,4 +1455,17 @@
                 psKernelTruncate(kernel, covarFrac);
                 covars->data[i] = psImageCovarianceCalculate(kernel, ro1->covariance);
+		if (0) {
+		    char name[128];
+		    snprintf (name, 128, "covar.sample1.%03d.fits", nOut1);
+		    psKernel *cov = covars->data[i];
+		    psFitsWriteImageSimple (name, cov->image, NULL);
+
+		    snprintf (name, 128, "incovar.sample1.%03d.fits", nOut1);
+		    psFitsWriteImageSimple (name, ro1->covariance->image, NULL);
+
+		    snprintf (name, 128, "conv.sample1.%03d.fits", nOut1);
+		    psFitsWriteImageSimple (name, kernel->image, NULL);
+		    nOut1 ++;
+		}
                 psFree(kernel);
             }
@@ -1471,4 +1489,17 @@
                 psKernelTruncate(kernel, covarFrac);
                 covars->data[i] = psImageCovarianceCalculate(kernel, ro2->covariance);
+		if (0) {
+		    char name[128];
+		    snprintf (name, 128, "covar.sample2.%03d.fits", nOut2);
+		    psKernel *cov = covars->data[i];
+		    psFitsWriteImageSimple (name, cov->image, NULL);
+
+		    snprintf (name, 128, "incovar.sample2.%03d.fits", nOut1);
+		    psFitsWriteImageSimple (name, ro2->covariance->image, NULL);
+
+		    snprintf (name, 128, "conv.sample2.%03d.fits", nOut1);
+		    psFitsWriteImageSimple (name, kernel->image, NULL);
+		    nOut2 ++;
+		}
                 psFree(kernel);
             }