[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