I found this page on the Japanese povray webring.
http://www.cs-w.com/~earth/material/isopattern.htm
There is here the use of the cells pattern in a way that does not give
diffuse=0 on the edges of step functions. I will study the source code,
but the reason why it looks so good hasn't hit me yet. Any clues?