Hi Dave,

You can use the field_cuts keyword for just such a thing. You can supply multiple strings to be eval'ed in a list. For example,

proj = pf.h.proj(0, 'Density', field_cuts=['grid["Density"] > 0.1', 'grid["Density"] < 0.2'])

Of course, you could also use one long string connected with ands. This may not be documented, but certainly should be.

Britton

On Mon, Apr 22, 2013 at 6:00 PM, David Collins <dcollins4096@gmail.com> wrote:

Hi, all--

Can I hand a cut_region to a projection?

I want to make a bunch of projections based on parameters, and don't want to have to instantiate 47 different fields. I'm hoping something like

region = pf.h.all_data().cut_region('grid["Density"] > 0.1 and grid["Density"] < 0.2')
proj = pf.h.proj(0,"Density",region=cut_region)

or something like that. Is that possible? I'm poking around the help and source, but no method immediately popped out at me.

Thanks,
d.

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org