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
I don't have a test dataset that looks like this. Any chance you can share one of your outputs (possibly off-list)? If so I can try to make an example over here.
-Nathan
On Tue, Mar 8, 2016 at 1:30 PM, Slavin, Jonathan jslavin@cfa.harvard.edu wrote:
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 ________________________________________________________
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org