new_contents: 
  # find outliers based on DVARS and FD
  
  def outlier(file_name, save_name):
      # load data
      data = np.loadtxt(file_name)
      outlier = []
      # set outlier values to 0
      for i in data:
          if i <= 0.5:
              outlier.append(0)
          else:
              outlier.append(i)
      # find outlier indices
      outlier_indices = np.nonzero(outlier)
      # write file to store outlier indices
      f = open(save_name, 'w')
      f.write(str(outlier_indices))
      f.close()
  
subject: add function to find outliers based on dvars and fd measurements