cosmo map2surface hdr

function s = cosmo_map2surface(ds, fn, varargin)