cosmo mask dim intersect hdr

function [indices_cell, ds_intersect_cell] = cosmo_mask_dim_intersect(ds_cell, dim, varargin)