MC59: Sort a sparse matrix to an ordering by columns

This subroutine performs an in-place sort on a sparse matrix to an ordering by columns. There is an option for ordering the entries within each column by increasing row indices and an option for checking for indices that are out-of-range or duplicated.