cosmo distatis hdr

function res = cosmo_distatis(ds, varargin)