Re: [yt-dev] problem with disk data container and related fields in yt-3.0

Hi Matt, Britton, Given a periodic box, a center at the bottom edge of the box, and a height that is nearly half the box height, I would expect to see some cells appearing at the top of your plot. That may have been a change from 2.X (since YTDisk accounts for periodicity), but I'm not sure it's incorrect behavior. There is definitely something strange going on that eats away that portion of the lower part of the cylinder, which must be something to do with the way the code computes the periodic radial distance, then subtracts the portion projected onto the cylinder normal to obtain the cylindrical radius. Douglas Rudd Scientific Computing Consultant Research Computing Center drudd@uchicago.edu On Aug 12, 2014, at 7:32 AM, Matthew Turk wrote:

Hi Britton,

On Tue, Aug 12, 2014 at 7:30 AM, Britton Smith wrote:

...

Hi all,

Some issues seem to have cropped up in the disk data container in yt-3.0. I made a disk object in a periodic simulation in the following way: http://paste.yt-project.org/show/5021/

I then made a project of the "cylindrical_z" field, which yielded the following: http://i.imgur.com/0XSv5Zs.png

For reference, the domain is ~45 Mpc in length on a side, so the sphere should be just the object at the bottom, extending a little less than halfway up. I don't know what that top part is about. My guess is that it's an issue with periodicity.

Yes, looks like it to me.

...

Secondly, I made a similar projection of the "cylindrical_r" field, which gave this: http://i.imgur.com/h8SEJp8.png

As you can see, it is not symmetric as one would expect. For reference, this is what I get from yt-2.x: http://i.imgur.com/O1xT9MI.png This is what I would expect.

Are you doing a MIP projection here? Could you do a MIP projection so we can compare?

It almost looks to me like the center is set incorrectly.

...

Does anyone have any insight on this?

Britton

_______________________________________________ yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org

_______________________________________________ yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org