Currently staggered grid uses a single index of the "other" variable from default implementation of valmaps. We should enrich this to use interpolations - lowest bar is centered average to achieve 2nd-order accuracy for non-homogeneous equations.
Specialize valmaps on discretization, map index to output of user-specified interpolation scheme.