We use the nonlocal linear hydrodynamic constitutive model, proposed by Evans and Morriss [Statistical Mechanics of Nonequilibrium Liquids (Academic, London, 1990)], for computing an effective spatially dependent shear viscosity of inhomogeneous nonequilibrium fluids. The model is applied to a simple atomic fluid undergoing planar Poiseuille flow in a confined channel of several atomic diameters width. We compare the spatially dependent viscosity with a local generalization of Newton's law of viscosity and the Navier-Stokes viscosity, both of which are known to suffer extreme inaccuracies for highly inhomogeneous systems. The nonlocal constitutive model calculates effective position dependent viscosities that are free from the notorious singularities experienced by applying the commonly used local constitutive model. It is simple, general, and has widespread applicability in nanofluidics where experimental measurement of position dependent transport coefficients is currently inaccessible. In principle the method can be used to predict approximate flow profiles of any arbitrary inhomogeneous system. We demonstrate this by predicting the flow profile for a simple fluid undergoing planar Couette flow in a confined channel of several atomic diameters width. (C) 2004 American Institute of Physics.