function [ds, params] = cosmo_normalize(ds, params, dim)
cosmo nfold partitioner hdr
cosmo norminv hdr