@pipeBase.timeMethod

-    def run(self, dataRef, selectDataList=[]):

+    def run(self, dataRef, selectDataList=[], tempExpRefList=None):

         """!

         \brief Assemble a coadd from a set of Warps

@@ -327,13 +327,14 @@ discussed in \ref pipeTasks_multiBand (but note that normally, one would use the

                  - nImage: exposure count image

         """

         skyInfo = self.getSkyInfo(dataRef)

-        calExpRefList = self.selectExposures(dataRef, skyInfo, selectDataList=selectDataList)

-        if len(calExpRefList) == 0:

-            self.log.warn("No exposures to coadd")

-            return

-        self.log.info("Coadding %d exposures", len(calExpRefList))

-

-        tempExpRefList = self.getTempExpRefList(dataRef, calExpRefList)

+        if tempExpRefList is None:

+            calExpRefList = self.selectExposures(dataRef, skyInfo, selectDataList=selectDataList)

+            if len(calExpRefList) == 0:

+                self.log.warn("No exposures to coadd")

+                return

+            self.log.info("Coadding %d exposures", len(calExpRefList))

+

+            tempExpRefList = self.getTempExpRefList(dataRef, calExpRefList)