Added utility function to retrive a submatrix of U that contains a certain species values in each voxel for all timepoints.