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 #############

def readin_temp(temp_file):
    """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 #############
    """
    with open(temp_file, 'r') as fid:
        temp = np.loadtxt(fid, skiprows=1, usecols=[2])

    return temp