  We use numerical simulations to study the flow of athermal, frictionless,
soft-core two dimensional spherocylinders driven by a uniform steady-state
simple shear applied at a fixed volume and a fixed finite strain rate
$\dot\gamma$. Energy dissipation is via a viscous drag with respect to a
uniformly sheared host fluid, giving a simple model for flow in a non-Brownian
suspension with Newtonian rheology. We study the resulting spatial structure of
the sheared system, and compute correlation functions of the velocity, the
particle density, the nematic order parameter, and the particle angular
velocity. Correlations of density, nematic order, and angular velocity are
shown to be short ranged both below and above jamming. We compare a system of
size-bidisperse particles with a system of size-monodisperse particles, and
argue how differences in spatial order as the packing increases leads to
differences in the global nematic order parameter. We consider the effect of
shearing on initially well ordered configurations, and show that in many cases
the shearing acts to destroy the order, leading to the same steady-state
ensemble as found when starting from random initial configurations.
