[yt-users] find densest point in sphere?

[Michael Kuhlen]

Hi

How would I go about finding the densest point in a region, say a sphere?

I'd like something like AMRHierarchy.find_max, but for AMRSphereBase:

pf = load(filename)
sphere = pf.h.sphere(center, radius)
v,c = sphere.find_max("Density")


I guess I could retrieve all high resolution grids with select_grids()
and then do a find_max() on each of those grids, but I'm hoping there's
something nicer.

Thanks,

Mike

-- 
*********************************************************************
*                                                                   *
*  Dr. Michael Kuhlen              Theoretical Astrophysics Center  *
*  email: mqk at astro.berkeley.edu   UC Berkeley                      *
*  cell phone: (831) 588-1468      601 Campbell Hall                *
*  skype username: mikekuhlen      Berkeley, CA 94720               *
*                                                                   *
*********************************************************************

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 262 bytes
Desc: OpenPGP digital signature
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20101206/8103088c/attachment.pgp>