cosmo phase itc hdr

function itc_ds = cosmo_phase_itc(ds, varargin)