def evalF(x,y,z): return dfgeom.evalField(shape,x,y,z)
def evalF(x,y,z): #return -1.0 return dfgeom.evalField(shape,x,y,z)