cosmo index unique hdr

function [cell_indices, unique_values] = cosmo_index_unique(values)