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