This paper proposes an efficient splitting algorithm to solve coupled liquid gas Darcy and free gas flows at the interface between a porous medium and a free-flow domain. This model is compared to the reduced model introduced in [6] using a 1D approximation of the gas free flow. For that purpose, the gas molar fraction diffusive flux at the interface in the free-flow domain is approximated by a two point flux approximation based on a low-frequency diagonal approximation of a Steklov-Poincaré type operator. The splitting algorithm and the reduced model are applied in particular to the modelling of the mass exchanges at the interface between the storage and the ventilation galleries in radioactive waste deposits.