[Amr] paraview-yt halo finding

[Jorge Poco]

Hi Matt

I will try to use HaloFinder from YT, the script would be:

# create Paraview Handler
pf = stream.StreamStaticOutput(handler)

# halos will be added to this polydata
spheres = vtkAppendPolyData()

halos = HaloFinder(pf)
for halo in halos:
    # create a vtkSphere using halo information
    s = vtkSphereSource()
    s.SetCenter(halo.get_sphere().center)
    s.SetRadius(halo.get_sphere().radius)
    s.Update()
    spheres.AddInput(s.GetOutput())
spheres.Update()


output.SetPoints(s1.GetOutput().GetPoints())

output.SetPolys(s1.GetOutput().GetPolys())



Using this script we will be able to recover the halos as a set of spheres
which

could be visualized using paraview.


Questions:

* HaloFinder needs the particle information right?

* Current stream.StreamHandler support the particle data?


Regards,



Att. Jorge Poco
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/amr/attachments/20110617/1660d5fd/attachment-0001.html>