cosmo chunkize hdr

function chunks = cosmo_chunkize(ds, nchunks_out)