Allow specification of which field to reflect for wall BCs.

For some systems, and depending on how the q fields are ordered, it may be useful to reflect different fields at "wall" boundaries.  This makes it possible to do so by setting e.g.

    solver.reflect_index = [2,3,1]