cosmo dataset slice fa hdr

function ds = cosmo_dataset_slice_fa(ds, to_select)