cosmo dataset slice sa hdr

function ds = cosmo_dataset_slice_sa(ds, to_select)