cosmo stat hdr

function stat_ds = cosmo_stat(ds, stat_name, output_stat_name)