function [nbrhood, vo, fo, out2in] = cosmo_surficial_neighborhood(ds, surfs, varargin)
cosmo surface dataset hdr
cosmo synthetic dataset hdr