cosmo surficial neighborhood hdr

function [nbrhood, vo, fo, out2in] = cosmo_surficial_neighborhood(ds, surfs, varargin)