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.


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.


yt-users mailing list