In this paper we extend nonlocal smoothing techniques for image regularization in [12] to surface regularization, with surfaces represented by level set functions. We test our algorithm on both phantom and observed surfaces, including city terrain and cortical surfaces.