import ceo
gmt = ceo.GMT_MX()
src = ceo.Source('V',rays_box_sampling=101,rays_box_size=25.5,rays_origin=[0,0,25])
src>>(gmt,)
+src
print(src.wavefront.rms(-9))