cosmo montecarlo phase stat hdr

function ds_stat = cosmo_montecarlo_phase_stat(ds, varargin)