The objective of this study is to present a numerical treatment of the pressure gradient when control volumes are sharing the interface between a homogeneous fluid and a porous medium. Two possible approaches, e.g. linear interpolation and extrapolation, are considered, and they are applied to the case of a steady and two-dimensional curved channel flow which is partially filled with a porous medium. It was found that the linear extrapolation produces a continuous velocity-field at the interface and thus is recommended. On the contrary, the linear interpolation entails a discontinuous velocity field at the interface, thereby warning its use in connection with the Brinkman-Forchheimer-extended Darcy flow model.