Ok, I solved the problems thanks to Mesh iterators. Still don't know what cells() does though. Alessio