get_function_grads is deprecated.
Saving shock indicator to avg.plt
Better checks for positivity limiter.