old_contents: 
  '''
  Create final data cubes
  '''
  
  import glob
  import os
  import sys
  
  from taskinit import ia
  
  output_type = sys.argv[-1]
  
  path = "/global/scratch/ekoch/combined/single_channels/"
  out_path = "/global/scratch/ekoch/combined/cube_images/"
  
  if output_type == "image":
      out_name = os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.image")
      search_string = "*.clean.image"
  elif output_type == "residual":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.residual")
      search_string = "*.clean.residual"
  elif output_type == "model":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.model")
      search_string = "*.clean.model"
  elif output_type == "mask":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.mask")
      search_string = "*.clean.mask"
  else:
      raise TypeError("Must provide 'image', 'residual', 'mask' or 'model' as argument.")
  
  channel_direcs = glob.glob(os.path.join(path, "channel_*"))
  
  # Sort by the channel number
  min_chan = 10000
  max_chan = 0
  
  for channel in channel_direcs:
      # Last 3 characters
      chan_num = int(channel.split("_")[-1])
      if chan_num < min_chan:
          min_chan = chan_num
  
      if chan_num > max_chan:
          max_chan = chan_num
  
  ordered_chans = []
  for chan in xrange(min_chan, max_chan+1):
      i = 0
      while True:
          chan_num = int(channel_direcs[i].split("_")[-1])
          if chan_num == chan:
              ordered_chans.append(channel_direcs[i])
              break
          i += 1
      else:
          print("No log found for channel "+str(chan))
  
  channel_direcs = ordered_chans
  
  images = []
  for channel in channel_direcs:
      image = glob.glob(os.path.join(channel, search_string))
      if not image:
          print("Cannot find log file in "+channel)
          continue
  
      images.append(image[0])
  
  ia.imageconcat(outfile=out_name, infiles=images)
  
new_contents: 
  '''
  Create final data cubes
  '''
  
  import glob
  import os
  import sys
  
  from taskinit import ia
  
  output_type = sys.argv[-1]
  
  path = "/global/scratch/ekoch/combined/single_channels/"
  out_path = "/global/scratch/ekoch/combined/cube_images/"
  
  if output_type == "image":
      out_name = os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.image")
      search_string = "*.clean.image"
  elif output_type == "residual":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.residual")
      search_string = "*.clean.residual"
  elif output_type == "model":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.model")
      search_string = "*.clean.model"
  elif output_type == "mask":
      out_name = \
          os.path.join(out_path, "M33_14B-088_AT0206_HI.clean.mask")
      search_string = "*.clean.mask"
  else:
      raise TypeError("Must provide 'image', 'residual', 'mask' or 'model' as argument.")
  
  channel_direcs = glob.glob(os.path.join(path, "channel_*"))
  
  # Sort by the channel number
  min_chan = 10000
  max_chan = 0
  
  for channel in channel_direcs:
      # Last 3 characters
      chan_num = int(channel.split("_")[-1])
      if chan_num < min_chan:
          min_chan = chan_num
  
      if chan_num > max_chan:
          max_chan = chan_num
  
  ordered_chans = []
  for chan in xrange(min_chan, max_chan+1):
      i = 0
      while True:
          chan_num = int(channel_direcs[i].split("_")[-1])
          if chan_num == chan:
              ordered_chans.append(channel_direcs[i])
              break
          i += 1
      else:
          print("No log found for channel "+str(chan))
  
  channel_direcs = ordered_chans
  
  images = []
  for channel in channel_direcs:
      image = glob.glob(os.path.join(channel, search_string))
      if not image:
          print("Cannot find log file in "+channel)
          continue
  
      images.append(image[0])
  
  print images
  
  ia.imageconcat(outfile=out_name, infiles=images)
  
subject: Print out the images