cosmo phase stat hdr

function stat_ds = cosmo_phase_stat(ds, varargin)