Self-consistent three-dimensional steady state solutions of the MHD equations with field-aligned incompressible flow