from hydrogen import common

def define():
    l = common('tmp/t.1.node', 3)
    return l

def funV( ts, coor, region, ig, mode = None ):
    from numpy import sqrt

    out = {}
    C = 0.5
    r = sqrt( coor[:,0]**2 + coor[:,1]**2 + coor[:,2]**2 )
    V = - C * 1.0 / r
    out['V'] = V
    return out

@@ -55,24 +55,12 @@
 mon(
-'tmp/t.1.node', 
+dim=
 3)%0A 

use tmp/mesh.vtk in hydrogen