[yt-users] slice of a part of domain

[Nathan Goldbaum]

Hi Chang-Goo,

yt's philosophy for presenting data is to do it exactly as it is
represented on disk. Since your data only fills up a fraction of the
domain, that is how you see it.

That said, I definitely see why you would want to see the duplicated
dataset covering the full domain.  To do that, I think you'll need to
modify the pixelizer or modify the pixelized representation produced by the
pixelizer.  You can take a look at how the pixelizers are set up in
spherical coordinates in yt/geometry/spherical_coordinates.py, near line
76.  Depending on the dimension you are slicing along, yt uses an
orthographic pixelizer or a pixelizer optimized for cylindrical coordinates.

You will have to somehow detect that your data only fills half the domain
and then mirror it before passing it to the pixelizer.  Alternatively you
could manipulate the pixelization its self before reutning it from the
`pixelize` function..

I think this would be a generally useful addition to yt its self.  If you
do go down this path, we would very much like to see a pull request with
these modifications so that others will be able to generate the pixelations
you want.

If you don't want to modify yt, you'll need to do the appropriate data
mirroring on the `data` array in your attached script.

Hope that helps,

Nathan


On Thu, Jul 10, 2014 at 11:10 AM, Chang-Goo Kim <cgkim at astro.princeton.edu>
wrote:

> Hi all,
>
> I have recently run a part of sphere by assuming mirror symmetry. So, the
> simulated domain only spans from 0 to L/2 in every direction, which should
> be -L/2 ~ L/2 for full domain. I then tried to make frb using slice and
> to_frb function. What I get is a slice in -L/4 ~L/4 domain with
> periodicity. How can I get correct slice? Here's a relevant part of script.
>
>         for iax in range(3):
>                 slc = pf.slice(iax,0.0,center=(0.0,0.0,0.0))
>                 slc_frb = slc.to_frb(pf.domain_width[0],512)
>                 frb.append(slc_frb)
>
>                 for i,field in enumerate(fields):
>                         ax = grid[iax*nrow+i]
>                         data=np.array(slc_frb[field])
>                         im = ax.imshow(data, origin='lower',
> norm=LogNorm())
>                         im.set_extent(np.array(slc_frb.bounds))
>
> Thanks,
>
> Chang-Goo
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140710/05e46476/attachment.htm>