Hello, I am trying to make a fixed resolution buffer of a projection of a region in one of my simulations. Here is the relevant part of the script: size=0.025 center=[0.48,0.545,1.003] right0=[center[0]+size,center[1]+size,center[2]+size] left0=[center[0]-size,center[1]-size,center[2]-size] region=pf.h.region(center,left0,right0) p=pf.h.proj(0,"Density",data_source=region) test=FixedResolutionBuffer(p,[0.52,0.57,0.978,1.028],(512,512)) f=h5py.File("test.h5","w") f.create_dataset("/density",data=test["Density"]) f.close When I write the buffer to a h5 file, it seems to be projecting through the whole domain, instead of only the 0.05x0.05x0.05 region I specify. Any ideas what I am doing wrong? Thanks, Mike
Hi Mike, Are you on 2.x or 3.0? 2.x has the nasty, nasty habit of accepting kwargs without validating them, *and* in 2.x, it's the "source" kwarg not "data_source". In 3.0 we've fixed both of these things. -Matt On Wed, Feb 26, 2014 at 5:29 AM, Mike Butler <mbutler.astro@gmail.com> wrote:
Hello,
I am trying to make a fixed resolution buffer of a projection of a region in one of my simulations. Here is the relevant part of the script:
size=0.025 center=[0.48,0.545,1.003] right0=[center[0]+size,center[1]+size,center[2]+size] left0=[center[0]-size,center[1]-size,center[2]-size] region=pf.h.region(center,left0,right0) p=pf.h.proj(0,"Density",data_source=region) test=FixedResolutionBuffer(p,[0.52,0.57,0.978,1.028],(512,512)) f=h5py.File("test.h5","w") f.create_dataset("/density",data=test["Density"]) f.close
When I write the buffer to a h5 file, it seems to be projecting through the whole domain, instead of only the 0.05x0.05x0.05 region I specify. Any ideas what I am doing wrong?
Thanks,
Mike
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Matt, I was using 2.x. If I try to run the same code in 3.0, I get a NotImplementedError at the "57 f.create_dataset("/density",data=test["Density"])" line. However, if I just switch data_source to source, it works fine in 2.x, so I guess that will work for now. Thanks, Mike On Wed, Feb 26, 2014 at 2:36 PM, Matthew Turk <matthewturk@gmail.com> wrote:
Hi Mike,
Are you on 2.x or 3.0? 2.x has the nasty, nasty habit of accepting kwargs without validating them, *and* in 2.x, it's the "source" kwarg not "data_source". In 3.0 we've fixed both of these things.
-Matt
Hello,
I am trying to make a fixed resolution buffer of a projection of a region in one of my simulations. Here is the relevant part of the script:
size=0.025 center=[0.48,0.545,1.003] right0=[center[0]+size,center[1]+size,center[2]+size] left0=[center[0]-size,center[1]-size,center[2]-size] region=pf.h.region(center,left0,right0) p=pf.h.proj(0,"Density",data_source=region) test=FixedResolutionBuffer(p,[0.52,0.57,0.978,1.028],(512,512)) f=h5py.File("test.h5","w") f.create_dataset("/density",data=test["Density"]) f.close
When I write the buffer to a h5 file, it seems to be projecting through
On Wed, Feb 26, 2014 at 5:29 AM, Mike Butler <mbutler.astro@gmail.com> wrote: the
whole domain, instead of only the 0.05x0.05x0.05 region I specify. Any ideas what I am doing wrong?
Thanks,
Mike
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
On Wed, Feb 26, 2014 at 7:47 AM, Mike Butler <mbutler.astro@gmail.com> wrote:
Hi Matt,
I was using 2.x. If I try to run the same code in 3.0, I get a NotImplementedError at the "57 f.create_dataset("/density",data=test["Density"])" line. However, if I just switch data_source to source, it works fine in 2.x, so I guess that will work for now.
Hi Mike, That's pretty weird. Without knowing what wasn't implemented it's hard to say, but that construction otherwise looks good to me. Can you file a bug with the full traceback? -Matt
Thanks,
Mike
On Wed, Feb 26, 2014 at 2:36 PM, Matthew Turk <matthewturk@gmail.com> wrote:
Hi Mike,
Are you on 2.x or 3.0? 2.x has the nasty, nasty habit of accepting kwargs without validating them, *and* in 2.x, it's the "source" kwarg not "data_source". In 3.0 we've fixed both of these things.
-Matt
On Wed, Feb 26, 2014 at 5:29 AM, Mike Butler <mbutler.astro@gmail.com> wrote:
Hello,
I am trying to make a fixed resolution buffer of a projection of a region in one of my simulations. Here is the relevant part of the script:
size=0.025 center=[0.48,0.545,1.003] right0=[center[0]+size,center[1]+size,center[2]+size] left0=[center[0]-size,center[1]-size,center[2]-size] region=pf.h.region(center,left0,right0) p=pf.h.proj(0,"Density",data_source=region) test=FixedResolutionBuffer(p,[0.52,0.57,0.978,1.028],(512,512)) f=h5py.File("test.h5","w") f.create_dataset("/density",data=test["Density"]) f.close
When I write the buffer to a h5 file, it seems to be projecting through the whole domain, instead of only the 0.05x0.05x0.05 region I specify. Any ideas what I am doing wrong?
Thanks,
Mike
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Matthew Turk -
Mike Butler