def METHOD_NAME(frequencies, path=os.path.dirname(os.path.realpath(__file__))):
    """
    Read out cable induced loss and phase. From standard 4 channel station.
    """
    cable_discrete = np.loadtxt(os.path.join(path, 'HardwareResponses/CableAntennuation_James2016.csv'), skiprows=1,
                                delimiter=',')
    max_frequency = 5000. * units.MHz
    if np.max(frequencies) > max_frequency:
        max_frequency = np.max(frequencies)
    # Convert to GHz
    cable_discrete[:, 0] *= units.Hz
    cable_discrete[0, 0] = 0.
    cable_discrete[-1, 0] = max_frequency
    cable_amp_db_f = interp1d(cable_discrete[:, 0], cable_discrete[:, 1])
    cable_amp_db = cable_amp_db_f(frequencies)
    cable_amp = 10 ** (cable_amp_db / 20.)
    cable_phase_f = interp1d(cable_discrete[:, 0], np.unwrap(np.deg2rad(cable_discrete[:, 2])))
    cable_phase = cable_phase_f(frequencies)
    cable_phase = np.exp(1j * cable_phase)
    return cable_amp * cable_phase