code: def readin_temp(temp_file):
      
      with open(temp_file, 'r') as fid:
          temp = np.loadtxt(fid, skiprows=1, usecols=[2])
  
      return temp
docs: The temperature file should be in mag-format: header + 3 columns with
      coordinates and value of temperature. The coordinates  have to be the same
      as from the resistivity data.
      Such a temperature file can be produced with #############