A numerical model of granular material at different packing fractions and under steady shear is submitted to a sudden shear reversal. We monitor consequences of the strong density and shear rate spatiotemporal heterogeneities, on the constitutive relations. We show that the dynamics can be decomposed into two subsequent regimes spanning a time scale inversely proportional to the shear rate. In the first regime, a nonlocal constitutive relation is satisfied, hence accounting for the spatial heterogeneity of the fluidity parameter. However at later time, we find that the local μ(I) constitutive relation can be applied, in spite of the fact that the fluidity parameter remains heterogeneous.