Add swirl example and necessary library changes

Adds aux array capabilities to grid_eval functions along with a Riemann
solver that averages speeds from the aux arrays.
Clean up advection demo to use FieldIndexers