This paper presents an alternative boundary conditions setup for the numerical simulations of stably stratified flow. The focus of the tested computational setup is on the pressure boundary conditions on the artificial boundaries of the computational domain. The simple three dimensional test case deals with the steady flow of an incompressible, variable density fluid over a low smooth model hill. The Boussinesq approximation model is solved by an in-house developed high-resolution numerical code, based on compact finite-difference discretization in space and Strong Stability Preserving Runge-Kutta method for (pseudo-) time stepping.