[yt-users] creating non-square frb
Slavin, Jonathan
jslavin at cfa.harvard.edu
Tue Mar 8 11:30:29 PST 2016
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 at cfa.harvard.edu 60 Garden Street, MS 83
phone: (617) 496-7981 Cambridge, MA 02138-1516
cell: (781) 363-0035 USA
________________________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160308/a805a816/attachment.htm>
More information about the yt-users
mailing list