Upper panel: phase folded radial velocity measurements of WASP-39 obtained combining data from SOPHIE (filled-circles) and CORALIE (open-squares) spectrographs. Superimposed is the best-fit model RV curve with parameters from Table 4. The centre-of-mass velocity for each data set was subtracted from the RVs (γ_SOPHIE =  −58.4826 km s^-1 and γ_CORALIE =  −58.4708 km s^-1). Lower panel: we show the bisector span measurements as a function of radial velocity, values are shifted to a zero-mean (⟨ V_span ⟩ _SOPHIE =  −0.032 km s^-1,  ⟨ V_span ⟩ _CORALIE =  −0.051 km s^-1). The bisector span shows no significant variation nor correlation with the RVs, suggesting that the signal is mainly due to Doppler shifts of the stellar lines rather than stellar profile variations due to stellar activity or a blended eclipsing binary.