I suspect the answer will be field cuts, britton can elaborate...

On Feb 3, 2009, at 11:53 AM, Stephen Skory wrote:

Hi all,

I'd like to make some parallel projections for a sub-region. What is the appropriate way to do that? I've tried
something like this, below, but it makes a projection for the whole
box, not the sub-region. I've tried a few other ways, something like
this <http://paste.enzotools.org/show/38/>, but it works in only serial but not in parallel, where it crashes, like this <http://paste.enzotools.org/show/39/> (that's abbreviated output).

i = 1
pc = PlotCollection(pf.h.region([0.5]*3,[i*.1,0.,0.],[i*.1 + .1,1.,1.]),center=[0.5]*3)
pc.add_projection("Density", 0)
pc.save("slice-0.%d" % i)


sskory@physics.ucsd.edu           o__  Stephen Skory
http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student

yt-users mailing list

Dr. Eric J. Hallman
NSF Astronomy and Astrophysics Postdoctoral Fellow                    
Center for Astrophysics and Space Astronomy 
University of Colorado at Boulder  
hallman (at) casa.colorado.edu 
office(s): (303) 735-0129 / (303) 492-7484