function [ds_useful, msk] = cosmo_remove_useless_data(ds, dim, type)
cosmo randperm hdr
cosmo run tests hdr