Boundary shape, particularly roughness, strongly controls the amount of wall slip in dense granular flows. In this paper, we aim to quantify and understand which aspects of a dense granular flow are controlled by the boundary condition, and to incorporate these observations into a cooperative nonlocal model characterizing slow granular flows. To examine the influence of boundary properties, we perform experiments on a quasi-2D annular shear cell with a rotating inner wall and a fixed outer wall; the later is selected from among 6 walls with various roughness, local concavity, and compliance. We find that we can successfully capture the full flow profile using a single set of empirically determined model parameters, with only the wall slip velocity set by direct observation. Through the use of photoelastic particles, we observe how the internal stresses fluctuate more for rougher boundaries, corresponding to lower wall slip, and connect this observation to the propagation of nonlocal effects originating at the wall.