Index: trunk/psModules/src/imcombine/pmSubtraction.c
===================================================================
--- trunk/psModules/src/imcombine/pmSubtraction.c	(revision 33089)
+++ trunk/psModules/src/imcombine/pmSubtraction.c	(revision 33425)
@@ -796,5 +796,5 @@
     pmSubtractionStamp *stamp = job->args->data[0]; // List of stamps
     pmSubtractionKernels *kernels = job->args->data[1]; // Kernels
-    int footprint = PS_SCALAR_VALUE(job->args->data[2], S32); // Stamp index
+    int footprint = PS_SCALAR_VALUE(job->args->data[2], S32); // Stamp index -- MEH - it is?
 
     return pmSubtractionConvolveStamp(stamp, kernels, footprint);
@@ -832,9 +832,11 @@
 
 #ifdef TESTING
+    //MEH - index conflict or changed in past?
     for (int j = 0; j < kernels->num; j++) {
         if (stamp->convolutions1) {
             psString convName = NULL;
-            psStringAppend(&convName, "conv1_%03d_%03d.fits", index, j);
-            psFits *fits = psFitsOpen(convName, "w");
+            //psStringAppend(&convName, "conv1_%03d_%03d.fits", index, j);
+            psStringAppend(&convName, "conv1_xxx_%03d.fits", j);
+	    psFits *fits = psFitsOpen(convName, "w");
             psFree(convName);
             psKernel *conv = stamp->convolutions1->data[j];
@@ -845,6 +847,7 @@
         if (stamp->convolutions2) {
             psString convName = NULL;
-            psStringAppend(&convName, "conv2_%03d_%03d.fits", index, j);
-            psFits *fits = psFitsOpen(convName, "w");
+            //psStringAppend(&convName, "conv2_%03d_%03d.fits", index, j);
+            psStringAppend(&convName, "conv2_xxx_%03d.fits", j);
+	    psFits *fits = psFitsOpen(convName, "w");
             psFree(convName);
             psKernel *conv = stamp->convolutions2->data[j];