docs: def load_rho(name, column):
      '''Load a datafile with rho structure like mag and phase
      '''
      try:
          content = np.loadtxt(name, skiprows=1, usecols=([column]))
      except:
          raise ValueError('Given column to open does not exist.')
  
      return content