Here we present a dataset of Receiver Functions from the South Line of the MEDUSA seismological array used in the studies of Pearce et al. (2012) and Millet et al. (2019).The geographical extend of the dataset is from (37.2N, 21.7E) to (39.0N, 24.6E).The data were aligned using multi-channel cross-correlation on the main P-wave arrival and obtained through frequency domain water-level type deconvolution.The data is three component (radial [R], transverse [T], and vertical [Z]) and the frequency band of the data is 0.03 to 0.3 Hz.It is written in a 4D matrix in netCDF format, with the first two dimensions corresponding to the source-receiver pairs grid, and the third and fourth to the time series data from the corresponding three component Receiver Function for all pairs.For source-receiver pairs where no recording is available the data is a string of zeroes.Metadata include information about the stations and events locations and timings, as well as recording parameters for the data.