Multiscale Gradient Computation for Subsurface Flow Models