add a XD recipe to combine wavelength solutions from 2 different sources, eg sky emission for some orders and lamp for the others.

Author: KathleenLabrie

geminidr/gnirs/recipes/sq/recipes_ARC_XD_SPECT.py

@@ -26,5 +26,20 @@ def makeProcessedArc(p):
     p.determineDistortion(spatial_order=1, step=4)
     p.storeProcessedArc()
 
+def combineWavelengthSolutions(p):
+    """
+    Combine wavelength solutions from two different sources, such as orders
+    from sky lines and the others from the lamp.
+
+    Inputs are:
+      * First processed arc
+      * Second processed arc
+      * List of extensions to take from the second processed arc
+    """
+    p.rejectInputs(at_start=1, outstream="second_arc")
+    p.rejectInputs(at_end=1)
+    p.combineSlices(from_stream="second_arc")
+    p.storeProcessedArc(suffix='_combinedArc')
+    p.writeOutputs()
 
 _default = makeProcessedArc