cosmo stack hdr

function ds_stacked = cosmo_stack(ds_cell, varargin)