Index: /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c
===================================================================
--- /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c	(revision 16585)
+++ /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.c	(revision 16586)
@@ -4,6 +4,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.79.2.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-02-22 02:00:30 $
+ *  @version $Revision: 1.79.2.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-02-22 02:27:40 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -636,5 +636,5 @@
 bool pmSubtractionConvolve(pmReadout *out1, pmReadout *out2, const pmReadout *ro1, const pmReadout *ro2,
                            const psImage *subMask, psMaskType blank, const psRegion *region,
-                           const pmSubtractionKernels *kernels, bool useFFT)
+                           const pmSubtractionKernels *kernels, bool doBG, bool useFFT)
 {
     PS_ASSERT_PTR_NON_NULL(out1, false);
@@ -802,6 +802,6 @@
             // (with the stamps) that it does not vary rapidly on this scale.
             polyValues = p_pmSubtractionPolynomial(polyValues, kernels->spatialOrder, xNorm, yNorm);
-            float background = p_pmSubtractionSolutionBackground(kernels, polyValues); // Background term
-            background = 0.0;
+            float background = doBG ? p_pmSubtractionSolutionBackground(kernels, polyValues) :
+                0.0; // Background term
             psRegion subRegion = psRegionSet(i, xSubMax, j, ySubMax); // Sub-region to convolve
 
Index: /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.h
===================================================================
--- /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.h	(revision 16585)
+++ /branches/pap_branch_080214/psModules/src/imcombine/pmSubtraction.h	(revision 16586)
@@ -6,6 +6,6 @@
  * @author GLG, MHPCC
  *
- * @version $Revision: 1.22 $ $Name: not supported by cvs2svn $
- * @date $Date: 2007-12-07 01:57:15 $
+ * @version $Revision: 1.22.6.1 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2008-02-22 02:27:40 $
  * Copyright 2004-207 Institute for Astronomy, University of Hawaii
  */
@@ -92,4 +92,5 @@
                            const psRegion *region, ///< Region to convolve (or NULL)
                            const pmSubtractionKernels *kernels, ///< Kernel parameters
+                           bool doBG,   ///< Apply background term?
                            bool useFFT  ///< Use Fast Fourier Transform for the convolution?
     );
Index: /branches/pap_branch_080214/psModules/src/imcombine/pmSubtractionMatch.c
===================================================================
--- /branches/pap_branch_080214/psModules/src/imcombine/pmSubtractionMatch.c	(revision 16585)
+++ /branches/pap_branch_080214/psModules/src/imcombine/pmSubtractionMatch.c	(revision 16586)
@@ -463,5 +463,6 @@
 
             psTrace("psModules.imcombine", 2, "Convolving...\n");
-            if (!pmSubtractionConvolve(conv1, conv2, ro1, ro2, subMask, maskBlank, region, kernels, useFFT)) {
+            if (!pmSubtractionConvolve(conv1, conv2, ro1, ro2, subMask, maskBlank, region, kernels,
+                                       true, useFFT)) {
                 psError(PS_ERR_UNKNOWN, false, "Unable to convolve image.");
                 goto MATCH_ERROR;
