
Feb. 26, 2014
11:29 a.m.
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