In yt-2.x, there were boolean data objects that could do this, but they have not been fully implemented in yt 3:

http://yt-project.org/docs/dev/analyzing/objects.html#combining-objects-boolean-data-objects

You could always use yt 2.x to do this for now...

Cameron

On Tue, Jun 16, 2015 at 9:18 PM, Stephanie Tonnesen <stonnes@gmail.com> wrote:

Hi yt-users!

I am trying to analyze a cylindrical loop in my dataset.  So what I want to do is to create two disks:

d1 = ds.disk([0.5, 0.5, 0.5], [0., 0., 1.], (8, 'kpc'), (1, 'kpc'))

d2 = ds.disk([0.5, 0.5, 0.5], [0., 0., 1.], (7, 'kpc'), (1, 'kpc'))

and have r1 be the region containing all the cells in d1 and NOT in d2.  I feel like there must be an easy way to do this, but can't think of it.  Also, searching the yt-users archive kept coming up with an error.

Thanks in advance!

Stephanie

--
Dr. Stephanie Tonnesen

Alvin E. Nashman Postdoctoral Fellow

Carnegie Observatories, Pasadena, CA

stonnes@gmail.com

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

--

Cameron Hummels

Postdoctoral Researcher

Steward Observatory

University of Arizona

http://chummels.org

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