ENH: Add t2star map to t2star workflow outputs

Author: mgxd

nibabies/workflows/bold/t2s.py

@@ -72,7 +72,7 @@ def init_bold_t2s_wf(echo_times, mem_gb, omp_nthreads, name="bold_t2s_wf"):
 
     inputnode = pe.Node(niu.IdentityInterface(fields=["bold_file", "bold_mask"]), name="inputnode")
 
-    outputnode = pe.Node(niu.IdentityInterface(fields=["bold"]), name="outputnode")
+    outputnode = pe.Node(niu.IdentityInterface(fields=["bold", "t2star_map"]), name="outputnode")
 
     LOGGER.log(25, "Generating T2* map and optimally combined ME-EPI time series.")
 
@@ -82,7 +82,8 @@ def init_bold_t2s_wf(echo_times, mem_gb, omp_nthreads, name="bold_t2s_wf"):
     workflow.connect([
         (inputnode, t2smap_node, [('bold_file', 'in_files'),
                                   ('bold_mask', 'mask_file')]),
-        (t2smap_node, outputnode, [('optimal_comb', 'bold')]),
+        (t2smap_node, outputnode, [('optimal_comb', 'bold'),
+                                   ('t2star_map', 't2star_map')]),
     ])
     # fmt: on