[yt-users] Creating a ring

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