Changes to help run the RS method for a 2D master equation.

The new version iterates over J when choosing the active-space matrix
bandwidth. Several subscript errors are also fixed (e.g. indexing by one
subscript when two are required).