Periodicity and cutting planes
Hi YT Users, I’m trying to figure out if there’s a way to get cutting planes to interact with periodicity correctly. I have an ORION data set, and I want to look at a cutting plane in a region that happens to be near a domain boundary. The calculation is periodic, and yt correctly flags it as such. However, when I take a cutting plane that goes off the edge of the computational domain, it doesn’t wrap back around on the other side as it should. Instead, the cells off the domain edge get filled with NaN’s. I have been unable to locate any information in the documentation about whether there is a way to have the cutting plane recognize the periodicity and wrap correctly. Is there a way to do this? Thanks. -- Mark Krumholz _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Mark, I think all of this will happen inside the CPixelize routine. The cutting plane coordinates (px, py, pz) are normalized internal to the pixelization routine to be -0.5*DomainWidth .. 0.5*DomainWidth, with 0,0,0 being the center. So the periodicity check would be off the edges there. I am not sure I can commit to adding this at the moment, but I think it's feasible. The way we do it in the orthogonal pixelizer is to set up a two-step system where we check if we run off the edge, and then check the periodicity loops. -Matt On Tue, Apr 14, 2015 at 12:26 PM, Mark Krumholz <mkrumhol@ucsc.edu> wrote:
Hi YT Users,
I’m trying to figure out if there’s a way to get cutting planes to interact with periodicity correctly. I have an ORION data set, and I want to look at a cutting plane in a region that happens to be near a domain boundary. The calculation is periodic, and yt correctly flags it as such. However, when I take a cutting plane that goes off the edge of the computational domain, it doesn’t wrap back around on the other side as it should. Instead, the cells off the domain edge get filled with NaN’s. I have been unable to locate any information in the documentation about whether there is a way to have the cutting plane recognize the periodicity and wrap correctly. Is there a way to do this? Thanks.
-- Mark Krumholz
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Mark Krumholz
-
Matthew Turk