Hi all, I've been creating square fixed resolution buffer objects to use for plotting using matplotlib and it's been working fine. Now I did a calculation that has rectangular output, i.e. 2D cylindrical going from r = 0 - 30 pc and z = -30 - +30 pc and for some reason I can't seem to create an frb to cover the whole computational domain: ds = yt.load(file) slc = ds.slice(2,0.) frb = slc.to_frb((width=(30, 'pc'), height=(60,'pc'),resolution=(512,1024)) In [91]: frb.limits['r'] Out[91]: (-1.63714434867e+14 code_length, 9.25701637144e+19 code_length) In [92]: frb.limits['z'] Out[92]: (-4.62851637144e+19 code_length, 4.62851637144e+19 code_length) So the frb is centered correctly, but only includes half the domain in the z direction. In [94]: ds.domain_left_edge Out[94]: YTArray([ 0.00000000e+00, -9.25700000e+19, 0.00000000e+00]) code_length In [96]: ds.domain_right_edge Out[96]: YTArray([ 9.25700000e+19, 9.25700000e+19, 6.28318531e+00]) code_length Can anyone tell me a way around this? Thanks, Jon -- ________________________________________________________ Jonathan D. Slavin Harvard-Smithsonian CfA jslavin@cfa.harvard.edu 60 Garden Street, MS 83 phone: (617) 496-7981 Cambridge, MA 02138-1516 cell: (781) 363-0035 USA ________________________________________________________