A buoyancy-induced stationary flow with viscous dissipation in a horizontal porous layer is investigated. The lower boundary surface is impermeable and subject to a uniform heat flux. The upper open boundary has a prescribed, linearly varying, temperature distribution. The buoyancy-induced basic velocity profile is parallel and non-uniform. The linear stability of this basic solution is analysed numerically by solving the disturbance equations for oblique rolls arbitrarily oriented with respect to the basic velocity field. The onset conditions of thermal instability are governed by the Rayleigh number associated with the prescribed wall heat flux at the lower boundary, by the horizontal Rayleigh number associated with the imposed temperature gradient on the upper open boundary, and by the Gebhart number associated with the effect of viscous dissipation. The critical value of the Rayleigh number for the onset of the thermal instability is evaluated as a function of the horizontal Rayleigh number and of the Gebhart number. It is shown that the longitudinal rolls, having axis parallel to the basic velocity, are the most unstable in all the cases examined. Moreover, the imposed horizontal temperature gradient tends to stabilise the basic flow, while the viscous dissipation turns out to have a destabilising effect.